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1.  INTRODUCTION 

The  magnetospheric  substorm  is  a  natural  phenomenon  associated  with  splendid 
visual  effects  and  complex  scientific  descriptions.  The  substorm  concept  came 
from  Akasofu  who  attempted  to  organize  a  wealth  of  all-sky  camera  data  into  a 
single  cohesive  picture.  Today,  Akasofu' s  original  auroral  substorm  concept  is 
known  to  be  manifest  in  many  complementary  data  sets  (magnetic  bay  phenomena, 
polar  magnetic  substorms,  ionospheric  substorms,  proton  aurora,  VLF  emissions,  etc.), 
and  it  has  become  apparent  that  a  full  understanding  of  the  substorm  process  will 
require  a  very  comprehensive  knowledge  of  the  magnetosphere,  including  its 
coupling  mechanisms  to  the  solar  wind  and  ionosphere. 

The  substorm  is  a  major  event  characterized  by  time  dependence,  vast  spatial 
scales,  and  probably  non  linear  physical  processes.  For  these  reasons,  a  complete 
observational  definition  of  a  substorm  is  not  yet  available.  In  this  research, 
we  addressed  certain  aspects  of  substorm  behavior  using  specific  data  sources 
and  well  defined  models. 

In  chapter  2,  we  report  on  a  study  made  of  the  feasibility  of  substorm  de¬ 
tection  at  subauroral  latitudes  that  would  utilize  the  AFGL  magnetometer  network. 

High  latitude  magnetic  disturbances  and  gradients  are  very  large  during  sub¬ 
storms  and  are  often  associated  essentially  with  overhead  auroras  and  currents. 

In  contrast,  mid  latitude  magnetic  signatures  are  smaller  but  distributed  more 
widely  geographically  than  at  high  latitudes.  We  have  studied  the  amplitude  and 
longitudinal  extent  of  mid  latitude  signatures  to  serve  as  indicators  of  sub¬ 
storm  occurrence.  The  reliability  of  this  approach,  evaluated  for  available  data 
samples,  depends  upon  such  factors  as  the  longitudinal  wedge  of  the  mid  latitude 
station  coverage,  the  local  time  of  this  wedge  during  a  substorm  event  and  the 
intensity  and  duration  of  the  substorm.  From  our  study,  a  mid  latitude  magnetic 
index  is  defined  and  an  operational  scheme  developed  for  substorm  detection. 


In  Chapter  3,  a  review  is  given  of  micropulsations  as  indicators  of  geo¬ 
magnetic  activity.  The  AFGL  sub  auroral  chain  has  a  digital  data  base  and  is  well 


suited  to  be  used  in  a  procedure  to  formulate  and  calibrate  a  micropulsation 
index  of  geomagnetic  activity  level. 

Chapters  4  and  5  deal  with  the  global  ionospheric  distribution  of  currents 
driven  by  field  aligned  currents,  and  the  properties  of  the  ground  magnetic 
signature  associated  with  such  currents.  The  model  developed  to  determine  global 
ionospheric  distributions  uses  input  variables  that  are  known  to  calculate 
selected  output  variables  quantitatively.  The  input  variables  are  the  driving 
currents  distributed  in  the  auroral  zone  and  the  ionospheric  conductivity.  These 
input  variables  are  defined  by  observational  evidence.  The  calculated  output 
variables  are  the  ionospheric  electric  field,  potential ,  and  currents.  Model 
limitations  and  desirable  generalizations  are  discussed  in  detail. 

Determination  of  the  ground  magnetic  signature  from  a  configuration  of  currents 
is,  in  principle,  a  straight-forward  matter  if  all  the  currents  are  known.  In 
our  case,  however,  we  have  a  firm  handle  only  on  the  ionospheric  portion  of  the 
currents.  The  crucial  missing  parts  of  the  current  network  are  induced  earth 
currents  and  the  upwards  (three  dimensional)  extension  and  closure  of  the  field 
aligned  currents  in  the  magnetosphere.  The  comparative  roles  of  Pederson  and  Hall 
currents,  the  shielding  effects  of  the  ionosphere,  and  other  physical  considerations 
necessary  for  a  reliable  determination  of  the  ground  magnetic  signature  are  dis¬ 
cussed  In  context . 

Chapters  6  and  7  relate  to  studies  made  on  auroral  boundary  locations.  Ob¬ 
servational  information  is  given  in  Chapter  6  on:  the  location  of  the  equator- 


ward  auroral  boundary  as  a  function  of 


and  AE;  and  the  location  of  the  pole- 


ward  boundary  as  a  function  of  Dst.  Boundary  locations  were  determined  from  DMSP 


images.  In  Chapter  7  a  model  is  developed  that  relates  properties  of  the  ring 


current  to  the  size  of  the  polar  cap.  Through  parametric  representation,  the 
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axially  symmetric  ring  current  evolves  into  a  storm  time  configuration  with  an 
associated  growth  in  Dst.  Ring  current  parameters  are  the  inner  edge  location 
and  the  distance  where  the  beta  of  the  plasma  is  unity.  The  ring  current  is 
shown  to  contribute  substantially  to  the  change  in  polar  cap  size  (taken  to  be 
determined  by  the  poleward  auroral  boundary)  during  magnetic  storms. 


I 

i 


I 


* 


2.  MID  LATITUDE  INDICATORS  OF  SUBSTORM  ACTIVITY 


A.  Introduction 

Magnetospher ic  suhstorms  influence  magnetic  fields  on  the  earth's  surface 
at  places  far  removed  from  the  site  of  spectacular  auroral  displays  and  ion¬ 
ospheric  disruptions.  Magnetic  bays  at  sub  auroral,  mid,  and  low  latitudes 
have  been  attributed  to  ionospheric  rt turn  currents  from  .auroral  electrojets 
(Ros^oker,  1966),  or  to  an  asymmetric  partial  ring  current  closing  in  the 
magnetosphere  (Cummings,  1966;  Akasolu  ami  Mong.  1964;  Crook or.  and  McPherrun 
197?;  Fnkunishi  jtnd  Kam ide ,  1971).  In  a  study  of  coordinated  satellite  and 
ground  based  data,  McPhcrron  (1471)  used  mid  la' itude  posit  tvr  II  havs  near  loeal 
midnight  as  a  prime  indicator  of  suhstorm  onset  .  lie  noted  that  midlatitude 
magnetic  signatures  were  simpler  than  auroral  zone  records  which  are  commonly 
affected  hv  Intense,  rapidlv  changing  overhead  currents.  Other  magnetic 
phenomena  which  have  been  connected  with  suhstorms  are  hursts  of  Pi  2  pulsations 
(Sa i to  et  ai . ,  1976)  and  long  period  oscillat  ions  (Wang  et  al  . ,  1977). 

The  purpose  of  this  study  is  to  explore  the  ut  i  1  i  ?.at  i  on  of  the  AFCL  sub- 
auroral,  five  station  magnetometer  network  to  detect  the  occurrence  and  strength 
of  magnet ospherlc  substorms.  These  AFtd.  stations  operate  automatically  to 
provide  high  time  resolution,  )  component  magnetic  data  (Knecht  and  Paaicli,  1976), 
and  are  located  in  Washington,  South  Dakota,  '.'isconsin,  Michigan,  and 
Massachusetts  (Figure  2.1),  Tile  station  network  spans  about  4  hours  of  local 
time  (60°  of  longitude)  at  66°  CGI.  (Corrected  Geomagnetic  Latitude).  (Two 
additional  stations  in  California  and  Florida  art'  now  operational  hut  are  not 
used  in  this  study.)  Lincoln  (1967)  and  Ro stoker  (1972)  have  reviewed  and 
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appraised  the  more  common  geomagnetic  indices.  Details  of  the  most  familiar 
indices,  and  AE  (Auroral  Electro/jet),  are  given  by  Mayaud  (1967),  Davis  and 
Suguira  (1966),  and  Allen  and  Kroehl  (1975).  An  implicit  requirement  for  the 
quality  of  any  geomagnetic  index  is  reasonably  complete  longitudinal  coverage, 
a  property  that  the  AFGL  network  does  not  offer.  Given  this  limitation,  we 
must  focus  on  the  local  time  dependence  of  magnetic  effects  accompanying  sub¬ 
storms  . 

Contour  maps  of  mid  latitude  magnetic  variations  in  local  time  vs.  universal 
time  were  constructed  by  Clauer  and  McPIierron  (1 974)  to  describe  the  time  de¬ 
velopment  of  substorms  as  seen  by  stations  at  low  latitudes.  They  emphasize  that 
mid  latitude  signatures  are  highly  dependent  on  location  relative  to  the  sub¬ 
storm's  central  meridian.  For  example,  an  idealized  magnetospheric-ionospheric 
current  system  centered  on  local  midnight  produces  positive  H  variations  about 
local  midnight  and  smaller  negative  H  changes  on  the  dayside  tending  to  zero  near 
local  noon.  Caan  et  al .  (1978)  and  Kokubun  and  lijlrna  (1975)  employed  superposed 
epoch  methods  to  depict  the  statistical  local  time  behavior  of  midlatitude  sig¬ 
natures  accompanying  selected  substorms.  The  selection  criteria  in  Caan  et  al . 
were  limited  to  mid  latitude  bays  within  5  hours  of  local  midnight  to  avoid  con¬ 
tamination  by  si  events.  World  wide  magnetic  changes  caused  by  an  si  could  be 
misinterpreted  as  indication  of  a  substorm  if  only  one  local  time  sector  is 
considered.  Substorm  triggering  by  an  si  or  ssc  has  been  analyzed  by  Kokubun 
et  al .  (1977),  who  found  that  high  AE  (>  lOOy)  and  southward  were  necessary 
preconditions  before  the  solar  wind  event.  They  suggested  that  a  favorable  tail 
configuration  allowed  triggering. 

Substorms  are  usually  more  complicated  than  a  simple  process  of  expansion 
and  recovery.  In  their  study  of  multiple  onset  substorms,  Pytte  et  al.  (1976) 
stated  that  Pi  2  bursts  seemed  to  accompany  each  brightening  of  an  equatorvard 


auroral  arc,  while  sub  auroral  and  mid  latitude  magnetic  activity  signaled  west¬ 
ward  travelling  surges.  However,  when  substorm  expansion  begins  during  the 
recovery  of  an  earlier  substorm  in  the  same  local  time  sector,  mid  latitude 
positive  H  bays  near  local  midnight  may  not  clearly  distinguish  the  two  events 
(Kamide  et  al . ,  1977),  Thus,  it  is  not  a  simple  matter  to  identify  unequivocally 
every  substorm  on  the  basis  of  mid  latitude  magnetic  bays,  especially  during 
extended  periods  of  magnetic  conditions  which  favor  continuous  substorm  activity. 

A  simple  current  system  designed  to  explain  high  and  low  latitude  magnetic 
disturbances  during  substorm  expansion  was  discussed  in  a  detailed  analysis  of 
a  single  event  (McPherron  et  al.,1973).  Rostoker  (1974)  suggested  a  modification 
of  this  system  featuring  field  aligned  currents  which  feed  the  east  and  west 
electrojets  at  dawn  and  dusk,  and  has  an  upward  field  aligned  current  at  the  Harang 
discontinuity  connecting  to  the  tail  current.  More  recent  treatments  of  current 
systems  seek  to  explain  conditions  during  moderate  activity  (Hughes  and  Rostoker, 
1977),  the  dynamics  of  dayside  auroral  and  cusp  movements  (Yashuhara  et  al . ,  1975; 
Eather  et  al . ,  1979),  and  seasonal  asymmetries  between  northern  and  southern 
hemispheres  (Cizlcr  et  al . ,  1976).  Finally,  there  is  a  class  of  events  reported 
by  Pytte  et  al .  (1978)  which  they  call  "convection  substorms".  These  events 
produce  sustained  negative  bays  in  the  auroral  zone  but  exhibit  hardly  any  mid 
latitude  indications  of  field  aligned  currents  originating  in  the  tail.  Instead 
of  impulsive  releases  of  stored  energy  symptomatic  of  ordinary  substorms,  it 
appears  that  a  steady  transfer  of  energy  from  the  solar  wind  to  the  ionosphere 
may  power  the  process. 

Park  and  Meng  (1973)  found  that  mid  latitude  magnetic  activity  did  not 
correlate  w<l!  with  ionospheric  electric  fields  deduced  from  F  layer  disturbances, 
and  concluded  that  the  currents  responsible  for  the  mid  latitude  activity  instead 
must  flow  In  the  magnetosphere.  Model  studies  by  Yashuhara  et  al.(1975)  with 
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field  aligned  currents  into  a  conducting  ionosphere  also  showed  remarkably  little 
current  spreading  to  mid  and  low  latitudes. 

Additional  studies  of  the  way  field  aligned  currents  distribute  globally  in 
the  ionosphere  have  been  conducted  quantitatively  by  Nopper  and  Carovl llano  (1978), 
Nisbet  et  al.(1978),  and  Kamide  and  Matsushita  (1979).  These  authors  use  field 
aligned  current  patterns  as  inputs,  like  those  derived  from  satellite  data  by 
Iijima  and  Pot etna  (1978),  and  solve  for  the  ionospheric  electric  field  using  tensor 
conductivity  models  and  Ohm's  law.  (This  topic  is  the  subject  of  chapter  4  of 
this  report.)  One  result  of  this  work  is  that  effects  traditionally  credited  to 
equivalent  ionospheric  currents  at  mid  latitudes,  which  represent  return  currents 
from  the  auroral  electrojets,  are  a  consequence  of  the  field  aligned  current  systems. 

This  brief  introduction  should  serve  to  highlight  some  of  the  complexities 
and  uncertainties  surrounding  the  topic  of  substorm  current  systems.  The  current 
effort  is  phenomenological  in  approach  and  concentrates  on  the  statistical  re¬ 
lationship  between  changes  in  the  midlatitude  H  component  (X  components  in  the 
AFGL  chain)  and  the  intensity  of  the  AE  index.  The  procedure  we  develop  has 
limitations,  but  it  also  has  the  virtues  of  simplicity  and  objectivity.  In 
addition,  our  method  should  be  applicable  even  when  all  data  is  interrupted 
momentarily,  or  when  one  or  more  stations  is  lost  for  an  extended  period. 

B.  The  Relationship  Between  Auroral  and  Mid  Latitude  Magnetic.  Signatures 

In  this  section  the  relationship  between  auroral  zone  and  mid  latitude  magnetic 
activity  is  examined  qualitatively.  The  AFGL  magnetometer  chain  is  our  main 
source  of  information  on  mid  latitude  activity;  magnetograms  from  Fredericksburg 
(FB)  and  Boulder  (BD)  supply  additional  data  in  the  local  time  sector  of  the  AFGL 
chain.  Auroral  zone  activity  indicating  substorms  is  selected  from  H  component 
magnetograms  recorded  at  9  high  latitude  stations  of  the  AE  network.  A  map 
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showing  the  stations  used  in  this  study  appears  in  Figure  2.1.  Table  2.1  lists 
the  station  names,  symbols,  and  Corrected  Geomagnetic  coordinates.  Figure  2.2 
displays  5  cases  of  magnetic  bay  onsets,  designated  by  vertical  lines,  which 
occurred  on  December  9,  1976  in  the  high  latitude  network.  As  expected  from  the 
extensive  literature  on  substorm  studies,  bay  onsets  usually  are  observed  first 
by  stations  located  near  local  midnight.  Thus,  from  00-12  UT,  onsets  generally 
are  seen  in  the  Canadian  sector,  while  from  12-24  UT  it  is  the  Russian  sector 
that  is  affected  the  most.  The  vertical  dashed  lines  in  Figure  2.2  correspond 
to  those  in  Figure  2.3,  which  shows  X  and  Y  components  from  the  AFGL  chain  on 
the  same  day,  December  9,  1976.  First,  we  note  that  in  Figure  2.2  significant 
changes  in  the  average  slope  of  the  X  and  Y  components  at  mid  latitude  coincide 
with  high  latitude  magnetic  bay  on  onsets  (solid  lines).  These  changes  occur  at 
most  or  all  of  the  AFGL  stations  and  are  roughly  simultaneous  with  the  bay  onsets 
(within  about  20  minutes).  On  this  particular  day,  it  appears  that  for  the  event 
at  1520  UT  the  Y  component  gives  a  somewhat  better  indication  of  substorm  onset, 
although  less  distinct  but  perceptible  changes  in  the  X  components  did  take  place 
about  20  minutes  earlier. 

In  addition  to  the  activity  in  Figures  2.3  and  2.4  accompanying  auroral  zone 
bay  onsets,  there  are  gross  patterns  at  other  times  which  also  indicate  some 
degree  of  high  latitude  disturbances.  This  "class"  of  events  typically  is  described 
at  mid  latitude  by  sharper  changes  in  the  Y  component  than  in  the  X  component; 
such  cases  are  marked  by  dashed  vertical  lines  in  both  Figures  2.3  and  2.4. 

Referring  again  to  the  high  latitude  magnetograms  in  Figure  2.2,  we  notice  that 
the  solid  lines  (by  definition)  indicate  onsets  of  major  magnetic  bays,  while  the 
dashed  lines,  selected  on  the  basis  of  mid  latitude  magnetic  activity,  occur  with 
less  pronounced  bays  taking  place  both  before  and  after  major  onsets. 


Northern  hemisphere  map  showing  locations  of  stations  used  in  this 
report.  Circles  defined  by  ticks  are  !0°  intervals  of  Corrected 
Geomai’net ic  Latitude;  tick  marks  at  one  hour  intervals  of  local  time 
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St at  ion 

Symbol 

CG  Latitude  (N) 

CG  Longitude 

Tixic  Bay 

TI 

65.6 

195.1 

Dixon  Island 

DI 

67.9 

154.7 

Ki runa 

KR 

64.3 

104  .8 

Le  irvogur 

LR 

66.6 

71.1 

Great  Whale  R 

GRW 

68.2 

353.7 

Fort  Churchill 

FC 

70.3 

325.9 

Baker  Lake 

BL 

75.1 

320.4 

Col  1 ege 

CO 

64.9 

260.3 

F  reder i ckshurg 

FB 

50 

352 

Boo  1  dor 

RD 

50 

315 

Al’Gi  Chain 

Sudbury 

MA 

•wn 

U"» 

OO 

1.9 

Mt .  Clemens 

MI 

55.8 

344  .8 

Camp  Deng  las 

WI 

56.3 

334.2 

Rapid  City 

SD 

54.1 

317.3 

Newport 

WA 

55.2 

299.6 

(able  I  1. 1st  of  stations  used  in  tills  report.  Symbols  appear  at  station  locations 
in  Figure  2.1.  Corrected  Geomagnetic  Coordinates  are  North  Latitude  and 
Last  longitude.  AFGL  station  symbols  indicate  the  state  in  which  the 
installation  is  located. 
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Figure  2.2  Stacked  H  component  magnetograms  (X  component  at  FC  and  BL)  from  9 
auroral  2ono  stations  on  December  9,  1976.  Bars  Indicate  2D0y  de¬ 
flections;  up  is  positive.  Dots  under  each  trace  mark  local  midnight. 
Solid  and  dashed  vertical  lines  are  malor  and  minor  events  discussed 
in  th»'  text. 


Stacked  mapnc top rams  l  rom  AI-'f 'I.  chain  st-U  ions  (MA,  Ml,  WI ,  SI),  nnd  VA) 
on  Dec.  •) ,  19  7b  (000(1-1 200UT)  .  X  components  abnvr ,  Y  components  below. 
Hots  under  each  (  race  at  local  midnipht.  Solid  vertical  line:' 
correspond  to  events  in  I'lpure  ?.?.  Amplitude  bar  in  50v  for  X 
r'OT;i|i(niL’iil  s ;  for  Y  components,  amplitude  bar  is  SOy  at  MA,  70y  at 
"i.  HO,  at  Wl  ami  WA,  and  1 20-,  at  SI). 
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The  purpose  of  this  discussion  is  to  demonstrate  the  range  of  possible 
inferences  to  be  drawn  from  a  mid  latitude  magnetic  data  set  spanning  4  hours 
of  local  time.  A  reasonable  conclusion  from  the  rather  complicated  example  in 
Figures  2.3  and  2.4  is  this: 

A  sustained,  rapid  change  in  the  X  and  Y  components 
occurring  almost  simultaneously  at  more  than  one 
midlatitude  station  indicates  that  significant 
magnetic  bay  activity  has  commenced  at  auroral 
zone  latitudes;  definite  changes  mainly  in  the 
Y  components  at  mid  latitudes  suggest  less  im¬ 
portant  high  latitude  activity  occurring  before 
or  after  major  onsets. 

Different  interpretations  are  possible  in  the  single  example  cited  in  the  figures. 

For  example,  it  could  be  argued  in  Figure  2.4  that  the  activity  at  1420  UT  is 
seen  also  in  the  X  components  at  MA  and  MI,  or  that  the  major  high  latitude  bay 
beginning  at  1520  UT  is  not  reflected  clearly  in  the  midlatitude  X  component  data. 

In  the  next  section  we  try  to  remove  or  reduce  such  ambiguities  and  uncertainties. 
This  is  done  by  making  detailed  comparisons  of  mid  latitude  data  from  Fredericksburg 
and  Boulder  (used  to  substitute  for  unavailable  AFGL  chain  data  at  the  time  of 
this  study)  with  AE  index  data  representing  auroral  zone  magnetic  activity.  Finally, 
a  simple  quantitative  scheme  using  AFGL  chain  data  is  proposed  and  examples  of  its 
application  given. 

C.  Correlations  between  Mid  Latitude  Activity  and  AE 

To  assess  the  applicability  of  mid  latitude  magnetic  activity  occurring  at 
various  local  times,  we  analyzed  4  months  of  magnetograms  from  Fredericksburg  (FB) 
and  Boulder  (BD).  The  magnetic  activity  at  these  two  stations  were  compared  with 
the  level  of  auroral  zone  disturbance  depicted  in  plots  of  2.5  minute  averaged 
Auroral  Electrojet  (AE)  index  values.  Auroral  zone  events,  defined  here  as 
significant  enhancements  in  AE  discernible  by  eye,  were  collected  into  200y 
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intervals  of  increasing  strength.  The  peak  value  of  AE  determined  both  the  time 
and  strength  of  an  event.  Examples  of  AE  events  occurring  in  an  8  day  period 
(October  9-16,  1972)  appear  in  Figure  2.5.  The  peak  values,  in  hundreds  of 
gammas,  are  given  above  each  enhancement.  H  component  magnetograms  from  FB 
over  the  same  time  period  are  shown  in  Figure  2.6;  AE  enhancements  and  their 
strength  are  indicated  at  appropriate  times  above  each  trace  in  the  same  manner 
as  Figure  2.5.  It  is  obvious  that  there  is  a  wide  range  of  mid  latitude  activity 
accompanying  auroral  zone  events,  but  that,  in  a  rough  way,  activity  at  FB  is 
stronger  with  larger  AE  events. 

In  order  to  uncover  a  more  definite  relationship  between  auroral  zone  and 
mid  latitude  magnetic  behavior,  a  better  characterization  of  mid  latitude  activity 
is  needed.  After  a  survey  of  FB  and  BD  magnetograms  it  was  decided  that  three 
descriptive  levels  of  activity  would  be  adequate;  these  levels  are: 

(1)  Strong:  a  definite,  sudden  change  in  the  average  slope  of  the 
H  or  D  component; 

(2)  Moderate:  distinct  irregular  activity  in  the  range  of  5-20y 
amplitude; 

(3)  Weak:  little  or  no  apparent  activity. 

The  strong  level  corresponds  to  the  mid  latitude  indication  of  auroral  zone  bay 
onset  given  in  Figures  2.3  and  2.4.  The  moderate  level  includes  lesser  activity 
which  is  nonetheless  obvious,  like  that  mentioned  in  Figures  2.3  and  2.4. 

In  our  procedure,  the  activity  level  at  FB  and  BD  is  assigned  by  the  above 
criteria  at  the  time  of  each  AE  event.  This  yields  4  seperate  activity  estimates 
for  an  event,  one  from  each  component  (H  and  D)  at  both  stations.  Then  the  4 
estimates  are  "superposed"  so  that  the  most  active  station  component  characterizes 
each  event.  Such  a  procedure  is  similar  to  creating  AE  indices,  where  the  most 
highly  disturbed  station  contributes  solely  to  the  index.  But,  unlike  the  entire 
circle  of  local  time  coverage  provided  by  the  AE  network,  FB  and  BD  span  less  than 


H  component  mngnetograms  From  Fredericksburg,  October  9  -  16,  1972. 
Events  marked  by  arrows  correspond  to  those  in  Figure  2.5,  with  the 
peak  AK  value  (In  hundreds  of  gammas)  indicated  above.  Solid  triangl 
on  UT  scale  Indicates  local  noon  at  Fredericksburg. 
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3  hours  of  local  time. 

Figure  2.7  shows  how  806  AE  events  occurring  from  October  2,  1972  through 
January,  1973  were  distributed  in  strength  and  time.  The  first  histogram  gives 
the  number  of  events  in  each  200y  interval  of  increasing  AE.  AE  ranges  for  each 
interval  are  listed  in  the  figures.  The  second  histogram  shows  the  diurnal 
occurrence  of  events,  regardless  of  strength,  collected  in  3  hour  intervals  of  UT. 
Corresponding  to  each  UT  interval  there  is  a  local  time  in  Figure  2.7  which  gives 
the  sector  that  FB  and  BD  occupied  at  the  time  of  each  event.  The  midpoint 
between  the  two  stations  defines  the  local  time  sector. 

The  distribution  of  events  by  strength  shows  that  almost  two-thirds  of  all 
events  are  less  than  600y  in  amplitude.  Above  600y  the  numbers  decrease  smoothly 
and  rapidly  to  about  25  cases  (37.)  above  1 000y .  The  diurnal  variation  is  less 
regular;  occurrences  peak  markedly  during  18-21  UT  and  decrease  by  about  45%  at 
06-09  UT. 

All  806  AF.  events  are  categorized  by  the  three  level  (qualitative)  description 
of  mid  latitude  activity  defined  earlier.  Using  clusters  of  3-bar  histograms. 

Figures  2.8  and  2.9  display  actual  numbers  of  events  falling  into  anv  combination 
of  AE  strength,  local  time  sector,  and  mill  latitude  activity.  The  first,  second, 
and  third  bars  in  each  3-bar  histogram  cluster  enumerate  cases  of  weak  (W) , 
moderate  (M) ,  and  strong  (S)  midlatitude  activity,  respectively.  The  3-bar  histo¬ 
grams  are  grouped  by  intervals  of  increasing  AK  strength  into  8  histograms  dis¬ 
tinguished  by  the  local  time  sector  in  which  mid  latitude  activity  was  recorded. 

This  somewhat  complicated  format  gives  the  actual  numbers  used  below  in  determining 
conditional  probabilities;  in  addition,  it  reveals  general  trends  and  permits  an 
appraisal  of  statistical  significance  which  would  be  lost  if  ratios  were  taken  • 

immediately . 

Considering  one  example  in  detail,  the  upper  left  histogram  in  Figure  2.8 
includes  all  AE  events  that  occurred  when  FB  and  BD  were  in  the  07-10  CGET 
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Figure  2.7  Distribution  of  806  AE  events,  identified  during  October  1972  -  January 
197  3,  according  to  strength  (left  histogram)  in  200y  intervals  of  AF., 
and  Ul  of  occurrence  (right  histogram),  in  8  three-hour  intervals  of 
UT  and  CGLT  at  FB  and  BD.  Midpoint  between  FB  and  BD  defines  the 
CGI.T  interval.  Intervals  span  the  beginning  of  the  first  hour  to  the 
beginning  of  the  second  (00  -  03  UT  =  0000  -  0300  Ul ) . 
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.rt  breakdown  ol  mid  latitude  responses  In  AK  invents  oi'r  irrinj.’  In  October 
1  ')72  -  danuarv  197  1.  l-aclt  l- bar  cluster  I'.ivus  tin-  number  of  weak  (W) , 
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(Corrected  Geomagnetic  Local  Time)  sector.  Weak  mid  latitude  activity  dominated 
for  events  in  the  two  lowest  strength  intervals  (AE  <  400y).  Strong  mid  latitude 
activity  prevailed  when  the  AE  of  an  event  was  greater  than  400y.  For  all  day- 
side  sectors  (Figure  2.8)  the  majority  condition  at  mid  latitudes  changes  from 
weak  to  strong  activity  with  AE  =  400-600y.  When  FB  and  BD  were  on  the  night- 
side  (Figure  2.9)  the  same  change  to  strong  mid  latitude  activity  occurs  with 
weaker  substorm  events  (AE  *  200-400y). 

Percentages  for  mid  latitude  magnetic  activity  calculated  from  Figures 
2.8  and  2.9  are  presented  in  Figure  2.10.  Again,  each  histogram  in  Figure  2.10 
represents  all  events  that  occurred  while  FB  and  BD  were  in  the  indicated  local 
time  sector.  Percentages  are  plotted  vs.  AE  intervals  with  the  same  ranges  as 
in  Figures  2. 7-2.9.  The  set  labeled  strong  shows  relative  percentages  of  strong 
activity,  compared  with  moderate  and  weak,  for  increasing  levels  of  AE  strength. 
The  set  labeled  "strong  and  moderate"  combines  events  with  strong  or  moderate 
mid  latitude  activity. 

Figure  2.10  depicts  the  data  in  Figures  2.8  and  2.9  compactly  and  quantifies 
the  properties  of  the  three  level  actlvitv  index  described  above.  For  example, 
approximately  50%  of  the  mid  latitude  responses  in  10-13  CGLT  sector  are  strong 
for  an  AE  enhancement  of  400-600y  (interval  3).  The  percentage  increases  to 
nearly  70%  in  that  local  time  sector  for  the  same  AE  category  when  moderate  re¬ 
sponses  are  added  (lower  set  of  histograms).  A  careful  comparison  of  the  upper 
two  rows  of  histograms  reveals  that  strong  responses  reach  or  exceed  50%  (48%  at 
10-H  CGLT)  when  AE  is  400-600y  or  more  for  davt ime  local  sectors,  plotted  in  the 
top  row.  In  the  second  row  containing  nighttime  sectors,  r'0%  is  reached  when 
AF.  is  only  200-400y  (interval  2). 

In  general,  the  occurrence  of  strong  activity  at  mid  latitudes  suggests 
different  probabilities  of  high  latitude  activity,  depending  on  whether  the  ob¬ 
servation  is  made  on  the  dav  or  night  side  of  the  earth.  On  the  day  side. 
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discrimination  between  relatively  weak  and  strong  AE  activity  is  quite  good.  That 
is,  strong  responses  are  rather  infrequent  until  AE  exceeds  400y.  On  the  night 
side  discrimination  is  poorer  because  weaker  AE  enhancements  elicit  a  large  share 
of  strong  responses,  but  the  reliability  is  greater  because  almost  all  very  large 
AE  enhancements  are  detected;  all  responses  are  strong  for  AE  greater  than  600y 
in  the  19-04  CGLT  sector,  and  nearly  100%  in  the  04-07  CGLT  sector. 

There  are  two  exceptions  to  the  general  day  and  night  behavior.  First, 
in  the  lowest  AE  interval  of  the  04-07  CGLT  sector,  strong  responses  occur  pro¬ 
portionately  less  compared  with  the  other  nighttime  sectors.  Second,  and  more 
interestingly,  the  noon  sector  (10-13  CGLT)  behaves  more  like  the  night  sectors 
in  that  strong  responses  exceed  90%  during  large  AE  enhancements  (>  600y) .  This 
reverses  the  trend  beginning  in  the  dawn  sector  (04-07  CGLT)  of  decreased  sensitivity 
to  large  AE  events.  After  noon  (13-19  CGLT)  the  proportion  of  strong  responses 
decreases  again  until  a  sharp  recovery  occurs  after  dusk  (19-22  CGLT). 

In  summary,  two  general  statements  can  be  made  concerning  the  effectiveness 
of  strong  mid  latitude  magnetic  response,  defined  earlier  in  this  section,  as  an 
indicator  of  high  latitude  magnetic  activity  characterized  by  the  AE  index. 

(1)  At  stations  located  in  the  night  sectors,  strong 
responses  accompany  over  50%  of  AE  enhancements  in  the 
200-400y  range.  With  AE  >  600y  the  percentage  is  almost 
always  100%. 

(2)  Stations  in  the  day  sectors  have  50%  strong  responses 
when  AE  is  400-600y.  When  AE  exceeds  600y,  the  propor¬ 
tion  of  strong  responses  occurring  in  pre  noon  sectors 

is  above  80%.  After  noontime,  the  proportion  is 
somewhat  lower,  70%,  for  AE  in  the  600-800y  interval. 

Combining  strong  and  moderate  responses  (lower  panel  in  Figure  2.10)  does  not 
alter  these  results  appreciably.  It  assures  the  detection  of  extremely  large  AE 
enhancements  (>  800y)  at  all  local  times,  but  greatly  reduces  the  ability  to  dis¬ 


criminate  between  large  and  small  enhancements  in  the  night  sectors.  As  a 
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practical  tool  it  appears  that  a  strong  response  is  the  least  ambiguous  to 
determine  and  the  cleanest  to  interpret.  In  the  next  section,  a  more  quantitative 
approach  is  adopted  from  the  start.  Although  the  technique  is  more  refined,  the 
overall  results  are  similar  to  those  already  stated  above. 

D.  Description  and  Examples  of  Mid  Latitude  Variation  Procedure 

One  direct  indicator  of  magnetic  activity  is  the  incremental  change  in 
a  magnetic  component  over  a  specified  time  period.  In  this  section,  we  apply 
a  simple,  quantitative  procedure  to  AFGL  chain  data  that  is  based  on  the 
qualitative  descriptions  of  mid  latitude  activity  presented  earlier.  Our 
procedure  may  be  likened  to  the  production  of  Kp  indices  but  omits  elaborate 
corrections  and  scaling  devices,  and,  in  particular,  uses  shorter  time  scales. 

Our  time  scales  are  1  hour  or  less  compared  to  3  hours  for  Kp  indices,  a  factor 
which  should  eliminate  much  of  the  need  for  quiet  time  corrections.  For  instance, 
a  typical  substorm  variation  at  55°  CGL  is  about  30y  -[n  ]  hour,  while  a  diurnal 
variation  may  amount  to  10-20y  occurring  smoothly  over  3  hours.  During  a  3  hour 
time  interval,  the  substorm  contributes  only  10-20y  above  the  diurnal  variation; 
in  contrast,  over  a  1  hour  interval,  the  substorm  variation  would  again  amount 
to  about  30y,  but  the  quiet  time  change  only  to  about  3-7y. 

Examples  of  AFGL  chain  data  are  given  in  the  remainder  of  this  section. 

Seven  days  of  data  recorded  in  December,  1976  were  scaled,  averaged,  and  plotted 
along  with  a  provisional  AE  index,  as  shown  in  Figures  2.11-2.13.  X  and  Y 
component  data  from  the  AFGL  chain  appear  in  the  upper  two  panels.  Three 
quantities  per  component  are  plotted  in  each  hour:  the  hourly  variation 
averaged  across  the  chain,  plus  the  maximum  and  minimum  hourly  variations.  An 
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Compar  i  son  of  AKCL  chain  data  (upper  two  panels)  with  a  modified  AM 
index  and  Kp  index  on  December  9-10,  197b.  Crosses  connected  hv  solid 
lines  in  the  AX  and  A Y  panels  indicate  the  average  of  the  hourlv 
variations  from  the  chain  (procedure  described  in  the  text).  Dots 
above  and  below  each  cross  indicate  the  maximum  and  minimum  hourlv 
variation  from  the  chain.  Crosses  in  t  lie  lower  panel  give  an  hourlv 
"AM"  value,  described  in  the  text.  values  shown  bv  dots. 
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hourly  variation  is  defined  as  the  difference  between  the  maximum  and  minimum 
values  attained  during  an  hour  at  a  station.  A  cross  depicts  the  chain  averaged 
variation,  and  a  bar  indicates  the  range  of  variations  measured  across  the 
chain  during  1  hour.  The  range  tells  how  uniform  activity  was  in  the  chain 
sector. 

To  gauge  the  effectiveness  of  this  procedure,  we  constructed  an  hourly 
index  similar  to  the  AE  index.  The  longitudinal  coverage  provided  by  9  auroral 
zone  stations  (see  Figure  2.1  for  locations)  is  not  ideal  but  should  nermit  adequate 
detection  of  most  high  latitude  activity.  First,  the  maximum  hourly  deviation 
of  the  H  component  from  an  estimated  baseline  was  determined  at  each  station. 

Then,  out  of  the  9  hourly  deviations  generated  every  hour,  the  largest  was 
selected  as  an  appropriate  value,  an  adaptation  of  the  envelope  of  superposed 
magnetograms  used  in  true  AE  determinations. 

The  general  behavior  of  the  mid  latitude  components  in  Figures  2.11-2.13 
is  quite  similar,  except  for  somewhat  greater  variability  in  the  Y  component 
during  disturbed  periods.  Also,  mid  latitude  variations  during  the  night  (00-12  UT) 
are  clearly  different  than  those  during  the  day  (12-24  UT)  .  The  midpoint  of 
the  chain,  near  WI,  passes  local  midnight  at  approximately  06  UT  and  local  noon 
around  18  UT.  A  quiet  level  of  X  component  variations  appears  to  lie  near  5-10y. 
Large  auroral  zone  disturbances  (>400y)  tend  to  produce  X  variations  of  40y 
or  more  when  the  chain  sector  is  at  night;  in  the  day  sector,  X  variations  are 
smaller  but  still  noticeable. 

Although  a  general  correlation  obviously  exists  between  mid  latitude 
variations  and  high  latitude  magnetic  activity,  there  is  some  disagreement  in 
the  examples  given.  Particularly  outstanding  is  the  rapid  decrease  in  mid 
latitude  variations  after  12  UT  on  December  29  (Figure  2.13).  Also,  the  high 
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latitude  event  at  09  UT  on  December  30  is  reflected  poorly  in  the  AFGL  variations, 
while  a  smaller  event  one  day  later  at  03  UT  is  marked  by  distinct  mid  latitude 
variations.  Such  discrepencies  underscore  the  fact  that  the  quantities  plotted 
here  do  not  provide  unambiguous  descriptions  of  simultaneous  auroral  zone 
and  mid  latitude  magnetic  activity.  Improved  descriptions  might  result  from 
more  complicated  considerations,  such  as  use  of  more  complex  functions  involving 
X  and  Y  components,  different  time  scales,  or  analysis  of  frequency  spectra. 

At  this  stage,  however,  the  natural  next  step  is  to  explore  the  effectiveness 
of  using  shorter  time  intervals. 

Before  investigating  finer  time  scales,  we  present  data  from  FB  and  BD 
which  was  treated  in  the  same  way  as  the  AFGL  data  in  Figures  2.11-2.13.  Hourly 
H  component  variations  from  FB  and  BD  are  shown  in  Figures  2.14  a-d  (upper  panels) 
along  with  2.5  minute  AE  plots  for  the  entire  month  of  December,  1972.  During 
periods  of  prolonged,  extremely  low  AE,  the  hourly  variations  were  generally 
less  than  lOy  and  often  less  than  5y  (Dec.  5-6,  10,  12,  25,  27-28).  However, 
even  with  1  hour  time  intervals,  variations  also  resulted  from  smooth  changes 
associated  with  the  Sq  focus  fixed  near  local  noon.  Clear  cases  are  seen  around 
18  UT  on  December  1,  6,  21,  and  27.  A  characteristic  double  wave  appearance 
is  evident  as  first  FB  and  then  BD  rotated  through  the  H  component  depression. 

At  the  other  extreme,  periods  when  AF,  was  above  500y  are  reflected  clearly  by 
mid  latitude  variations.  There  were  no  cases  of  very  high  AE  when  one  or  both 
mid  latitude  variations  were  under  10y.  Intervals  of  sustained  AE  activitv 
produced  the  greatest  mid  latitude  activity,  as  seen  on  December  13,  15-16, 
and  22-24.  AE  episodes  reaching  700-1000y  generally  were  accompanied  by 
variations  of  20y  or  more,  regardless  of  the  local  time  sector  being  monitored 
at  mid  latitudes.  Marked  variations  in  the  day  sector  occurred  with  activitv 
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during  the  latter  half  of  the  UT  day  on  December  13,  15-16,  22-23,  and  30. 
Finally,  just  as  our  qualitative  studies  suggested,  AE  events  of  lesser  magnitude 
(<500y)  are  not  always  detectable  in  the  mid  latitude  day  sector.  A  comparison 
between  data  on  December  5  and  December  18  reveals  similar  looking  mid  latitude 
variations  occurring  between  12  and  24  UT,  in  spite  of  the  very  dissimilar 
auroral  zone  conditions  on  the  two  days . 

E.  Comparison  of  Short  Time  Scale  AFGL  Data  with  Auroral  Zone  Magnetograms 


AFGL  date  is  now  readily  accessible  in  digital  form  on  a  two-minute  averaged 
time  scale.  These  data  can  be  used  to  calculate  magnetic  variations  over  any 
longer  time  scale  desired.  Although  1  hour  intervals  were  sufficient  for 
preliminary  work  with  plotted  data  and  normal  magnetograms,  such  a  time  scale 
is  too  coarse  to  exploit  the  full  potential  of  good  quality  AFGL  data.  To  move 
toward  a  finer  time  scale,  we  show  in  figures  2.15  a-e  the  X  component  variations 
every  10  minutes  during  June  2-5,  1978.  Except  for  the  finer  time  scale, 
variations  are  calculated  in  exactly  the  same  way  as  those  appearing  in  Figures 
2.11-2.13.  It  is  interesting  to  observe  the  improvement  in  time  resolution  if 
one  compares  Figures  2.15  a-e  with  Figures  2.16  and  2.17  a,  b  which  treat  the 
data  on  June  2-3,  1978,  with  1  hour  and  \  hour  intervals,  respectively. 

Auroral  zone  magnetic  activity  during  June  2-3  is  shown  in  Figure  2.18. 

A  provisional  AE  index  gives  the  maximum  value  achieved  in  each  hour  and  is 
derived  from  4  stations,  TI,  DI,  GWR  and  FC.  As  seen  earlier  in  this  study, 
agreement  between  peak  hourly  AE  and  hourly  mid  latitude  variations  (Figure  2.16) 
is  fairly  good.  Note  the  correspondence  between  peaks  in  Figures  2.16  and  2.18 
near  10,  14  and  19  UT  on  June  2  and  near  00  UT  on  June  3.  A  large  AE  enhancement 
at  19  UT  on  June  3  is  clearly  registered  in  Figure  2.16  with  smaller  variations 
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.  1  8  Maximum  AK  neb  ieved  in  each  hour,  Juno  2  (0300  l!T)  -  June  3  (2100  UT)  ,  1*178. 
Four  auroral  zone  stations  employed,  Tixie,  Dixon,  Oreat  Whale  River, 
and  Fort  Churrhil 1  . 
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than  one  day  before. 

Two  notable  cases  of  improved  time  resolution  are  apparent  in  Figure  2.15  a-e. 
A  100y  variation  occurring  between  09  and  10  UT  on  June  2  (Figure  2.16)  actually 
commenced  sometime  between  0910  and  0920  UT  in  Figure  2.15a.  Later  in  the 
same  day,  a  large  hourly  variation  between  1800  and  1900  UT  did  not  take  place 
until  1840  UT. 

On  a  10  minute  time  scale,  it  becomes  profitable  to  refer  directly  to 
auroral  zone  magnetograms  (Figures  2.19  a-e).  Although  two  stations  alone 
cannot  observe  everything  taking  place,  TI  and  GWR  are  suitably  located  nearly 
12  hours  apart  in  local  time  and  should  give  at  least  an  indication  of  major 
activity.  Scanning  the  AFGL  variation  data  in  Figures  2.15  a-e,  we  see  that 
activity  rose  on  June  2  at  0910,  1000,  1350,  and  1840  UT;  on  June  3,  at  00  and 
0320  UT;  on  June  4,  at  2230  UT;  and  on  June  5,  at  0500  UT.  Each  of  these 
occasions  was  accompanied  by  an  identifiable  high  latitude  event  ranging  in 
activity  from  relatively  mild  (June  3,  00-03  UT)  to  severe  (June  2,  09-21  UT; 

June  4-5,  20-12  UT)  as  seen  in  Figures  2.19  a-e. 

Inevitable  exceptions  in  the  overall  agreement  between  Figures  2.15  and  2.19 
cannot  be  overlooked.  On  June  3,  at  about  0330  UT,  the  AFGL  data  was  quite 
disturbed,  but  there  w,.s  curiously  little  activity  at  GWR  and  at  Fort  Churchill, 
even  though  both  stations  were  in  the  evening  sector.  In  fact,  there  was  more 
activity  at  TI  and  DI  than  in  the  morning  sector.  On  June  4,  beginning  at  03 
UT,  the  AFGL  data  became  very  quiet  while  precisely  at  that  time  the  Canadian 
stations  began  to  detect  negative  bays  of  about  250y  in  the  midnight  sector. 

We  have  presented  examples  demonstrating  the  ability  of  magnetic  variations 
at  mid  and  subauroral  latitudes  to  indicate  the  state  of  auroral  zone  magnetic 
activity.  Even  when  observed  in  a  restricted  longitude  sector  of  the  size  spanned 
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by  the  AFGL  chain,  large  substorms  with  AE>400y  generally  produce  substantial 


variations  on  1  hour,  \  hour  and  10  minute  time  scales.  However,  weaker  substorms 
yield  less  consistent  results,  ranging  from  barely  detectable  to  unexpectedly 
large.  Some  failures  undoubtedly  stem  from  gaps  in  the  AE  network  or  in  the 
particular  array  of  auroral  zone  stations  used  to  identify  substorms.  Unless 
high  latitude  coverage  is  extensive,  an  intense  but  highly  localized  disturbance 
may  affect  some  favorably  located  stations  and  give  a  false  impression  of 
activity  in  the  auroral  zone,  and  consequently  raise  the  expectation  of  mid 
latitude  activity.  Insufficient  station  coverage  does  not  explain  satisfactorily 
the  converse  problem  of  large  mid  latitude  magnetic  variations  but  little  or  no 
auroral  zone  activity.  In  this  case,  the  procedure  shows  a  lack  of  discrimination 
between  substorm  effects  and  variations  coming  from  other  sources.  Fortunately, 
based  on  the  limited  number  of  examples  in  this  study,  occurrences  of  this 
problem  appear  to  be  relatively  rare. 

F.  Summary  and  Estimate  of  Reliability 

In  the  preceeding  sections  of  this  chapter,  two  related  methods  for  deducing 
high  latitude  magnetic  activity  from  mid  latitude  data  were  developed  and  evaluated 
The  first  procedure  employs  a  qualitative  3-level  scale  of  magnetic  response  which 
was  applied  to  magnetograms  from  Fredericksburg  and  Boulder  over  a  4  month  period. 
Results  from  this  preliminary  study  demonstrated  that  such  a  scheme,  although  quite 
unsophisticated,  can  detect  most  substorm  activity  producing  AE  enhancements  above 
AOOy  if  the  mid  latitude  stations  are  in  the  night  sector,  or  above  600y  if  the 
stations  are  in  the  day  sector. 

The  second  method  involves  net  excursions  of  the  magnetic  X  component  recorded 
at  AFGL  stations.  Figure  2.20  illustrates  the  procedure  with  data  from  the 
Wisconsin  station  taken  over  two  8-hour  intervals  in  December,  1976;  boxes  de¬ 
lineating  hourly  excursions  are  transferred  to  the  AX  plot  directly  below.  AX  plot 
from  the  other  AFGL  stations  also  can  be  constructed,  so  that  in  each  time  interval 
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up  to  five  basic  measurements  (AX^,  AX^ ,  AX^,  AXgD,  and  AX^)  can  be  computed, 

depending  on  data  availability.  From  these  quantities  average  AX,  maximum  AX, 

and  minimum  AX  (AX  ,  AX  ,  AX  )  can  be  determined  easily  and  plotted,  as 

in  Figures  2.15  -  2.17  (Section  E).  Treating  in  this  manner  the  limited  amount  of 

AFGL  data  available  to  us  at  the  time  of  this  study,  it  was  found  that  10  minute 

time  intervals  vield  appreciable  values  of  AX  on  the  order  of  lOy  or  more  when 

ave 

AE  exceeds  400v.  This  performance  holds  at  all  local  times. 

Although  not  used  explicitly  in  our  analysis,  the  extreme  values  of  AX 

(AX  and  AX  .  )  indicate  the  uniformity  of  disturbance  across  the  chain  sector, 
max  min  ’ 

and  also  provide  some  protection  against  bad  data  from  a  single  station.  For 

example,  a  situation  where  AX^  =  lOy,  AX ^  -  20y,  AX^  =  30y,  AX^  =  40,,  and 

AX,  =  50>  vields  AX  =  )0> ,  AX  ,  =  10y,  and  AX  =  50v .  If  a  bad  value  from 

WA  ave  min  max 

WA  makes  AX  =  500y,  instead  of  50y,  then  AX  =  120>,  AX  ,  =  10y,  and  AX  = 

WA  ave  min  max 

5Q0y.  The  high  value  of  "'^ave »  120y,  may  not  be  unusual,  but  the  large  difference 

between  AX  -  AX  and  AX  -  AX  ,  11 0y  vs.  480y,  would  warn  that  a  station 

ave  min  max  ave 

value  might  be  in  error. 

Although  there  are  not  enough  examples  to  determine  the  statistical  reliability 
of  the  quantitative  procedure  outlined  above,  data  from  the  AFGI.  chain  should  equal, 
if  not  exceed,  the  performance  indicated  by  data  from  Fredericksburg  and  Boulder 
(Figures  2.10  and  2.14  a-d).  At  55°  CGL,  the  AFGL  chain  is  closer  to  substorm 
current  systems  than  FB  and  BD  and  spans  an  additional  1.5  hours  in  local  time. 

Both  factors  make  it  reasonable  to  expect  AFGL  chain  data  to  be  even  more  sensitive 
and  reliable  than  the  FB  and  BD  results.  More  extensive  investigations  with  chain 
data  are  needed  to  assess  seasonal  effects  and  wider  ranges  of  magnetic  activity. 

As  experience  is  gained,  the  procedure  can  be  modified,  such  as  by  including  tin 
Y  component,  varyini  the  time  interval,  or  utillzying  the  sign  of  the  excursion. 

Table  2.2  summarizes  the  procedure  and  lists  some  of  its  advantages  and  dis¬ 
advantages.  Sensitivity  and  reliability  estimates  are  based  on  qualitative  and 
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and  quantitative  results  gathered  from  FB  and  BD.  Examples  from  the  AFGL  chain 
in  December,  1976  and  June,  1978  are  consistent  with  these  estimates. 
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Pro<  edure:  1.  Choose  time  interval  (e.g.,1  hour,  h  hour,  10  minutes) 

2.  h.-t.  rm i in-  X  component  excursion  during  l  Line  interval  (e. 

“VlA  "'max-'  'min)’  whcre  W  andXmin  are  the  maximu,n 
aim  minimum  values  attained  (.may  occur  at  beginning  or  end 

points  of  interval,  e.f..  Figure  2.20). 

1.  Calculate  AXave  (average  of  AX  from  stations  reporting  during 
the  interval. 

4.  Plot  'X.,ve  against  time.  Range  of  A X  given  by  AXmax  and 
4  min  from  stations  reporting  during  each  interval  (c.f.. 
Figures  2.15  a-e,  10  minute  intervals). 

Interpretat ion:  AFGL  chain  (55°  CGL) ,  Sections  D  and  E 

aX'vyillO)  (10  minute  interval)  generally  implies 
AE-'v400y  (c.f..  Figures  2.15  n-e.  Figure  2.16). 

FB  and  BD  (50°  CGL),  Sections  B  and  C 
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fable  2.2  Summary  of  procedures  and  interpretations  discussed  in  Chapter  2 
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3.  MICROPULSATIONS  AS  INDICATORS  OF  GEOMAGNETIC  ACTIVITY 

A.  Introduction 

The  question  of  correlation  between  geomagnetic  activity  and  micro¬ 
pulsations  (ultra-low-frequency  oscillations  in  the  geomagnetic  field)  is 
an  intriguing  problem  in  magnetospiieric  physics.  I t  is  generally  accepted 
that  geomagnetic  activity  is  driven  bv  the  solar  wind  (cf.,  Burch ,  1974; 

Crooker,  1975;  and  reviews  therein,  for  background).  In  addition,  micro¬ 
pulsation  occurrence  correlates  reasonably  well  with  solar  wind  parameters 
(cf.,  Raspopov  and  Lanzerotti,  1976;  Lanzerotti  and  Fukunishi,  1974  for 
theoretical  background;  and  Troitskaya,  1967;  Sai to ,  1969;  Jacobs ,  1970; 
and  Or r ,  1973;  for  observational  inf orm.it  ion) .  Thus,  since  both  geomagnetic 
activity  and  micropulsations  are  strongly  related  to  the  solar  wind  and  correlate 
well  with  many  of  its  parameters,  it  follows  that  the  two  phenomena  should 
correlate  with  each  other.  As  will  be  shown  later,  reasonably  good  correlation 
between  micropulsation  activity  and  geomagnetic  activity  has  been  found  in  a 
number  of  studies. 

In  this  summary,  we  will  first  review  the  background  literature  on  (i) 
correlations  between  solar  wind  parameters  and  continuous  davside  micropul¬ 
sations;  (ii)  relationships  between  the  solar  wind  and  geomagnetic  activity; 
and  (iii)  the  relation  between  micropulsations  and  geomagnetic  activity.  Then 
we  will  discuss  three  recent  sets  of  observations  which  lend  strong  support 
to  the  relationships  suggested  earlier  in  the  literature.  Finally,  we  will 
suggest  ways  in  which  the  data  from  the  AFGL  magnetometer  network  may  be  used 
to  formulate  a  micropulsation  index  which  would  reflect  the  state  of 
magnetospher i c  activity. 


Correlative  studies  have  shown  strong  relationships  between  continuous 
micropulsations  and  solar  wind  parameters.  These  include  solar  wind  velocity 


(Troitskaya ,  1967),  magnitude  of  the  IMF  (Gul'yel’mi  and  Bol'shakova,  1973), 
and  the  angle  between  the  IMF  and  the  sun-earth  line  (Plyasova-Bakunina,  1972) 
The  work  by  Plyasova-Bakunina  is  in  agreement  with  external  wave 
generation  at  the  bow  shock,  as  predicted  by  Greenstadt's  mechanism;  the 
correlation  is  shown  in  Figure  3. 1.  Note  that  Plyasova-Bakunina  shows  a  high 
occurrence  frequency  of  Pc  2-4  when  the  projection  of  the  angle  between  the 
IMF  and  the  sun-earth  line  is  small.  This  correlation  was  supported  by  the 
observations  of  Greenstadt  and  Olson  (1976)  using  the  amplitudes  of  Pc  3-4 
measured  at  Calgary  (L  =  3.6)  and  the  angle  between  the  IMF  and  the  sun-earth 
line. 


m*. 


Figure  3.1.  From  Plyasova-Bakunina  (1972):  Occurrence  frequency  of 

pulsations  for  intervals  of  angle  y  measured  in  the  solar 
wind  at  Imp  3,4.  y  is  the  projection  onto  the  ecliptic 
plane  of  the  angle  between  the  IMF  and  sun-earth  line. 


B .  General  Background 

1)  The  Solar  Wind  and  Micropulsations.  It  was  originally  suggested 
by  Dungey  (1954)  that  micropulsations  are  the  result  of  hydromagnetic  oscil¬ 
lations  of  the  earth's  geomagnetic  field  in  the  magnetosphere.  Following 
that  work,  micropulsation  theory  developed  around  the  ideas  of  field  line 
resonance  and  resonance  of  the  magnetospher ic  cavity.  Observations  from 
latitudinal  chains  of  magnetic  observatories  stimulated  the  extension  of 
those  resonance  theories  to  coupled  waves  to  take  account  of  non-uniformities 
in  the  plasma  and  field  distributions,  and  to  account  for  latitude  localization 
of  micropulsation  activity  (Radoski ,  1974;  Southwood ,  1974;  Lanzerotti  et  al . , 
1974;  Chen  and  Hasegawa,  1974a, b).  Subsequent:  observations  supported  the 
theory  that  Pc  1-5  micropulsations  are  the  result  of  excitation  of  resonant 
field  lines  inside  the  magnetosphere  by  oscillations  at  the  magnetopause 
(Ar tjhu r__e t_ al . ,  1977;  H ughcs  et  al . ,  1977). 

If  micropulsations  are  generated  by  oscillations  on  the  magnetopause, 
what  is  the  original  driving  force?  Two  mechanisms  have  been  proposed  for 
this  driving  force.  The  first  is  the  Kel vin-Helmhol tz  instability  generated 
at  the  magnetopause  by  the  flow  of  magnetoshcath  plasma  along  the  flanks  of 
the  magnetosphere  (Southwood,  1968).  Under  the  Kolvin-Helmhol tz  mechanism, 
waves  would  lie  generated  given  a  high  solar  wind  speed  and  Large  angle  between 
the  mngnotosheath  flow  direction  and  the  magnetic  field  directions  on  either 
side  of  the  magnetopause.  The  second  magnetopause  wave  generation  mechanism 
is  the  direct  transmission  of  waves  through  the  mugnetosheath  from  the  bow 
shock  (Grecnstadt,  1972).  Under  tiiis  mechanism,  waves  would  be  generated 
preferentially  for  small  angles  between  the  sun-earth  line  and  thp  direction 
ol  the  solar  wind  magnetic  field. 
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ii)  The  Solar  Wind  and  Geomagnetic  Activity.  It  is  generally  accepted 
that  geomagnetic  activity  is  driven  by  the  solar  wind.  Whereas  the  mechanisms 
by  which  energy  is  transferred  are  not  well  understood,  a  number  of  correlation 
studies  (cf.,  Burch,  1974;  Crooker ,  1975)  strongly  indicate  that  the  transfer 
mechanism  involves  merging  between  the  IMF  and  the  geomagnetic  field.  The 
southward  component  of  the  IMF  and  solar  wind  speed  give  excellent  correlations 
with  geomagnetic  activity  (Crooker  et  al.,  1977;  Burlaga  and  Lepping,  1977; 
and  others) . 

iii)  Micropulsations  and  Geomagnetic  Activity.  Both  dayside  continuous 
nicropulsations  and  geomagnetic  activity  correlate  well  with  many  solar  wind 
parameters.  However,  the  two  correlate  differently  and  with  different  sets  of 
parameters,  so  that  a  direct  theoretical  relationship  between  geomagnetic  activity 
and  micropulsations  is  difficult  to  formulate. 

Micropulsations  made  at  Borok  (L  =  2.9)  have  shown  a  good  correlation 
between  micropulsation  period  and  Kp  (Troitskaya,  1967).  Figure  3.2  shows 
those  results.  60 
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Figure  3.2  From  Troi tskaya  (1967):  Dependence  of  Pc  2-4  periods  on  Kp  at  Borok. 

C-,  Recent  Results 

Three  recent  pieces  of  work  have  given  significant  credibility  to  the 
mechanisms  previously  proposed  for  pulsation  generation,  and  have  given 
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support  to  the  possibility  of  producing  a  comprehensive  micropulsation 
Index  which  would  indicate  the  level  of  geomagnetic  activity.  Two  of  the 
studies  were  based  upon  satellite  observations  and  the  third  was  based  upon 
ground  micropulsation  observations.  Most  importantly,  all  three  studies  used 
digital  data,  which  allowed  flexible  data  manipulation  and  sophisticated 
computer  analysis. 

Recent  simultaneous  observations  of  geomagnetic  pulsations  using  three 
synchronous  satellites  (Hughes  et  al.,  1978)  were  in  strong  agreement  with 
generation  of  the  waves  by  the  Kelvin-Helmholtz  mechanism.  The  data  were 
collected  in  February,  1975  when,  for  a  few  days,  the  three  geosynchronous 
satellites  ATS  6,  SMS  1,  and  SMS  2  were  operating  in  close  proximity.  The 
data  were  cross-spectral-analyzed  to  obtain  phase  differences  between  the 
spacecraft.  All  periods  of  pulsation  activity  observed  on  the  dayside  of  the 
magnetosphere  during  one  week  were  analyzed.  Figure  3.3  shows  the  results 
of  power  spectral  analysis  upon  a  typical  micropulsation  event  observed  during 
the  local  morning.  It  should  be  noted  that  the  same  spectral  peak  was  observed 
at  all  three  spacecraft  (the  separation  between  adjacent  satellites  was  about 
10°)  and  that  the  signals  were  highly  coherent. 

During  the  local  morning  and  early  afternoon  hours,  it  was  found  that 
when  there  was  a  significant  peak  in  the  power  spectrum  at  two  satellites, 
there  was  a  fair  degree  of  coherence  between  the  two  signals  to  enable  calculation 
of  an  m  number  (representing  azimuthal  wavelength).  In  local  afternoon  after 
1500,  however,  this  coherence  often  did  not  apply. 
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Figure  3.3.  Spectra  and  cross-spectral  coherence  of  the  radial  and  azimuthal 

magnetic  fields  obtained  by  the  three  space  craft  during  a  typical 
local  morning  pulsation.  From  Hughes  et  al.,  (1978). 

Figure  3.4  shows  the  variation  of  m  number  obtained  as  a  function  of  time. 

The  results  provide  support  for  the  Kelvin-Helmholtz  wave  generation  mechanism. 

In  the  vast  majority  of  the  morning  events,  the  more  easterly  satellite 

leads  in  phase,  while  around  noon  this  trend  is  reversed,  and  in  the  early 

afternoon  the  more  westerly  satellite  leads.  This  indicates  wave  propagation 

away  from  noon,  in  the  same  direction  as  the  flow  of  the  solar  wind  plasma 

around  the  magnetopause. 
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Figure  3.4.  The  variation  of  the  azimuthal  wave  number  (m)  as  a  function  of 
time,  ‘'can  local  noon  is  computed  at  the  midpoint  between  ATS6 
and  SMS  1 .  From  Hushes  et  al. ,  (1978) . 

The  lack  of  coherence  between  the  signals  observed  at  two  different 
satellites  in  the  late  afternoon  indicates  a  strong  localization  in  longitude. 
A  possible  source  of  those  pulsations  could  be  hot  protons  injected  from 
the  nightside.  Thus,  such  pulsations  would  be  expected  to  be  associated 
with  individual  substorm  events,  and  the  power  spectral  density  of  those 
waves  would  be  an  indicator  of  geomagnetic  activitv  level. 

A  recent  study  by  Creenstadt  et  al .  (1979)  indicates  that  both  the 
Kelvin-Helmholtz  mechanism  and  direct  injection  of  waves  through  the  magneto¬ 
pause  are  important  sources  of  micropulsation  energy.  In  that  study,  hourly 
distribut ions  of  micropulsation  activity  orientation  angle  were  correlated 
both  with  the  solar  wind  alignment  of  the  IMF  and  with  solar  wind  speed. 

As  shown  in  Figure  3.3,  the  results  showed  positive  point  correlation  with 
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both  speed  and  orientation.  In  addition,  the  joint  correlations  showed 
that  enhanced  magnetic  pulsation  amplitudes  were  associated  with  solar  wind 
speeds  above  200-300  km/sec  and  angles  between  the  IMF  and  sun-earth  line  of 
less  than  50°-60^.  The  threshold  effects  further  support  both  the  bow 
shock  origin  and  the  Kelvin-Helmholtz  generation  of  daytime  micropulsation 
signals . 


'‘BM  Ft  *i 


Figure  3.5.  Diagrammatical  summary  of  joint  SB  versus  cosO^g  and  Vgw 
correlations.  From  Greenstadt  et  al.  (1979). 

In  another  pulsation  study  using  data  from  the  geosynchronous  satellite 
ATS-6 ,  Arthur  et  al.  (1978)  calculated  the  power  spectral  density  of  geomagnetic 
field  variations  for  the  month  of  August,  1974.  Two  segments  of  approximately 
9  hours  each,  centered  at  local  noon  and  at  local  midnight,  were  used  for 
the  analysis.  Figure  3.6  shows  the  time  variation  of  the  power  spectral 
density  as  well  as  the  7Kp  for  the  month  of  data.  These  results  have 
striking  implications  on  the  correlation  between  micropulsation  power  content 
and  geomagnetic  activity.  It  is  significant  that  the  power  spectral  density 
is  directly  related  to  EFp  on  both  the  dayside  and  the  nightside.  Smaller 
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(larger)  i'Kp  corresponds  to  smaller  (larger)  values  of  power  spectral  density. 

Arthur  et  al .  also  compared  the  average  characteristics  of  power  spectra 
at  synchronous  orbit  with  observations  from  other  regions  of  Interest  in 
magnetospher ic  physics:  the  interplanetary  medium,  the  magnetosheath,  the 
outer  magnetosphere  ne.ir  the  magnetopause,  and  the  surface  of  the  earth. 

Figure  3.7  shows  the  average  spectrum  for  each  of  the  five  regions.  The 


Figure  3.7.  Average  power  spectra  for  the  five  regions,  interplanetary, 
magnetosheath,  magnetopause,  synchronous  orbit,  and  ground. 
From  Arthur  et  al.  (1978). 

results  in  Figure  3.7  indicate  that,  except  for  the  magnetosheath ,  which  is 
known  to  he  verv  disturbed,  the  power  levels  in  each  region  increase  from 
the  interplanetary  medium  inward  to  the  surface  of  the  earth.  Thus,  given 
a  good  source  of  comprehensive  digital  data,  ground  observations  should 
be  at  least  as  sensitive  an  indicator  of  geomagnetic  activity  as  these 
reported  by  Arthur  et  a  1 . 

D.  Summary  and  Reeommendat ions 
The  results  indicated  herein  suggested  that: 

i)  both  the  Kelvin-Helmhol tz  mechanism  and  direct  infection  of 

waves  from  the  how  shock  are  important  sources  of  micropulsat ion 
energy  (Hughes _ et  al.,  19  78:  t'.reenstadt  and  Olson,  1979) 


i  i )  m  i  eropulsa  t  ion  period  is  a  possible  Indicator  ol  geomagnet  i  c 
activity  level  (Troitskava,  19b7), 
iiL)  and  the  power  spectral  density  oi  micropu  1  sat  ions  :nav  he  used 

as  an  indicator  of  geomagnetic  activitv  level  (Art  luir  el_  al  .  ,  1 11 7  S)  . 

The  AFCl.  magnetometer  network  is  an  ideal  tool  for  formulating  a  micro- 
pulsation  index  for  the  level  of  geomagnetic  activitv.  The  AFt'd  data  is  also 
ideal  for  comparison  with  observations  from  other  chains  of  magnetometers. 

Since  the  Al'lll,  data  is  digital,  it  would  hi’  straightforward  to  perform  power 
spectral  analysis.  And  since  all  facilities  used  for  the  network  are  dedicated 
to  the  network,  the  data  are  available  uninterrupted  for  long  periods  of  time. 

Two  possible  approaches  may  be  used  to  formulaic  .and  calibrate  a  micro¬ 
pulsation  index  from  the  AFt'd.  magnetometer  network  data. 

The  first,  and  simplest,  approach  would  he  to  pass  the  data  through 
a  suite  of  tillers,  each  passband  carefully  chosen  to  cover  a  particular 
part  o!  the  mi crepu 1 sa t  i on  spectrum.  The  output  of  each  filter  would  then 
he  time  averaged  .and  correlated  with  geomagnetic  activity,  o.g.,  Kp .  After 
appropriate  studv.  the  filter  hand  passes  could  he  precisely  defined  so  as 
to  cover  the  spectrum  in  a  way  that  the  most  significant  channels  could  he 
isolated.  At  that  point,  the  super; luotis  channels  could  be  eliminated  (although 
they  probably  should  be  retained  for  future  studvl. 

A  second,  mere  soph i st i cared ,  approach  would  hi'  as  follows.  Compute  the 
power  spectril  densit\  matrix  tor  at.  appropriate  time  interval  (one  hour, 
as  i  start  1.  Ketuin  the  power  :  r.)  1  density  for  each  component  oi  the 

m  ic.net  ie  i  iold.  Caretul  comparison  of  shifts  ot  spectral  peaks  with  changing 
y  eotr.agnot  ic  activity  would  verify  or  refute  the  possibility  ot  using,  micropul  sat  ion 
period  as  the  kev  indicator  ot  activity  level. 
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Following  tile  technique  of  Arthur  et  al .  (1978),  the  level  of  tiie  spectrum 
(characterized  by  key  spectral  estimates')  could  be  correlated  with  geomagnetic 
activity.  Based  on  the  success  of  this  approach  t  synchronous  orbit,  it  is 
likely  that  a  micropulsation  index  using  this  technique  on  AFGL  magnetometer 
network  data  could  be  a  reliable  indicator  of  geomagnetic,  activity. 

In  this  report,  we  have  emphasized  dayside  continuous  micropulsations  as 
possible  indicators  of  geomagnetic  activity.  However,  it  is  apparent  that 
irregular  micropulsations,  especially  Pi2,  may  also  be  good  indicators  of  the 
level  of  activity.  Pi2  are  transient  pulsations  associated  with  substorms,  with 
maximum  effect  at  middle  and  low  latitudes.  Pi2  are  known  to  correlate  with  Kp 
(Troitskaya,  1967).  A  possible  future  study  of  the  potential  of  Pi  2  occurrence 
as  a  geomagnetic  activity  index  using  the  AFGL  network  data  should  be  considered 
at  some  point. 
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4.  IONOSPHERIC  ELECTRIC  FIELDS  AND  CURRENTS  DRIVEN  BY  BIRR ELAND  CURRENTS 

A.  Introduct ion 

A  significant  observational  finding  of  the  past  decade  is  the  discovery 
of  magnetospheric  electric  currents  that  flow  along  magnetic  field  lines 
and  enter  the  ionosphere  in  the  auroral  regions.  These  so-called  Birkeland, 
or  f ield-aligned  currents,  constitute  a  major  coupling  between  the  magnetosphere 
and  the  ionosphere.  Other  recent  observations  suggest  that  certain  ionospheric 
phenomena,  that  may  occur  at  all  latitudes  from  the  polar  cap  to  the  equator, 
are  causal Ly connected  to  the  electric  field,  current,  and  conductivity 
variations  that  occur  in  the  auroral  regions.  For  the  purpose  of  eventually 
relating  observed  mid  latitude  field  signatures  to  magnetospheric- ionospheric 
current  configurations,  we  have  developed  a  model  for  the  global  distribution 
of  ionospheric  currents  generated  by  a  specified  input-output  Birkeland 
current  configuration.  The  model  developed  is  discussed  in  this  chapter  and 
some  principal  results  are  included  in  the  three  publications  reprinted  at 
the  end  of  this  chapter. 

The  model  developed  of  global  distributions  of  electric  fields  and 
currents  driven  by  field  aligned  currents  is  useful  in  several  ways.  To  cite 
three  reasons:  (1)  Within  the  model,  the  ionospheric  distribution  of  electric 
currents  that  result  from  various  models  of  geomagnetic  events,  sucli  as 
suhstorm.s,  are  determined,  and  from  these  currents,  ultimately,  the  magnetic 
disturbance  patterns  can  then  he  calculated  for  such  events.  These,  modeled 
events  are  then  used  in  conjunction  with  actual  records  (magnetograms,  say , 
at  middle  latitude)  to  infer  the  behavior  of  the  high  latitude  ionosphere. 

{•a rrying  out  this  procedure  is  discussed  in  more  detail  later  in  this  chapter.) 
(.')  Tin-  electric  fields  in  the  ionosphere  are  determined  globally  for  model 
events.  Since  a  great  many  plasma  phenomena,  c.g.,  at  the  dip  equator,  depend 
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upon  the  local  electric  field,  one  could  use  such  a  model  to  determine  the 
relationship  between,  say,  equatorial  ionospheric  plasma  observations  and 
Birkeland  currents.  (3)  This  gives  quantitative  insight  into  the  relative  con¬ 
tributions  of  various  electrodynamic  effects  that  are  operative  (such  as 
Birkeland  currents,  neutral  winds,  conductivity  structure,  etc.)  in  determining 
the  resultant  distributions  of  ionospheric  electric  fields  and  currents. 

Thus,  in  summary,  the  kind  of  modeling  developed  serves  to  enhance  our  under¬ 
standing  of  the  connection  between  magnetospheric  dynamics  and  various  electro¬ 
dynamic  effects  in  the  global  ionosphere.  This  improved  understanding,  in  turn, 
will  enable  us  to  develop  ways  to  exploit  the  connection  for  practical  purposes. 

B.  Definition  of  the  Model 

The  definition  of  our  model  is  as  follows.  The  ionospheric  height-integrated 
tensor  conductivity  and  the  distribution  of  field-aligned  currents  into  and 
out  of  the  ionosphere  are  specified  (i.e.,  taken  as  known  variables)  for  each 
model  event  considered.  The  basic  condition  imposed  at  each  point  in  the 
ionosphere  is  that  charge  is  conserved.  This  is  done  through  the  continuity 
equation,  which,  using  Ohm's  law,  has  the  form 

I(field-aligned)  =  V 2  •  L(0, y)  •  V2<I>(0,y) 

The  left-hand  side  of  this  equation  is  provided  by  the  given  Birkeland  current 
distributional  the  ionosphere.  <J>(8,y)  is  the  electric  potential  in  the  (height- 
integrated)  ionosphere  at  the  colatitudinal  and  azimuthal  angles  0  and  y.  The 
ionospheric  electric  field  and  current  distribution  are  determined  by  the 
relations : 

E(0,y)  =  -V 2  <5>(0,y)  , 

J(6,y)  =  1(6, y)  •  E(6 ,y) • 

The  subscript  on  the  gradient  indicates  that  only  transverse  derivatives  (tangential 
to  the  spherical  ionosphere)  are  operative. 

Tiie  conductivity  tensor  |L(0,y)  assumed  known  and  includes  the  principal 
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spatial  variations:  a  largo  increase  of  dayslde  over  nights Ide  conductivity 
due  to  solar  photoionization  (i.e.,  the  solar-zenith  angle  variation);  and 
a  large  increase  in  the  auroral  regions  due  to  energetic  particle  precipitation. 
These  variations  are  depicted  schematically  in  Figure  4.1.  In  addition,  the 
layer  conductivity  fonnuLation  employed  leads  to  enhancements  of  the  diagonal 
conductivities  and  a  zero  of  the  off-diagonal  conductivities  at  the  equator 
(as  shown  in  Figure  1  of  the  third  reprinted  paper  in  this  chapter  (Nopner  and 
Carovillano,  1979b).  The  assumed  distribution  of  driving  currents  is  adapted 
from  the  Triad  satellite  observations  (c.f..  Figure  2  of  the  second  reprinted 
paper  in  this  chapter  (Nopper  and  Carovillano,  1979a).  With  these  modeling 
parameters,  we  are  able  to  represent  a  range  of  magnetospheric  conditions 
from  very  quiet  to  quite  disturbed. 

The  numerical  scheme  employed  to  solve  the  continuity  equation  is  an 
iterative  method  using  successive  relaxations.  Iterations  are  made  to  a  finite 
difference  representation  of  the  continuity  equation  on  a  hemispherical  grid. 
This  well-known  technique  (Forsyth  and  Wasow ,  1960;  Binns  and  Lawrenson, 

1971)  vields  reliable  numerical  results,  as  tested  independently  using  Gauss's 
Law  in  two  dimensions  on  the  final  output,  and  gives  accurate  electric  fields 
and  currents  to  within  a  few  percent  almost  everywhere. 

C.  Applications  of  the  Model 

In  this  section,  we  briefly  describe  the  situations  modeled  and  the 
determined  cqui potential  patterns  and  electric  field  distributions.  We  also 
summarize  here  the  results  contained  In  the  three  reprinted  papers  included 
in  section  D  of  tills  chapter.  We  shall  refer  to  those  three  papers  as  Paper  I 
(Nopper  and  Carovillano,  1978),  Paper  II  (Nopner  and  Carovillano,  1979a), 
and  Paper  III  (Nopper  and  Carovillano,  1979b). 
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PHYSICAL  IONOSPHERIC 
CONDUCTIVITIES 


.  I  Schemat  i  c  representation  of  the  distribution  of  ionospheric  conductivitv 
used  in  the  model.  Conductivity  is  proportional  to  density  of  stippling. 
Features  shown  include  the  solar-zenith  (i.e.,  dav-nifiht)  variation, 
and  enhanced  conductivitv  in  the  auroral  ovals. 


-  ,#S. 
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The  basLc  quiet-time  (or  weakly  disturbed)  pattern  of  field-aligned  currents 
discerned  by  the  Triad  satellite  includes  a  "Region  1 "  (pol eward)  pair,  with 
current  inti'  the  dawnside  and  out  of  the  deskside  ionosphere,  at  the  high 
latitude  edge  of  the  auroral  oval;  and  a  "Region  2"  (equatorward)  pair,  with 
opposite  polarity  to  the  Region  1  set,  at  the  equatorward  edge  of  the  oval; 

Region  1  currents,  at  lO^A,  are  typically  twice  the  strength  of  the  Region  2 
currents  during  quiet  periods.  (See  Figure  2a  of  Paper  I  (Nopper  and  Carovi llano, 
1978)  for  our  idealization  of  this  structure;  this  figure  also  shows  the  equi- 
potentials  for  this  set  of  source  currents).  Systematic  findings  arc?  that 
ttie  duskiike  equipotent ials  are  pushed  up  toward  the  pole,  while  the  dawnside 
equ Lpotent ial s  are  drawn  down  toward  the  equator;  in  addition,  the  polar  cap 
equipotent ials  are  pushed  over  toward  the  dawnside  of  the  polar  cap.  These 
distortions  are  understood  (Wol_f,  1970;  Nopper  and  Carovillano,  1979b)  to 
result  from  the  contrast  in  conduct ivitv  between  the  day  and  night  sides  of  the 
earth,  and  are  present  in  all  models  with  such  a  conductivity  structure.  A 
consequence  of  this  asymmetry  appears  in  the  electric  field  distribution.  Thus, 
Paper  I  (Nopper  and  Carovillano,  1978)  shows  the  quiet-time  east-west  electric 
field  at  the  equator  as  a  function  of  local  time;  one  sees  that  the  field 
is  eastward  over  three-quarters  of  the  day  and  westward  for  only  about  six 
hours  from  2100-0100  I.T,  and  that  the  field  is  stronger  during  its  westward 
orientation  than  its  eastward  (a  consequence  of  the  requirement.  that^E-dJ  =  0). 
Quiet-time  electric  field  components  at  latitudes  54°  and  30°,  as  functions 
of  local  time,  are  displayed  in  Figures  4  and  5  of  Paper  II  (Nopper  and  Carovi 1 1 ano , 
1979a)  where  the  local-time  asymmetries  are  also  evident. 

During  disturbed  periods,  observations  indicate  that  the  topology  of  the 
field-aligned  currents  remains  the  same  as  for  quiet  times,  hut  that  the 
strengths  of  the  currents  increase  to  roughly  3x1 0^A  in  both  the  Region  1 
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and  Region  2  pairs.  Thus,  not  only  do  the  currents  increase  in  strength  during 
active  conditions,  but  also  the  Region  2  current  becomes  approximately  as 
strong  as  in  Region  1  (in  quiet  time.  Region  1  dominates).  The  general 
effect  on  the  equipotentia]  pattern  can  be  seen  in  various  figures  in  Papers 
I,  II,  and  III.  Thus,  changes  in  the  electric  fields  are  evidenL  in  Figures 
2a  and  3  of  Paper  I,  and  Figures  8  and  9  of  Paper  II.  The  change  in  equipotentiaL 
pattern  is  seen  in  Figures  2a  and  2c  of  Paper  III.  A  striking  result  is  that, 
during  disturbed  periods,  the  east-west  component  of  the  electric  field  reverses 
its  quiet-time  direction  at  the  equator  (and  other  latitudes)  at  most  local  times. 

Using  models  of  the  field-aligned  currents  during  quiet  and  disturbed 
periods  as  just  discussed,  we  show  in  Paper  I  that  the  source  currents  as 
observed  during  active  times  generate  equatorial  electric  fields  of  sufficient 
strength  and  in  the  proper  direction  to  account  for  the  fluctuations  and 
reversals  of  the  electric  field  seen  in  equatorial  records.  Originally  inferred 
from  magnetic  observations,  these  disturbances  in  the  equatorial  electric 
field  have  now  been  seen  directly  as  actual  ionospheric  effects  by  radar  and 
other  sensing  techniques,  and  often  occur  simultaneously  with  high- latitude 
disturbances.  Thus,  Paper  I  establishes  the  ionospheric  conductivity  as  a 
contributing  mechanism  to  the  coupling  of  electric  fields  from  high  latitudes, 
where  driving  forces  of  magnetospheric  origin  are  applied,  to  low  ionospheric 
latitudes. 

The  consequences  of  global  spreading  of  electric  fields  of  magnetospheric 
origin  are  examined  further  in  Paper  II.  Here,  in  addition  to  the  polar-equatorial 
coupling  discussed  above,  we  conclude  that  the  electric  fields  due  to  Birkeland 
currents  are  comparable  or  exceed  in  magnitude  those  generated  in  neutral 
wind  models  at  all  latitudes,  and  generally  have  large  subdiumal  (12,  8,  ...  hour) 
components.  Tlius,  even  during  quiet  periods,  neutral  wind  models  are  insufficient 
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to  motif  I  ionospheric  electric-  fields  and  currents,  and  the  effects  of 
magnotospherie  sources  must  be  included.  II  was  also  observed  in  Paper  II 
that  during  disturbed  times,  the  electric  field  in  the  polar  cap  deviates  from 
its  canonical  dawn-to-dusk  direction  and  re-orients  toward  the  sun. 

This  last  effect  is  the  topic  of  Paper  III  where  the  polar  cap  electric 
field  is  examined  under  various  disturbance  conditions.  The  principal  result 
is  that  the  structure  of  equ i potent  ial s ,  global Iv  and  in  the  polar  cap,  depends 
sensitively  on  the  relative  strengths  of  the  field-aligned  current  pairs.  For 
Region  1  currents  much  stronger  than  Region  2  (a  very  "quiet"  condition),  tin- 
average  electric  field  direction  in  the  polar  cap  is  about  (i0°  east  of  noon. 

As  the  Region  2  currents  increase  in  strength,  corresponding  to  magnetic 
conditions  becoming  increasingly  active,  the  orientation  of  the  polar  cap 
electric,  field  rotates  clockwise  (as  seen  from  above  the  pole)  towards  noon 
and  would  actually  point  into  the  pre-noon  sector  should  the  intensity  of  Region 
currents  slightly  exceed  Region  )  currents.  In  Paper  111,  Figures  2n-o 
show  tin-  smooth  progression  of  equipotent i  patterns  for  Region  1/Region  2 
intensity  ratios  of  2:1  (quiet  conditions)  through  1:1  (disturbed  conditions) 
to  1:2  ( intensely  and  probably  unreal istical ly  disturbed  conditions).  The 
electric  field  direction,  measured  counterclockwise  f rom  noon,  is  given  in 
Figure  i.2  as  a  function  of  the  intensity  ratio.  Because  the  intensity  of 
Region  2  currents  is  directly  dependent  upon  geomagnetic  activity  conditions, 
an  important  implication  of  our  work  is  that  tin’  electric  field  orientation  in 
the  polar  cap  should  also  correlate  with  geomagnot  ic  activity. 

I).  Reprinted  Publications  (Papers  I,  1),  and  111) 

I  he  reprints  mentioned  in  the  foregoing  text  and  briefly  discussed  in 
Si  cl  ion  C  as  Papers  1,  II,  and  11T  are  presented  in  this  section. 


CURRENT  RATIO  (region  i /region  2) 


Hit-  orientation  of  the  mode]  electric  field  at  the  north  pole, 
measured  eoimterclockwl.se  from  noon,  as  a  function  of  the  ratio 
•’I  total  Uei'ion  1  current  to  total  Region  2  current. 
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POLAR-EQUATORIAL  COUPLING  DURING  MAGNK i ICALLY  ACTIVE  PERIODS 
R.  W.  Nopper  ,  Jr.  and  R.  L.  Carovillanc 
Department  of  Physics,  Boston  College,  Chestnut  Hill,  Massachusetts  02167 


Abstract.  Observational  evidence  Indicates 
that  high^lat i tude  magnetic  activity  can  be 
causally  related  to  fluctuations  and  reversals 
of  the  equatorial  electric  field.  Using  Birke- 
1  and  currents  as  the  driving  forces,  we  have 
made  numerical  calculations  «>t  the  ionospheric 
electric  fields  and  currents  that  distribute 
globally.  Our  calculations  utilize  a  fairly 
realistic  empirical  ionospheric  conductivity 
model  and  the  field-aligned  current  distribution 
Inferred  from  TRIAD  measurements.  The  results 
presented  here  demonstrate  that  the  field- 
aligned  currents  observed  during  disturbed 
periods  o-m  account  for  the  equatorial  fluctua¬ 
tions  and  reversals  in  electric  field  and 
associated  magnetic  disturbances  noted  during 
such  times.  In  this  way,  magretospher ic  dyna¬ 
mics  has  a  direct  and  significant  Influence  on 
the  equatorial  ionosphere  through  the  iono¬ 
spheric  conductivity. 

Introdut  t ion 

A  body  of  ohservat lonal  evidence  has  been 
accumulating  tor  decades  chat  high-latitude  and 
equatorial  ionospheric  electric  fields  and 
currents  are  coupled  during  active  times.  The 
first  evidence  involved  ground  uia gne t omc t e r 
measurements  t  rom  equatorial  and  hlgh-lat  it.ude 
observator  i»T-  which  showed  variations  closely 
corr*  ‘.tied  In  time  lAkasofu  and  Chapman,  1963, 
1967;  Onvunipohl  1  i  et  al.,  1973).  These  fluc¬ 
tuations  were  often  attributed  to  a  coherent, 
extensive  ionospheric  current  system  generated 
during  point  Mibstorms  with  "it?  turn  currents" 
that  at  fee  ted  the  diyside  equatorial  electrojet. 
However,  there  was  ro  immediate  way  to  rule  out 
the  exp}an.«t  ion  thn:  such  ground  magnetic 
•fleers  were  of  ext i a-ionospher ic  origin,  e.g., 
iron  rino  currents.  Thus  i  causative  polar- 
equ.ifori.il  ionospheric  coupling  covl<i  not  hr 
established  from  magnetic  records  alone.  With 
the  rise  of  ionospheric  sounding  techniques 
such  as  incoherent  scatter  and  VHF  Doppler  radar, 
in  conjunction  with  ground  magnetic  records, 

•  hserv.it  Inna  1  evidence  now  indicates  that  hlgh- 
lat  It. nle  and  equal oiial  magnetic  activity  often 

•r re  late  well  with  each  other  and  with  details 
*1  equatorial  electric,  ti  ’  J  variations.  For 
example,  Matsushita  «aml_  Bnlsjev  |ld77)  found  that 
U>0  *  changes  in  ge* magnet  ic  H  at  equatorial 
Huam.avo  correlated  well  with  radar-measured 
■hi •  roast's  in  westward  electron  drift  velocity 
duiing  a  period  whet'  the  ring  iiitient  was  too 
w ■■  ik  to  iicc.'unt  for  H.  Using  VHF  Doppler  radar, 
carter  el  al.  jlAVij  found  eases  vuero  equa- 
’  *  r  i  a  ;  tin  f  u.if  ii»ns  1  ri  i  ■  m  vt  I i  t  v  v  i  t  h  t  l  me 
ales  of  ..O  '»')  mlie  tes  correlated  with  ge<- 
■a.;r.fi.  M  and  iri  ti.ru  with  tin*  h  i  gh- 1  at  I  l  udr 

*  as  i  n. i L  1  mis  •  t  en  or  the  iH.son  I  ;i.md  magneto- 
c  .ri.  :  i  jet  -  r  ii.  f  1  (f  '*■ i  1 1 :  n!  that  fast 
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reversals  in  the  equatorial  ele> trojet  current 
occurred  both  lay  and  night  during  times  of 
high  Kp,  when  polar  magnetic  activity  is  sta¬ 
tistically  more  frequent.  Thus,  at  least  for 
certain  disturbances,  radar  measurements 
support  the  inference  from  magnetic  records 
that  the  high-latitude  and  equatorial  iono¬ 
spheres  are  coupled  electrically. 

we  will  use*  the  term  "po) ar-equaLorial 
coupling"  to  refer  to  such  coupling  of  iono¬ 
spheric  eltctric  fields  and  currents.  In  our 
approach,  the  coupling  results  from  driving 
forces  in  the  hi gh- l*t  it ude  ionosphere,  and 
the  connection  with  mngnetospheric  dynamics 
is  from  Birkeland  currents.  The  results 
reported  below  essentially  establish  the  role 
of  the  ionospheric  conductivity  in  explaining 
the  observed  polar -equal •  r  la  i  coupling.  The 
major  improvements  In  >ur  results  over  past 
efforts  in  the  area  [Iwasakl  and  Nish Ida,  1967; 
Wolf  ,  1970;  . i  and  ”w off  T 9_7 3 ;  Haeda  _and 
Haekawa,  197  3;  ley  a  t  s  ky  *  t  al  .  ,  1974;  YyuBuhara 
et  al..  197S]  are  in  our  Improved  represeuta- 
! Ions  of  field-aligned  current  distributions 
aria  ;onducr.  iviry,  and  especially  in  our  care¬ 
ful  treatment  of  the  equatorial  conductivity. 
Because  of  these  improvements  we  are  able  to 
determine  the  magnitudes  of  low-  and  equa¬ 
tor  ial-lat  i  t ude  efiects  accurately.  Modeling 
parameters  can  be  better  specified  than  in  the 
past  mainly  because  of  the  synoptic  pictures 
provided  by  the  TRIAD  Satellite  of  the  driving 
field-aligned  currents  under  different  condi¬ 
tions,  and  the  improved  crip  irJ  cal  knowledge  of 
the  ionospheric  co  ’activity  provided  by 
techniques  such  as  incoherent  scatter.  In  the 
following  sections  we  describe  our  models  and 
the  calculated  results. 

The  Models 

'nr  nppt'iach  is  to  solve  the  equation  of  con¬ 
tinuity  fot  the  he ight -integrated  ionospheric 
nii rent  ,  given  the  high- 1st  it ude  distribution  of 
the  driving  ( iel  <1- a  l  l  currents  and  the  iono- 

spherli  layer  c»  mlm  t  i  v  i  t  \  rs .  In  those  calcula¬ 
tions  wo  only  obtain  the  Ionospheric  electric 
fit- his  that  re  suit,  from  magnet  emptier  Ic  coupling, 
and  t  bet  c  fore  do  not  take  in*o  account  other 
dynamical  effects  such  as  neutral  winds.  This 
Is  probably  not  a  serious  deficiency  In  modeling 
the  present  effects,  for  which  the  global  coher¬ 
ence  in  time  occurs  scales  considerably  shorter 
than  the  response  time  of  the  global  wind  system, 
which  is  perhaps  a  «ouple  of  hours  1 Fodder  and 
hanks,  19721. 

The  pi.  run-  of  the  Mold  aligned  currents 
adopted  is  taken  from  the  TRIAD  results  f  li^imn 
and  Futemm,  I'd  761.  Biieflv,  TRIAD  reports  two 
distinct  sots  -  r  large-scale  field-aligned  cur¬ 
rents:  a  poleward  set  (region  1)  which  is  per¬ 

manent,  ami  ;.n  equa f  •' r  i a  1  set  (region  2),  which 
f*«  more  transient.  7he  J  i  s  1 1  l  hut  f  c-ntt  arc  roughly 
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Fig.  1*.  Polar  plot  of  equlpotenc lal  pattern 
(In  kV)  due  to  a  region  1  current  aet  of  10^  A  at 
72“ :  Into  the  Ionosphere  (cross-hatching)  on  the 

dawnaide,  outward  (stippling)  on  the  duakside. 

Noon  is  at  the  top,  dusk  to  the  left,  etc.;  the 
outer  circle  represents  the  equator,  and  the  inner 
circles  are  at  30°  and  60°  latitude. 


symmetrical  northern  to  southern  hemisphere. 
Region  1  currents  flow  into  the  ionosphere  on 
the  daunslde  and  out  on  the  duskside;  region  2 
currents  are  of  opposite  polarity.  The  regions 
have  a  large  spread  in  local  time  and  overlap 
neat  midnight.  During  quiet  times,  region  1 
currents  are  stronger  than  the  region  2  currents 
iwhich  may  even  be  iibsenti  and  amount  to  =  10*5  A. 
During  active  times  the  intensities  of  both  sets 
Increase  and  become  equal  to  about  3  x  10^  A 
flljima  end  Fotemra,  1978).  At  all  times,  evi¬ 
dence  is  that  the  net  current  into  the  ionosphere 
equals  the  net  current  out.  For  our  present 
tiles  1 1  rat  Ion  of  the  TRIAD  synoptic  observations, 
we  placed  the  current  sources  along  arcs  of  con¬ 
stant  latitude  with  a  longitudinal  extent  of  132“ 
and  distributed  the  current  uniformly  along  these 
arcs.  This  structure  is  illustrated  in  Figures 
1  and  2a. 


Fig.  2a.  Equipotential  pattern  of  region  1  cur¬ 
rents  of  10"  A  at  72“  and  region  2  currents  of 
5  x  10^  A  at  66“.  This  scenario  approximates  the 
quiet-time  TRIAD  observations. 


The  direct,  Pedersen,  and  Hall  conductivities 
were  represented  bv  simple  functions  that  repro¬ 
duce  the  observed  major  structures  of  the  iono¬ 
spheric  conductivity,  the  solar-renith  angle 
variation  and  the  auroral  enhancement.  The 
coefficients  were  selected  as  typical  or  roost 
characteristic  values  from  tneasur oment s  of 
Rowe  and  Mathews  [  1973],  Brek_ke  _el_al.  [1974), 
Kennel  sn^  Rees  {1972) ,  and  others.  With  0  the 
colntltude  and  C  the  longitude  measured  eastward 
from  noon,  the  he lght- Integrated  conductivities 
for  the  calculations  presented  here  are  (in  mhos): 

Pedersen:  ct  (8, <>)  •  0.3  +  3.0  exp{-[cos  1 

(sin?  cosif) )  Vl .  804 ) 

+  3.0  expf-l?  -  0.33)^/0.011 

Hall:  cu(0,*)  -  2.0  h  (0,$) 

n  p 

direct:  o,,  (P,*>  -  31.62  a  (0,*) 

0  p 

The  first  two  terms  In  Op  provide  a  low  value  at 
night  and  a  smooth  solar-zenith  angle  variation 
on  the  dayslde  which  closely  approximates  the 
measured  form  of  Van  Zand t  11967];  the  dav-to- 


f*ig.  dh.  K! re  trie  I  IcMs  for  the  -,m  i  et  - 1  j  me 
model  (.same  format  as  Figure  lb>. 
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Fig.  >.  Klee  trie  fields  for  the  vase  ot  equal 
region  l  and  region  currents  (10^  Al  .  The  for¬ 
mat  i:.  the  name  as  in  Figure  lb  and  2b.  Compari¬ 
son  with  Figure  2a  shows  that  the  equatorial 
electric  field  is  oppositely  directed  over  most 
local  times. 

night  conductivity  contrast  is  about  20  to  1. 

The  last  term  provides  a  broad  auroral  enhance¬ 
ment  zone  of  IS0  width  peaking  at  b9°.  These 
functions  represent  a  compromise  between  the 
"best"  observed  or  phenomenological  models  and 
the  requirements  of  suitable  smoothness  for  the 
numerical  ca 1 culat ional  scheme.  The  above 
expression  for  the  conductivities  were  vised  in 
the  usual  lover  formulas  vet.,  R.shbelh  and 
•  •arilott  [19b4)))  to  obtain  the  :  wo-d  linens  ional 
ot  thin -shell  conductivity  tensor  ,  taking  the 
t.irf  h’.i  magnetic  field  to  be  j  ’  a: .  I  he 
ii-ntimiitv  equat Ion 

.1  ( source!  »  V  ,  •  ,  •  ‘  ,  M 

is  solved  numerically  using  difference  formulas 
n  a  Hi  v  10  grid.  The  resulting  algebraic 
i  ;uat  ions  art'  solved  tor  the  ionospheric  poten¬ 
tial  *  row  vhhh  the  electric  field  and 
v  i. rrent  density  art  determined.  Our  model  has 
north-south  svtnmetrv  so  that  ca ! *  ul at  ions  are 
•t'.ade  til  the  northern  hemisphere  with  the  bound¬ 
ary  condition  ot  zero  cut  rent  classing  the 
i  4uator  (Wolf,  19  2  >].  North-south  symmetry 
i  liminates  th.e  cons  ide t  at  ton  ot  coupl  ing 
ofween  conjugate  hemispheres,  e.g.,  along 
n! nsCMwphet  lc  field  lines.  burner  leal  accu¬ 
racy  ot  the  solutions  was  carefully  monitored, 
mJ  an  addtfu’nal  Independent  oust  r  tint  was 
ii”i'«':.ed  tv'  inqiove  the  quality  of  the  results, 
hi.  . ontratnt  is  derived  from  dauss*s  law  in 
'.wo  d  imen*.  ions  and  basically  expresses  the 
’•'cal  requirement  ot  continuity  ot  ionospheric 
rent  that  the  correct  solution  must  satisfy. 

''h  largest  numerical  errors  in  the  results 
•:*-.*nfed  are  a  tew  percent. 

Results  and  PI  sens:,  ion 

d,  i •  1  quiet  and  aitive  von,,  .t  t v  the 

>  1 1  ,-w  i  "  c.  wav.  For  vet  v  juiet  tinea,  the  steady 
.  huat.!  i  t  e  Id- a  l  i  gtved  curvet?:  set  is  given  an 
;ntr::-.rfv  el  1  Pf‘  \  and  fhe  eq  via  t  o  r  we  t  d  -.et  Is 
..•a'  weak*  r  lit  int  *■  a  •,  1 t  \  .  ,'u  l  ivr  t  ires  at** 

. -thing  hot1.  *.et  •.  equal  v\:t>  an  intensity 
.  «■  times  !  dr  \.  it'  t}-.--.,-  ,  a *  a  1  a  t  t  s 


tent  la  Is  are  pulled  poleward  (pushed  equat  or 'war  d ) 
on  the  dusk  (dawn!  side  is  due  to  the  conduct i- 
vitv  contrast  between  day  and  night  (Wolf,  1  *■«  2  C  ]  . 
The  tendency  et  the  subaurotsl  equ  ipot enl  i a ! s  t o 
rotate  toward  later  loeal  fnr.es  was  attributed 
by  Wo  It  [  l  9  ’0]  ami  Vasvl  tunas  ]1970]  to  the 
enhanci'd  mode;  auroral  couduc t i vi ty .  Figure  lb 
shows  the  resulting  elect! le  field  morphology 
at  the  equator  as  a  t unction  of  local  time. 

Note  that  the  elect  tic  Meld  at  the  equat  ot  is 
westward  over  just  a  brief  sector  of  lovul  l  ln«' 
(essentially  the  mi du i ght -dawn  quadrant!.  The 
equatorial  electric  Meld  ts  about  1-2  raV  in. 

This  large  magnitude  suggests  that  impulsive 
activation  or  enhancement  of  field-aligned 
currents  during  disturbed  times  could  contri¬ 
bute  electric  fields  which  are  of  proper 
polarity  and  gioat  enough  strength  at  low 
latitudes  tv’  pro.  ice  significant  fluctuations 
in  (and  even  reversals  ot )  the  quiet-time 
pattern.  Fox  a  quiet-time  equatorial  mor¬ 
phology,  we  regard  the  observational  picture 
o!  Woodman  l  lll  70]  to  be  typical:  eastward 
field  o‘  0.  *>  iiiV’/m  over  most  of  the  dayside, 
with  a  sharp  increase  to  l  mV/m  near  dusk, 
followed  by  a  reversal  to  0 .  *>  mV/m  in  a 
westward  direction  over  most  of  the  nightsidr. 

Figure  2a  gives  the  equ Ipotent ial  pattern 
that  results  from  adding  region  2  currents  of 
S  "  10s  A  at  bbp  tv'  the  current  pattern  of 
Figure  la.  Th.is  nodding  situation  corresponds 
t*'  the  quiet-time  TRIAD  measurements.  Here,  one 
associates  the  region  l  uric-  ls  with  the  rangne- 
tospherlz  convection  electric  field  and  the 
region  2  currents  with  the  shielding  currents 
coming  from  the  Inner  edge  of  the  plasma  street 
v«'.g.t  Wolf  fll)7S]).  Figure  2a  shows  the  pene¬ 
tration  of  the  equipol eut  ial a  to  low  latitudes 
to  be  greatly  restricted,  and  Figure  2b  shows 
cori espond  Ing l v  that  the  electric  fields  at  low 
latitude's  arc  sign  i  J  leant  lv  reduced  by  the 
region  2  source,  tin' ugh  their  directions  are 
more  or  less  the  same  as  in  Figure  lb.  Our 
quiet -lime  results  at  the  equator  resemble  the 
quiet 'tine  morphology  ot  Woodman  119/0].  To 
represent  active  conditions,  the  strength  of 
the  region  2  source  is  Increased  to  match  that 
of  region  1,  namely,  1  a.  Under  these  condi¬ 
tions  the  direction  ot  the  equatorial  electric 
field  is  reversed  compared  to  quiet  conditions 
over  most  local  times,  as  shown  in  Figure  ’. 
it^iraa  and  Poterora  (197b]  associated  the 
strength  of  region  2  currents  with  the  level 
ot  auroral  elcittolet  activity.  The  region  2 
currents  pt**smaablv  originate  In  the  inner 
edge  of  the  plasm.*  sheet /ring  current,  which  is 
often  th*'  source  ot  high- latitude  ionospheric 
activity.  As  suggested  bv  the  TRl/M  measure¬ 
ments,  during  ncfivi  perioJs  current  intensi¬ 
ties  tench  i  '  li'f’  and  the  scales  ot  Figure 
*  are  to  be  multiplied  bv  abt'ul  three.  For 
these  venditions,  the  equatorial  reversal  ot 
the  quiet  -t  1  me  morphology  of  Woodman  ‘19'  P 1 
is  quite  pronoun.  ed .  If  active  periods  at  -■ 

■tode  1  *'d  h\  ■.  lt'ptii  si\v  ,  out -of -phase  tine  fan 
t  ions  ot  tig  lor.  1  tad  no  /  ,  intent  set- 
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In  summary,  the  present  calculations  indi¬ 
cate  the  degree  to  which  the  ionospheric  con¬ 
ductivity  couples  electric  fields  driven  by 
high-latitude  field-aligned  currents  to  the 
equatorial  ionosphere.  The  coupling  is 
sufficient  to  account  for  the  fluctuations  and 
reversals  in  the  equatorial  electric  field 
which  occur  in  active  periods.  Our  model  gives 
<i  consistent  picture  of  magnetically  quiet  and 
active  periods  by  incorporating  essential  fea¬ 
tures  of  the  observed  behavior  of  field-aligned 
currents  to  infer  associated  equatorial  effects. 
Future  applications  will  deal  with  selected 
modeling  conditions  associated  with  several 
classes  of  magnetic  activity. 
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Abstract.  Field-aligned  currents  provide  an  Important  coupling 
between  magnetospherlc  and  ionospheric  processes,  driving  Ionospheric 
electric  fields  on  a  global  scale.  We  have  developed  a  model  In  which 
global  distributions  of  electric  fields  and  currents  are  calculated 
using  mathematical  boundary  conditions  on  field-aligned  currents  that 
are  chosen  to  correspond  to  observations  made  on  the  Triad  satellite. 

The  ionospheric  conductivity  tensor  is  represented  empirically  and 
Includes  spatial  inhomogeneities  resulting  from  the  solar-zenith  angle 
variation  and  a  local  auroral  zone  enhancement.  Numerical  results 
are  accurate  over  the  entire  latitude  range,  pole  to  equator. 

One  application  made  was  to  examine  r '>e  role  of  the  ionospheric 
conductivity  in  coupling  electric  fields  from  high  to  low  latitude. 
Ground  magnetic  observations  and  direct  measurements  of  equatorial 
ionospheric  phenomena  have  indicated  that  magnetospherlc  coupling 
effects  penetrate  at  least  to  L  *  3.2  during  quiet  times  and  to  the 
equator  during  disturbed  times.  Calculations  using  quiet-time  sources 
as  observed  by  Triad  yield  electric  fields  at  all  latitudes  which  are 
at  least  comparable  In  strength  to  those  produced  by  dynamo  wind  models 
and  which  have  non-diurnal  morphologies.  In  our  calculations  for  dis¬ 
turbed  times,  imposed  fluctuations  in  the  strength  of  the  field-aligned 
(Region  2)  currents  are  found  to  produce  large  variations  in  the  mid- 
and  low-latitude  electric  field  and  account  for  reversal  from  its  quiet 
time  configuration,  as  is  often  observed.  The  electric  field  in  the 
polar  cap  is  shown  to  vary  significantly  in  response  to  substorra  condi¬ 
tions  when  Region  2  currents  are  appreciable.  We  conclude  that  under 
both  quiet  and  disturbed  conditions  magnetospher ic  sources  produce 
ionospheric  effects  on  a  global  scale. 

Introduction 

Most  past  efforts  at  modeling  the  coupling  between  the  magnetosphere 
and  ionosphere  have  excluded  proper  consideration  of  the  low-latitude 
ionosphere.  The  importance  of  auroral  zone  effects  on  the  coupling 
has  recently  been  emphasized  with  the  observations  that  have  established 
Klrkeland  currents  to  be  a  permanent  feature  at  high  latitudes.  Impor¬ 
tant  low-latitude  effects  occur  essentially  by  leakage  of  current  from 
high  latitudes.  Thus,  calculated  Ionospheric  results  of  Wolf  [1970], 
Jaggl  and  Wolf  [1973],  Maeda  and  Maekawa  [1973],  Lyatsky  et  al.  [1974], 
Yasuhara  et  al.  [1975],  and  others  are  generally  adequate  at  high  lati¬ 
tudes  but  are  not  designed  to  apply  at  middle  and  low  latitudes. 
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Ionospheric  evidence  Is  now  available  that  suggests  a  direct  rela¬ 
tionship  between  high-  and  low-latitude  phenomena.  Thus,  during  dis¬ 
turbed  times,  rapid  magnetic  variations  closely  correlated  in  time  are 
observed  at  auroral  zone  and  equatorial  stations  (e.g.,  Akasofu  and 
Chapman  [1963];  Onwumechlll  et  al.  [1973]).  Ionospheric  electric 
fields  and  currents  sometimes  fluctuate  and  reverse  in  concert  with 
this  class  of  magnetic  fluctuations  as  shown  by  correlative  studies  of 
magnetograms  and  Tadar  measurements  of  mid-  and  low-latitude  ionospheric 
plasma  drift  episodes  (e.g.,  Balsley  [1966],  Van  Zandt  et  al.  [1971], 
Carter  et  al.  [1976],  Fe.ler  et  al.  [1976]).  and  other  ionospheric  pheno¬ 
mena,  such  as  the  sudden  disappearance  of  equatorial  sporadic-E  and  vhf 
radar  echoes  (Cohen  and  Bowles  [1963],  Rastogl  [1977],  Rastogl  et  al. 
[1977],  Patel  [1978]).  Thus  the  existence  of  a  polar-equatorial  iono¬ 
spheric  coupling  during  disturbed  times  is  documented  observationally. 

At  quiet  times,  Matsushita  [1971]  and  Leont'yev  and  Lyatskly  [1975] 
proposed  that  the  Sq  equivalent  current  system  might  be  driven  in  part 
by  magneto8pheric  sources,  the  conventional  generating  mechanism  being 
the  E-region  dynamo.  Observationally,  F-region  plasma  drift  measure¬ 
ments  at  Millstone  Hill  by  Carpenter  and  Kirchhoff  [1976,  1975]  corrob¬ 
orate  a  magnetospheric  origin  for  the  quiet-time  electric  field  down  to 
L  »  3.2.  More  recently,  Gonzales  et  al.  [1978]  have  shown  that  only  in 
extremely  quiet  times  is  the  classical  dynamo  effect  dominant  over  mag¬ 
netospheric  sources  in  plasma  drifts  seen  at  Millstone  Hill.  It  remains 
an  open  question  how  far  equatorward  the  quiet-time  configuration  of 
magnetospheric  sources  exerts  a  measurable  effect. 

Our  motivation  to  calculate  the  distribution  of  ionospheric  electric 
fields  was  to  evaluate  the  degree  of  coupling  between  the  high-  and 
low-latitude  ionospheres.  As  detailed  in  the  next  section,  our  modeling 
procedure  utilizes  a  realistic  model  of  the  ionospheric  conductivity 
and  the  observed  field-aligned  current  distribution;  careful  numerical 
treatment  of  dip-angle  effects  in  the  equatorial  regions  renders  the 
calculation  accurate  globally.  These  inputs  (conductivity  and  Birke- 
land  current)  can  be  modified  to  conform  to  Improved  empirical  require¬ 
ments  or  special  physical  circumstances.  Model  results  are  presented 
both  for  quiet  and  for  disturbed  periods.  We  find  that  magnetospheric 
sources  contribute  on  an  equal  footing  with  calculated  dynamo  sources 
to  the  quiet-time  electric  field  morphologies  at  all  latitudes;  on  the 
other  hand,  magnetospheric  sources  alone  produce  large  enough  electric 
field  variations  at  middle  and  low  latitudes  to  account  for  observed 
disturbed-time  fluctuations  and  reversals.  Thus  effects  of  magneto¬ 
spheric  coupling  to  the  ionosphere  are  significant  at  all  latitudes 
during  both  quiet  and  disturbed  times. 

Definition  of  the  Model 

Our  approach  is  to  solve  the  equation  of  continuity  of  the  height- 
integrated  ionospheric  current,  given  the  ionospheric  conductivity 
distribution  £(0, 6)  and  the  driving  field-aligned  currents  J(source). 

The  electric  fields  we  calculate  may  not  be  the  total  ionospheric  field; 
f or  example,  we  do  not  include  neutral  wind  effects.  We  impose  boundary 
>  d it  ions  that  require  north-south  symmetry  in  the  conductivities  and 
sources.  If  (0,6)  denote  the  co-latitude  and  azimuth  and  V  is  the 
transverse  gradient,  we  obtain  the  ionospheric  potential  6(0,6)  from 
the  equation  of  continuity 
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J (source)  »  *  £,(Q>$)  '  ly  *(0,9) 

by  numerically  relaxing  its  finite  difference  representation  on  a 
30  x  30  grid  of  nodal  points  covering  the  northern  hemisphere  iono¬ 
sphere.  Southern  hemisphere  results  are  determined  by  the  symmetry 
condition  of  zero  current  crossing  the  equator.  The  numerical  accuracy 
is  verified  by  checking  that  the  calculated  current  is  locally  con¬ 
served  among  adjacent  nodal  points  to  a  good  approximation.  The  calcu¬ 
lated  electric  fields  are  accurate  solutions  to  the  continuity  equation 
(with  assumed  source  and  conductivity  models)  within  a  few  percent 
everywhere  on  the  grid. 

The  ionospheric  height-integrated  conductivities  are  expressed  by 
simple  functions  multiplied  by  adjustable  coefficients.  These  expres¬ 
sions  give  a  good  fit  to  the  principal  structures  observed,  namely, 
a  base  level  constant  conductivity  (denoted  by  subscript  "K")  a  solar- 
zenith  angle  variation  (subscript  "X"),  and  an  auroral  enhancement 
(subscript  "A") : 

Pedersen  :.p(8,<t>)  “  P.pK  +  Epx  exp  {- [cos-1  (sine  cosji) ] 2/ 1 .804 ) 

+  £  exp  f-[0  -  0.35J2/0.01)  mhos, 

I  A 

Hall  ?:H(0 , <p)  «  +  £hx  exp  {- [cos' 1  (sine  cos*) ] 2/ 1 .804} 

+  exp  1-10  -  0.35]2/0.01)  mhos, 
direct  £^(0,41)  «  31.62  T.p(0,4>)  mhos. 

Here,  1)  and  4  (measured  eastward  from  noon)  are  in  radians.  The  coeffi¬ 
cients  are  chosen  to  reproduce  the  most  characteristic  or  typical  values 
Iron  Kennel  and  Kees  (1972],  Rowe  and  Mathews  [1973],  Brekke  et  al. 
[1974],  and  others.  Those  employed  in  the  present  modeling  are  as 
fo 1  lows : 


LPX 

i?A 

UK 

^HX 

rHA 

Quiet 

0.3 

5.0 

3.0 

0.6 

10.0 

6.0 

Dist  urbed 

0.3 

5.0 

6.0 

0.6 

10.0 

18.0 

These  conductivities  are  then  adjusted  for  dip  angle  effects  in  the 
usual  way  (Rlshbeth  and  Garrlott,  1969]  to  yield  the  elements  of  the 
tensor  £(',1)  used  in  the  continuity  equation.  Figure  1  illustrates 
the  variation  in  as  a  function  of  0  along  the  noon-midnight  meri¬ 
dian.  Note  the  different  auroral  enhancements,  peaking  at  69°  latitude, 
for  quiet  and  disturbed  times,  and  the  20  to  1  contrast  between  day  and 
night  values. 

For  driving  currents  we  utilize  the  Triad  results  reported  by  Iljlma 
and  Polemra  [1976,  1978],  shown  In  Figure  2.  These  consist  of  two 
large-scale  current  structures:  Region  1,  often  associated  with 
magnetospher lc  convection,  which  flows  into  the  ionosphere  on  the  dawn- 
side  and  out  on  the  duskside;  Region  2,  often  correlated  with  auroral 
electro  jet  activity,  which  has  the  opposite  polarity.  In  our  idealiza- 
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Fig.  1.  Variation  of  conductivity  tensor 
element  I0e  along  the  noon-midnight  meri¬ 
dian  from  the  dayside  to  the  pole  to  the 
nightslde  equator.  The  dashed  (dotted) 
portions  give  the  model  auroral  conduc¬ 
tivity  values  during  quiet  (disturbed) 
times. 

tions.  Region  1  currents  are  located  at  latitude  72°,  and  Region  2  at 
66°.  We  abstract  the  quiet-  and  disturbed-time  current  strengths  from 
the  Triad  results  to  be: 

Region  1  Region  2 

Quiet  106  A  5  *  105  A 

Disturbed  106  A  106  A 

It  Is  characteristic  of  quiet  times  that  the  Region  1  currents  are 
roughly  twice  as  strong  as  the  Region  2,  whereas  during  more  active 
periods  the  two  sets  carry  equal  currents,  from  106  A  (as  tabulated 
here)  to  several  times  106  a. 

At  present  there  Is  no  definitive  esiplrical  evidence  available  that 
specifies  the  origin  or  the  time-dependent  behavior  of  the  field-aligned 
currents.  Likewise,  theoretical  simulation  efforts  of  the  dynamical 
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Fig.  2.  The  Hold-aligned  current  distributions  determined  using  the 
Triad  satellite  [Iljlma  and  Potemra.  19781 . 

processes  that  couple  the  magnetosphere  and  Ionosphere  are  in  their 
early  stages  of  development.  (See  the  paper  by  Harel  et  al . ,  this 
volume.)  The  Triad  results  adopted  In  our  modeling  scheme  have  only  a 
statistical  significance;  thus  our  calculated  results  would  only  have 
statistical  significance  and  would  not  necessarily  be  representative 
of  pirticular  events  on  a  case  by  case  basis.  Rather,  our  results 
should  characterize  typical  distributions  of  ionospheric  electric  fields 
and  currents  that  are  generated  by  typically  observed  field-aligned 
currents. 


Quiet-Time  Kesui Is 

Oar  main  interest  in  this  first  application  is  to  determine  the 
contribution  made  by  the  field-aligned  currents  to  the  quiet  time  elec¬ 
tric  field  morphology  at  middle  and  low  latitudes.  In  Figure  3  the 
cqui potent ia!  curves  are  given  for  the  quiet-time  modeling  parameters 
presented  above.  It  is  seen  that  a  two-cell  pattern  results  from  the 
quiet  time  pattern  of  Region  1  and  Region  2  currents.  This  agrees  with 
the  two-cell  convection  patterns  of  Heppner  [1977).  The  polar-cap 
Held  is  directed  roughly  dawn  to  dusk  with  magnitude  23  mV/m.  There 
is  a  striking  dusk-dawn  asymmetry,  in  which  duskside  equipotentlals  are 
pushed  toward  higher  latitudes  and  those  on  the  dawnside  are  pulled  to 
lower  latitudes.  This  effect  was  explained  by  Wolf  [1970)  as  due  to 
the  difference  in  conductivity  between  day  and  night.  It  is  interest¬ 
ing  to  note,  since  the  ionospheric  currents  are  likewise  asynmetr ical , 
that  even  for  a  symmetrical  disposition  of  field-aligned  currents  (as 
taken  here)  the  resulting  ground  magnetic  effect  will  be  asymmetrical 
dusk  to  dawn. 

Figure  4  shows  the  meridional  (or  0 )  component  of  the  electric  field 
as  a  function  of  local  time  at  the  latitudes  of  Millstone  Hill  (54°) 
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Fig.  3.  Polar  plot  of  the  equlpotent lal 
curves,  It:  kV,  for  quiet  times.  Here,  the 
Region  1  currents  are  twice  as  strong  as 
the  Region  2,  and  the  auroral  conductivity 
Is  moderate.  Circles  are  centered  at  the 
pole  at  colatitude  30°,  60°,  and  90°  (the 
equator) .  The  total  cross-polar-cap 
potential  drop  is  about  50  kV  In  this 
model. 

and  Areclbo  (30°);  Eg  at  the  equator  la  Identically  aero  In  these 
models.  This  field  component  would  correspond  to  an  eastward  component 
of  drift  of  Ionospheric  plasma.  The  calculated  magnitudes,  5  mV/m  at 
54°  and  1.5  mV/m  at  3D°,  are  approximately  as  large  as  those  generated 
in  neutral  wind  models  (Matsushita  and  Tarpley.  1970]  .  Thus,  magneto- 
spheric  sources  and  neutral  winds  mist  both  be  Included  and  play  com¬ 
parable  roles  in  modeling  quiet  time  effects  at  low  latitudes.  Figure 
4  displays  a  predominantly  diurnal  morphology  at  54°,  where  maximum  and 
minimum  electric  field  values  occur  about  twelve  hours  apart.  In  our 
model,  the  premidnight  field  is  northwards  and  this  would  generate  a 
westward  F-reglon  plasma  drift  at  Millstone  Hill;  this  drift  is 
observed  by  Gonzales  et  al.  (1978]  on  all  but  the  quietest  days.  The 
classical  dynamo  effect,  which  gives  southward  fields  in  this  sector, 
would  be  dominant  over  magnetospher lc  driving  effects,  according  to 
our  results,  only  during  very  quiet  times. 

In  Figure  5  the  azimuthal  (or  f,  positive  eastwards)  qulct-tlme 
electric  fields  are  given  at  54°,  30°,  and  the  equator,  corresponding 
to  northward  plasma  drifts.  Again  the  calculated  magnitudes  are  com- 
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Fig.  4.  Meridional  component  of  the 
electric  field,  in  mV/m,  at  54°  and 
30°  latitude,  as  a  function  of  local 
time,  for  the  quiet-time  model  (Q)  . 

Positive  Eg  corresponds  to  a  south¬ 
ward  directed  electric  field. 

parable  to  those  from  neutral  wind  models.  A  new  result  is  that  this 
component  has  strong  non-diurnal  Fourier  components,  with  principal 
maxima  and  minima  occurring  about  six  hours  apart.  This  results  as  a 
beating  effect  between  the  diurnal  field-aligned  source  structure  and 
the  diurnal  but  out-of-phase  conductivity  distribution.  To  see  the 
magnitude  of  the  non-diurnal  components,  we  idealize  the  morphologies 
of  Figure  5  and  show  the  Fourier  amplitude  spectrum  of  this  signal  in 
Figure  6.  We  find  that  the  12  hour  component  dominates  over  the 
diurnal,  and  that  the  8  hour  component  is  also  quite  strong.  Thus  the 
subdiurnal  harmonics  in  observed  quiet-time  plasma  drift  morphologies 
(e.g.,  Blanc  et  al.  11977])  may  have  significant  contributions  from 
magnetospheric  sources  (as  well  as  from  the  atmospheric  dynamo) . 

Disturbed-Time  Results 

Our  goal  in  modeling  active  times  is  to  investigate  the  fluctuations 
and  reversals  observed  in  the  otherwise  quiet-time  electric  field  at 
low  latitudes.  In  order  to  do  this  it  is  necessary  to  Impose  a  time- 
dependence  in  our  modeling  parameters.  According  to  observations, 

Regl  on  2  currents  tend  to  be  variable  and  roughly  equal  in  intensity 
lo  the  relatively  fixed  Region  1  currents  during  disturbed  times 
llijtma  and  Potemra.  1976,  1978].  In  addition,  the  auroral  conduc¬ 
tivities  Increase,  due  to  enhanced  particle  precipitation,  during  times 
when  the  auroral  electrojets  are  active.  Thus  we  adopt  the  basic 
disturbed-time  model  parameters  as  tabulated  above  in  section  11. 
Fluctuations  are  considered  to  be  alternating  sequences  of  the  quiet 
and  disturbed  model-results. 

The  disturbed  time  equipotential  pattern  Is  given  in  Figure  7,  The 
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Fig.  5.  Azimuthal  component  of  the 
electric  field  (mV/m)  at  54°,  30°, 
and  0°,  as  a  function  of  local  time, 
for  quiet  times  (Q) .  Positive 
corresponds  to  an  eastward  directed 
field. 


enhanced  Region  2  source  has  moved  the  formerly  prominent  duskside 
vortex  to  the  noon  sector  and  reduced  it  in  size.  The  equipotent lal 
curves  tend  to  be  distorted  poleward  (equatorward)  on  the  dusk  (dawn) 
side,  as  in  the  quiet-time  model.  An  interesting  feature  of  this  case 
is  the  sunward-directed  electric  field  over  the  polar  cap.  The  magni¬ 
tude  and  intensity  of  the  polar  cap  electric  field  is  a  sensitive 
function  of  characteristics  of  the  Region  2  current  system.  The  sunward 
direction  of  the  electric  field  shown  in  Figure  7  results  because  in 
this  case  the  Rec.lon  1  and  Region  2  currents  have  equal  intensities. 

The  direction  would  rotate  towards  dusk  if  the  Region  2  Intensity  were 
decreased,  and  towards  dawn  if  Increased.  The  polar  cap  electric  field 
is  weak  (8  mV/m  her  because  of  the  opposite  polarity  and  equal 
strength  of  the  Region  1  and  Region  2  currents.  Polar  cap  convection  in 
this  case  is  directed  dawn-to-dusk.  Such  deviations  from  quiet-time 
antisunward  convection  are  occasionally  seen  during  disturbed  periods 
(e.g.,  Ileppner  (19771). 

In  Figure  8  a  comparison  is  made  of  the  azimuthal  electric  field 
components  at  30°,  between  quiet-  and  disturbed-time  results,  curves 
Q  and  D’.  Over  most  of  local  times,  the  strengthened  Region  2  currents 
have  reversed  the  electric  field  at  this  latitude.  (The  added  auroral 
conductivity  during  disturbed  times  has  little  effect  upon  the  reversal, 
as  can  be  seen  from  curve  D,  which  uses  the  sources  of  D'  but  the  con¬ 
ductivity  of  Q.)  The  equatorial  electric  field  is  similar  to  that  at 
30°,  but  of  roughly  half  the  magnitude,  as  shown  in  Figure  9.  Since 
the  characteristic  Triad  active  field-aligned  current  strengths  often 
tend  to  be  several  times  the  value  (106  A)  used  in  obtaining  Figures  7 
through  9,  it  follows  that  the  equatorial  electric  field  can  be  several 
times  tliat  shown  In  these  figures.  Thus  strengths  up  to  2  mV/m  can 
sometimes  be  obtained  and  these  will,  for  most  local  times,  be  directed 
oppositely  to  the  quiet-time  electric  fields.  Hence,  variations  In  the 
Region  2  currents  may  be  responsible  for  the  fluctuations  and  reversals 
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Fig.  8.  Azimuthal  electric  field  com¬ 
ponents  at  30°  for  quiet  (Q)  and  dis¬ 
turbed  (D1)  models.  (Model  D  utilizes 
sources  of  D'  and  conductivities  of 
Q.)  Fluctuations  in  the  strength  of  the 
Region  2  currents  and  auroral  precipi¬ 
tation  would  cause  this  electric  field 
component  to  vary  between  Q  and  D'. 

in  the  equatorial  electric  field  and  associated  ionospheric  phenomena, 
as  referenced  earlier.  The  observed  rapid  variations  in  the  equatorial 
field  cannot  be  accounted  for  by  the  neutral  wind  system  with  its  large 
inertia  and  slower  response  time.  Further  discussion  of  this  case  is 
given  in  Nopper  and  Carovillano  [1978]. 

Summary  and  Conclusions 

The  principal  outcome  of  our  modeling  study  is  the  calculation  of 
relatively  large  electric  fields  due  to  magnetospheric  sources  at  all 
ionospheric  latitudes.  During  quiet  times,  the  observed  configuration 
of  field-aligned  currents  contributes  large  magnitudes  and  non-diurnal 
spectral  terms  to  mid-  and  low-latitude  electric  fields  and,  conse¬ 
quently,  to  plasma  drifts.  These  effects  of  magnetospheric  origin 
should  not  be  excluded  from  models  of  the  quiet  ionosphere  since  they 
compete  with  the  traditional  source  of  quiet-time  variations,  dynamo 
winds.  Adopting  a  plausible  model  for  time  variations  in  the  Region  2 
currents  leads  to  an  explanation  of  the  fluctuations  and  reversals  of 
the  low-latitude  electric  field  observed  in  disturbed  times.  Region  2 
currents  also  affect  the  direction  and  magnitude  of  polar  cap  electric 
fields  in  an  important  way.  In  summary,  the  dynamics  of  the  magneto¬ 
sphere/ionosphere  coupled  system,  as  manifest  in  field-aligned  currents. 


Fig.  9.  Same  as  Figure  8  for  the 
equator . 


NOITEK  AND  I  AK0V1LLAN0 


667 


exerts  a  global  Influence  on  the  ionosphere  for  all  but  the  quietest 
levels  of  geomagnetic  activity. 
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On  the  Orientation  of  the  Polar  Cap  Electric  Field 

R.  W.  Nopper,  Jr.,  and  R.  L.  Carovillano 

Department  of  Physics,  Boston  College,  Chestnut  Hill.  Massachusetts  07/47 

Using  realistic  models  of  the  ionospheric  conductivity  and  the  field-aligned  currents,  we  have  deter¬ 
mined  how  the  distribution  of  the  electric  field  in  the  polar  cap  ionosphere  is  controlled  by  the  day-night 
contrast  in  conductivity  and  by  the  relative  strengths  of  the  region  I  and  region  2  field-aligned  currents 
The  sunward  directed  conductivity  gradient  acts  to  set  up  a  space  charge  in  the  polar  cap  which  crowds 
the  equipotentials  toward  the  dawn  sector  for  current  sources  of  both  region  I  and  region  2  polarity,  this 
effectively  shifts  the  polar  cap  convection  pattern  toward  dawn  Our  results  show  further  that  for  a  given 
conduc.ivity  distribution  the  orientation  of  the  electric  held  in  the  central  polar  cap  depends  sensitively 
on  the  relative  strengths  of  the  Birkeland  current  pairs:  for  very  weak  region  2  currents  (quiet  limes)  the 
polar  electric  Held  is  directed  =60°  east  of  noon,  for  equal  region  I  and  2  currents  (disturbed  times)  the 
direction  is  - 10"  east  of  noon,  and  for  stronger  region  2  than  I  currents  (which  may  happen  on  occasion) 
the  electric  Held  points  into  the  prenoon  sector.  These  findings  imply  that  the  orientation  of  the  polar  cap 
electric  Held  should  serve  as  a  measure  or  index  of  the  ratio  of  region  I  to  region  2  net  current  intensity 
and  should  possess  correlations  with  geomagnetic  activity  similar  to  those  of  this  ratio.  This  analysis  does 
not  include  effects  of  possible  source  currents  in  the  region  of  the  polat  cusp 


I  INTRODUCTION 

Although  the  electric  field  in  the  polar  cap  is  often  de¬ 
scribed  as  being  directed  dawn  to  dusk,  it  is  well  known  that 
only  rarely  is  this  precisely  the  case  In  fact,  both  the  electric 
field  magnitude  and  the  electric  field  direction  across  the  polar 
cap  are  known  to  change  in  response  to  variations  in  a  num¬ 
ber  of  effects,  including  the  season  of  the  year  and  the  values 
ol  various  components  of  the  IMF  (interplanetary  magnetic 
field)  (eg.  I  .angel,  1974:  llrppner.  1977;  Fairfield,  1977)  Such 
observations  indicate  that  the  patterns  of  plasma  convection 
(or  electrostatic  cqinpoiential  lines)  frequently  differ  from 
uniform  anltsunward  flow. 

In  recent  years  the  Triad.  Isis  2,  and  S3-2  spacecraft  have 
measured  vector  magnetic  field  fluctuations  that  indicate  the 
presence  of  held -aligned  currents,  and  a  synoptic  pattern  for 
these  currents  has  been  defined  |/i/imn  and  Poiemra.  1978).  In 
addition  atmospheric  models  are  constantly  being  improved 
\Jacchia.  |977|,  and  the  tv.  'mque  of  incoherent  scatter  radar 
has  been  useful  tn  obtaining  measurements  of  ionospheric 
conductivities,  thus  allowing  fajrly  realistic  conductivity  mod¬ 
els  to  be  developed  T (eating  the  observed  field-aligned  cur¬ 
rents  as  driving  terms  and  given  the  ionospheric  conductivity 
distribution,  one  can  calculate  the  distributions  of  ionospheric 
electric  licit!  and  current  which  would  result  in  that  model  sit¬ 
uation  Recent  efforts  of  this  sort  have  been  made  by  Ya.su 
horn  ct  al  ||975|.  Maekawa  and  Maeda  11978],  Ntshvt  el  a  I 
|  I97K|,  Sapper  and  Carovillano  (1978,  I979|,  and  Kamtde  and 
Matsushita  |I979]  The  Iasi  three  references  utilized  con¬ 
ductivity  models  that  lepresent  the  principal  structure  ob¬ 
served  (variation  with  solar  zenith  angle  and  auroral  enhance¬ 
ment)  and  are  intended  to  apply  globally 

I  he  purpose  of  this  brief  note  is  to  discuss  properties  of  the 
electric  field  in  the  polar  cap.  We  note  that  even  for  a  symmet- 
inal  disposition  of  field-aligned  currents,  an  idealized  model 
situation,  the  electric  field  at  the  pole  is  usually  not  directed 
dawn  to  dusk  In  addition,  our  calculations  show  that  the  di- 
icciion  of  the  polar  cap  electric  field  is  a  sensitive  function  of 
the  relative  strengths  of  the  field-aligned  current  sheets  (re¬ 
gion  I  and  2  current  pairs)  Since  the  relative  strengths  of 


these  input  currents  are  known  lo  vary  with  geomagnetic  ac¬ 
tivity,  one  would  expect  the  polar  cap  electric  field  direction 
to  show  a  similar  dependence 

2.  Tut  Modeling  Scheme 

We  have  developed  a  numerical  procedure  to  solve  the 
equation  for  the  continuity  of  ionospheric  current, 

/(source)  -  V, •£(<?.  it >).Vjff>(0,  <>) 

using  finite  difference  methods  and  an  iterative  technique. 
Here.  <1>  is  the  ionospheric  electrostatic  potential,  from  which 
the  ionospheric  electric  field  and  current  distributions  are  ob¬ 
tained.  The  assumed  representations  for  the  height-integrated 
tensor  conductivity  £  and  for  the  driving  field-aligned  cur¬ 
rents  ./(source)  are  discussed  briefly  below.  We  ignorr  dynamo 
winds,  which  are  expected  lo  play  a  secondary  role  in  the 
high-latitude  region  of  interest  here.  The  accuracy  of  numeri¬ 
cal  results  is  checked  by  applying  a  test  based  upon  Causs's 
law  in  Iwo  dimensions  lo  various  areas  of  the  numerical  grid. 
This  check  certifies  that  ionospheric  current  is  locally  con¬ 
served  everywhere  on  the  grid  to  within  a  few  percent. 

Our  empirical  conductivity  model  is  given  by  simple  func¬ 
tions  of  position  multiplied  by  adjustable  coefficients.  The 
functions  are  designed  to  give  reasonable  approximations  to 
the  solar  zenith  angle  and  auroral  variations  of  conductivity, 
and  the  coefficients  are  selected  to  give  agreement  to  pub¬ 
lished  observed  values  |e.g.,  Brekke  et  at ,  1974;  Rowe  and 
Mathews.  1973)  With  0  the  colatitude  and  <>  the  azimuth  mea¬ 
sured  eastward  from  noon  the  conductivity  model  employed 
in  the  present  study  is,  in  mhos, 

Pedersen 

03  ♦  VO  exp  j  -|cos  '(sin  (>cos$))!/l.804) 

+  3  0  exp  {-(*-0.35)70.0l| 


Hall 


Direct 
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Fig  |.  The  variations  in  tensor  conductivity  elements  Em, 
and  along  the  noon-midnight  meridian  (i.e..  from  the  noon  equa¬ 
tor  to  the  pole  to  the  midnight  equator). 

These  expressions  are  used  in  the  usual  formulas  |e.g.,  Kish- 
belli  and  Harriott.  1969]  for  the  elements  of  the  layer  con¬ 
ductivity  tensor  X(8,  <>)  appearing  in  the  continuity  equation, 
f  igure  I  shows  the  variation  of  these  elements  along  the 
meridian  from  the  equator  at  noon  to  the  pole  to  the  equator 
at  midnight.  Azimuthally  symmetrical  auroral  zone  enhance¬ 
ments  peaked  at  69°  latitude  are  superposed  on  the  back¬ 
ground  conductivity,  which  has  a  20: 1  day-to-night  contrast. 
This  conductivity  model  is  symmetrical  between  northern  and 
southern  hemispheres  We  calculate  northern  hemisphere  re¬ 
sults  using  the  boundary  condition  of  zero  meridional  current 
at  the  equator. 

The  driving  currents  used  are  based  upon  the  statistical 
1  riad  results  reported  by  lijima  and  Potemra  [1978],  These  re¬ 
sults  indicate  the  existence  of  two  distinct  sets  of  large-scale 


4  A 


Fig.  2 a 

f  ig  2  Polar  plou  of  the  cquipotential  patterns  for  the  con¬ 
ductivity  model  and  field-aligned  currents  given  in  section  2.  In  each 
figure  the  same  five  equipotennal  curves  (labeled  in  kilovolts)  are 
shown  I  he  outer  circle  is  the  equator,  the  inner  ones  are  30*  and  60° 
latitude  The  ratios  of  region  I  and  region  2  current  strength  art  (<r)  2, 
i hi  I  V  (r)  I.  (if)  0  67.  and  (r)0.s  In  all  cases  the  region  I  current  was 
Ilf  A  1  he  resulting  polar  cap  electric  field  orientations  are  about  (a) 
♦  'll".  (6|  40°.  (c)  12".  Id)  -411°.  and  (e)  -72°.  measured  counter¬ 
clockwise  from  noon 


field-aligned  current:  region  I  currents  flow  into  the  iono¬ 
sphere  at  the  high-latitude  edge  of  the  auroral  oval  on  the 
dawnside  and  out  of  the  ionosphere  on  the  duskstde.  region  2 
currents  have  the  opposite  polarity  and  are  located  somewhat 
cquatorward  of  the  region  I  set.  Both  regions  have  a  large 
spread  in  azimuth  and  overlap  near  midnight.  Characteristi¬ 
cally,  in  weakly  disturbed  periods  the  region  I  set  carries 
about  10*  A,  and  the  region  2  half  to  three  quarters  of  that 
value.  (In  quiet  periods  the  region  2  set  is  weaker  and  may 
even  be  absent.)  During  active  periods,  however,  the  strengths 
of  both  sets  increase  to  perhaps  3  x  10*  A,  and  the  region  2  set 
may  even  exceed  the  region  I  set  in  net  amperage.  For  our 
present  purposes  we  have  idealized  the  above  as  follows.  The 
region  1  currents  arc  applied  in  two  bands,  extending  azimu¬ 
thally  from  18°  to  150°  (duskstde)  and  from  210°  to  342° 
(dawnside)  at  72°  latitude:  the  region  2  currents  are  in  bands 
with  the  same  azimuthal  extents  but  at  66°  latitude.  Since  one 
purpose  of  this  note  is  to  illustrate  the  effect  of  changing  the 
ratio  of  region  I  and  region  2  current  strength,  the  region  I 
source  is  held  fixed  at  10*  A.  and  the  region  2  source  is  varied 
to  give  region  l/region  2  ratios  of  2 : 1,  1.5 : 1,  1:1,  1 : 1.5,  and 
1 : 2.  The  last  two  ratios  are  included  here  for  pedagogical  pur¬ 
poses  and  probably  do  not  occur  in  practice,  since  the  current 
intensity  in  region  2  is  rarely  reported  to  exceed  that  in  region 
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Fig  2d 

I.  The  effects  of  input  currents  in  the  region  of  the  polar  cusp 
are  not  included  in  this  analysis. 

3.  Discission  ano  Results 

Before  considering  our  calculated  results  it  is  useful  to  re¬ 
view  polar  cap  equipotential  patterns  under  simpler  circum¬ 
stances.  In  the  presence  of  only  symmetrically  located  region 
1  sources,  if  the  Pedersen  conductivity  were  uniform  and  the 
Mall  conductivity  absent  then  the  equipotentials  would  pos¬ 
sess  dawn-io-dusk  symmetry,  and  the  electric  field  would  es¬ 
sentially  point  from  dawn  to  dusk.  The  effect  of  a  nonzero 
Hall  conductivity  is  demonstrated  in  the  axisymmetrical 
model  of  Vasyliunas  1 1 970],  this  model  has  a  driving  source 
that  is  sinusoidal  in  azimuth,  neglects  day-night  variations  in 
conductivity,  and  includes  an  enhanced  auroral  conductivity 
band.  For  this  case,  polar  cap  equipotentials  remain  uniform 
and  sun  aligned,  although  auroral  and  subauroral  equipoten- 
tials  are  rotated  because  of  the  band  of  enhanced  auroral  con¬ 
ductivity.  Thus  such  axisymmetrical  modeling  variations  have 
no  effect  on  the  orientation  of  the  electric  field  or  the  distribu¬ 
tion  of  equipotentials  in  the  polar  cap. 

The  condition  of  axisymmetry  was  broken  by  Wolf  |I970], 
who  included  a  realistic  variation  with  solar  zenith  angle  in 
the  conductivity  tensor;  however,  his  model  was  not  applied  to 
the  polar  cap.  Atkinson  and  Hutchison  J 1978}  have  recently 
presented  analytical  results  for  Pedersen  and  Halt  con- 
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i 
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ductivities  that  decrease  exponentially  across  the  polar  cap. 
They  find  that  the  steeper  the  day-night  conductivity  decrease 
front  noon  to  midnight  is,  the  more  the  equipotential  lines  are 
squeezed  into  the  dawn  sector  of  the  polar  cap.  The  corre¬ 
sponding  effect  was  essentially  seen  and  explained  earlier  by 
Wolf  [  I970J  in  his  analysis  of  mid-latitude  equipotentials.  In 
the  polar  cap  the  combined  effect  of  the  polarity  of  the  Birke- 
land  currents  fin  at  dawn,  out  at  dusk,  or  vice  versa)  and  the 
direction  of  the  conductivity  gradient  produces  an  electric 
field  mote  intense  at  dawn  than  at  dusk.  Equipotential  lines 
are  therefore  more  dense  near  dawn  than  at  dusk.  This  effect 
arises  physically  because  of  the  charge  density  (proportional  to 
-V  2,  •  I,  where  I  is  the  total  current)  generated  by  the  day- 
night  conductivity  gradient.  This  charge  density  is  negative 
for  region  I  polarity  and  positive  for  region  2  polarity,  be¬ 
cause  of  the  reversed  direction  of  the  Hall  current,  and  in¬ 
creases  in  magnitude  toward  the  dayside.  Thus  for  either 
polarity  the  density  of  the  equipotential  lines  is  greater  in  the 
dawn  sector  than  in  the  dusk  sector.  Of  course,  the  polar  cap 
electric  fields  generated  by  region  I  and  region  2  polarities  are 
oppositely  directed  and,  for  example,  give  rise  to  oppositely 
directed  convection  patterns. 

For  a  fixed  source  the  distortion  of  the  equipotentials  and 
the  direction  of  the  electric  field  are  controlled  by  the  con¬ 
ductivity  gradient  This  feature  can  be  seen  in  the  results  of 
Atkinson  and  Hutchison  1 1978]  for  symmetrical  sources,  where 
the  polar  cap  electric  field  direction  always  differs  from  dawn 
to  dusk  and  is  further  altered  by  the  steepness  of  the  con¬ 
ductivity  gradient.  The  effect  of  the  conductivity  gradient  can 
also  be  seen  in  our  earlier  work  (Hopper  and  Carovillano, 
1978,  I979J  and  in  the  calculations  of  Nisbet  el  al.  (1978)  and 
Kamidr  and  Matsushita  (1979)  In  our  present  calculations  we 
show,  in  addition,  how  the  direction  of  the  polar  cap  electric 
field  is  affected  by  the  relative  intensity  of  the  region  I  and  re¬ 
gion  2  current  sources. 

In  Figures  2 a  2e  we  show  the  global  patterns  of  the  same 
five  equipotential  lines  for  the  ratios  2,  I  S,  1.  2/3,  and  1/2  of 
region  I  to  region  2  current,  as  defined  in  section  2.  (Region  l 
currents  are  fixed  at  10*  A  for  these  calculations.) 

Figure  2 a  would  correspond  roughly  to  a  quiet  time  field- 
aligned  current  configuration.  The  mid-latitude  equipotential 
pattern  follows  the  ‘midlatitude  rules’  given  by  Wolf  [  1970); 
the  duskside  (dawnside)  equipotentials  are  distorted  toward 
(away  from)  the  pole  by  the  action  of  the  day-night  con¬ 
ductivity  asymmetry.  The  polar  cap  equipotentials  are  more 
dense  in  the  dawn  sector,  as  discussed  above.  The  direction  of 
electric  field  at  the  pole  is  about  +50°  (i.e.,  east  of  noon). 

The  region  2  source  intensity  increases  in  strength  in  the  se¬ 
quence  of  Figures  2 h-e  These  figures  show  the  negative  vor¬ 
tex  to  be  decreasing  in  size  and  rotating  about  180°  across  the 
dayside  toward  dawn,  while  the  positive  equipotentials  rotate 
in  the  same  sense  (clockwise)  through  midnight  toward  dusk, 
eventually  to  form  a  large  vortex  there.  (This  vortex  does  not 
appear  in  Figure  2r  because  of  our  choice  of  displayed 
equipotentials.)  In  all  cases  the  expected  behavior  of  the 
equipotential  shapes,  given  above,  is  seen.  In  the  sequence  of 
cases  in  Figure  2  the  direction  of  the  polar  cap  electric  field  is 
seen  to  rotate  clockwise  as  the  region  2  source  intensity  is  in¬ 
creased.  The  direction  goes  from  about  +50°  for  the  ratio  2;  I 
of  region  I  to  region  2  current  intensities  to  about  -70°  for 
the  ratio  1:2.  To  see  the  extremes  for  the  cases  represented, 
the  polar  cap  electric  field  direction  is  +60°  when  region  2 
currents  are  absent  and  -96°  when  region  I  currents  are  ab- 
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sent.  The  lack  of  symmetry  between  these  extreme  values  is 
because  the  region  1  and  2  currents  are  applied  at  different 
latitudes.  As  was  discussed  earlier,  these  extreme  values  de¬ 
pend  upon  parameters  in  the  conductivity  model.  (We  will 
publish  a  further  study  on  this  issue.) 

In  summary,  we  have  discussed  how  the  direction  of  the  po¬ 
lar  cap  electric  field  is  affected  by  the  day-night  conductivity 
gradient  and  by  the  relative  strengths  of  region  I  and  region  2 
currents.  The  effect  of  the  conductivity  gradient  is  to  increase 
the  density  of  etjuiftotentials  (and  hence  shiA  the  polar  cap 
convection  pattern)  in  the  dawn  sector,  as  shown  in  the  work 
of  several  authors  [Atkinson  and  Hutchison,  1978;  Hopper  and 
Carovillano,  1978,  .1979;  Kamide  and  Matsushita,  1979;  Hisbet 
et  al  ,  1978].  Our  principal  purpose  is  to  display  the  depen¬ 
dence  of  the  direction  of  the  polar  cap  electric  field  on  the  ra¬ 
tio  of  region  I  to  region  2  current  intensities.  This  ratio  is  ob¬ 
served  to  vary  from  slightly  greater  than  2  to  less  than  unity. 
In  our  calculations  the  direction  of  the  polar  cap  electric  field 
is  found  to  depend  sensitively  on  this  ratio.  Conversely,  this 
leads  us  to  suggest  that  measurement  of  the  orientation  of  the 
polar  cap  electric  field  could  serve  as  an  index  of  the  relative 
strength  of  region  1  /region  2  currents  and  should  possess  cor¬ 
relations  with  geomagnetic  activity  similar  to  those  of  the  ra¬ 
tio.  Our  results  do  not  include  effects  of  currents  in  the  polar 
cusp  region. 
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E.  Refinements  of  the  Model 

As  indicated  by  the  foregoing  discussions,  global  ionospheric  modeling 
is  a  worthy  means  to  gain  understanding  of  many  physical  effects  and  offers 
a  wide  range  of  possibilities  for  comparing  theory  and  observation.  To  advance 
the  utility  of  our  model,  a  number  of  refinements  are  needed  both  in  representing 
our  input  parameters  in  a  more  realistic  fashion  and  in  addressing  a  new  class 
of  physical  effects.  Here  we  shall  briefly  discuss  possible  improvements  to 
the  global  modeling  scheme. 

The  conductivity  model  used  is  adequate  for  general  studies  but  can  be 
refined  to  accord  with  the  latest  improved  ionospheric  models  (e.g.,  Mehta,  1979), 
and  particularly  with  auroral  region  incoherent  scatter  results  (e.g.,  Vondrak 
and  Vickrey,  1979;  de  la  Beaujardlere  et  al.,  1977)  which  provide  a  far  more 
realistic  picture  of  high-latitude  conductivities  for  both  quiet  and  active  times. 
In  order  to  maintain  numerical  accuracy,  the  high-latitude  portion  of  the  mesh 
will  need  to  be  refined,  as  observation  shows  the  highly  conducting  auroral 
regions  to  possess  abrupt  gradients.  Similar  numerical  techniques  used  for 
our  model  can  be  applied  to  the  modeling  of  individual  auroral  arcs.  Given 
the  field-aligned  currents  over  an  arc  and  the  ionospheric  conductivity  in 
that  region,  the  distributions  of  electric  field  and  current  (electrojets) 
can  be  calculated.  The  parameters  involved  in  such  modeling  can  be  determined 
in  case  studies  from  the  ground  and/or  from  rockets,  so  that  a  tight  control  on 
the  calculations  may  be  available.  In  addition  to  upgrading  the  auroral  region 
conductivity,  it  would  be  useful  to  include  a  seasonal  dependence  in  the  model, 
in  which  the  latitude  of  the  subsolar  point  (where  the  conductivity  is  largest) 
changes  with  time  of  year.  This  seasonal  effect  would  give  rise  to  very  large 
effects  in  the  electric  fields  and  currents  at  the  equator  and  polar  caps.  For 
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example,  the  summer  polar  cap  is  much  more  conducting  than  the  winter  polar 
cap,  and  during  equinox  there  is  a  strong  gradient  across  the  polar  cap  (Paper 
III).  Thus,  quantities  which  depend  significantly  upon  polar  cap  conductivities, 
such  as  electric  fields,  convection  patterns  (equipotentials) ,  and  currents, 
can  be  expected  to  depend  sensitively  upon  time  of  year. 

In  addition  to  the  two  major  pairs  of  field-aligned  currents  modeled  at 
present,  there  is  another  set  observed  in  the  cusp  region  (Iljima  and  Potemra, 
1976)  which  seem  to  depend  upon  the  y-component  (By)  of  the  interplanetary 
magnetic  field  (IMF)  (Iij ima  et  al . ,  1978;  McDiarmid  et  al.,  1977,  1978a,  b, 

1979;  Wilhjelm  et  al . ,  1978).  This  is  of  importance  because  evidence  has  been 
accumulating  that  the  interplanetary  magnetic  field  affects  the  global  ionosphere 
and  hence  could  profitably  be  introduced  into  our  modeling.  For  example, 

Maezawa  (1976)  and  Burke  et  al.  (1979)  have  postulated  a  correlation  between 
the  high  latitude  convection  pattern  and  IMF  Bz;  Heppner  (1972,  1977),  Friis- 
Christensen  and  Wilhjelm  (1975),  Wilhjelm  et  al.  (1978),  and  McDiarmid  et  al. 
(1978b,  1979)  have  studied  the  relation  of  high  latitude  convection  and  IMF  By. 

In  middle  latitudes,  Matsushita  (1975)  has  shown  the  dependence  of  the  Sq 
current  system  focus  on  IMF  By.  At  low  latitudes,  studies  relating  ionospheric 
electrodynamic  effects  and  IMF  Bz  have  been  published  by  Patel  (1978), 

Rastogi  (1977),  Rastogi  and  Chandra  (1974),  Rastogi  and  Kroehl  (1978),  Rastogi 
and  Patel  (1975),  and  Rastogi  et  al.  (1978);  correlation  between  ionsopheric 
effects  and  By  have  been  shown  by  Bhargava  and  Rangarajan  (1977)  and  Galperin 
et  al. (1978) .  A  schematic  summary  of  the  salient  observations  for  the  northern 
high  latitude  region  is  given  in  Figure  4.3.  Unfortunately,  the  inclusion  of 
cusp  currents  to  introduce  an  IMF-dependent  effect  into  the  modeling  would  greatly 
enlarge  the  computational  task.  These  currents  are  anti-symmetrical  (northern- 
to-southern  hemisphere)  whereas  the  Region  1  and  2  current  sources  used  In  <>ur 
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present  model  are  symmetrical.  With  symmetrical  sources  we  are  able  to 
impose  a  boundary  condition  at  the  equator  that  greatly  improves  numerical 
accuracy  and  reduces  the  size  of  the  calculational  effort.  With  cusp-currents 
included,  both  hemispheres  will  need  to  be  explicitly  modeled  and  solutions 
obtained  over  the  complete  globe. 

For  a  more  complete  picture  of  the  ionospheric  electric  fields  and 
currents  during  quiet  and  disturbed  times,  it  would  be  necessary  to  include 
neutral  winds  in  the  dynamo  (E)  region  of  the  ionosphere.  For  wind  field  u, 
magnetic  field  JB,  and  conductivity  this  source  term  has  the  form  V^£(uxB)  . 
Strictly  speaking,  the  neutral  wind  speed  u  to  be  used  in  our  formulation 
should  be  height-independent.  Unfortunately,  models  of  neutral  winds  do  not 
provide  such  modes  for  us  to  adopt  (Forbes  and  Lindzen,  1976);  one  exception 
is  the  (1,-2)  mode,  which  fortuitously  is  the  main  contributor  in  quiet  times 
(J.  Forbes,  pvt.  communication).  At  the  present  time,  there  is  actually  no 
comprehensive  dynamo  wind  model  that  has  general  acceptance.  This  problem 
limits  the  possibility  of  including  neutral  winds  in  our  modeling  scheme  with 
a  high  degree  of  confidence,  at  least  for  the  immediate  future. 

F.  Magnetic  Ground  Signatures 

The  culmination  of  our  model  would  be  to  relate  predictions  to  ground 
magnetic  signatures.  Difficulties  to  be  overcome  to  accomplish  this  task  are 
discussed  in  this  section.  The  problem  is  complicated  by  many  factors,  but 
significant  advances  are  certainly  possible  if  the  proper  time  and  effort  is 
devoted  to  the  task.  Ground  magnetic  signatures  are  treated  more  specifically 
in  chapter  5  within  a  well  defined  context. 

There  is  a  prodigious  amount  of  information  available  to  the  space  physics 
community  in  the  form  of  magnetograms  from  thousands  of  station-years.  Magnetic 
observatories  are  relatively  inexpensive  to  instrument  and  operate  and,  since 
recent  efforts  at  installing  high-quality  instruments  in  chains  and/or  networks 
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have  been  exerted,  magnetic  records  will  continue  to  form  the  core  of  our  ground 


data  base.  Two  ways  in  which  magne tog rams  have  been  used  for  many  years  are: 

(1)  the  records  from  one  or  a  small  number  of  stations  are  inspected  to  infer 
limited  information  about  nearby  current  systems  (e.g.,  the  use  of  auroral  zone 
magnetograms  to  investigate  electrojets);  (2)  the  records  from  a  global  collection 
of  stations  are  processed  into  "equivalent  current  systems"  (see  Stern  (1977) 
for  a  discussion  of  the  procedure).  Recently,  the  equivalent  current  system 
concept  has  been  greatly  upgraded  by  the  use  of  magnetometer  arrays,  which  by 
viitue  of  a  relatively  dense  distribution  of  instruments  on  a  regular  grid 
allows  more  information  to  be  extracted  from  magnetic  variation  data  than 
otherwise  possible.  Even  with  dense  global  coverage  on  the  ground,  however, 
it  is  not  possible  to  infer  the  current  system  uniquely  that  is  responsible 
for  a  given  magnetic  disturbance  distributed  over  the  globe.  Special  assumptions 
are  needed  to  obtain  a  current  system  that  corresponds  to  a  given  magnetic 
disturbance.  Thus,  to  derive  "equivalent  current  systems,"  the  ansatz  is  made 
that  all  currents  flow  in  the  ionosphere  (which  is  demonstrably  false). 

The  task  of  inferring  currents  from  a  given  magnetic  disturbance  is  an 
"inverse  problem"  to  the  usual  task  of  calculating  the  magnetic  disturbance 
that  results  from  a  given  set  of  currents  in  space.  For  specified  boundary 
conditions,  the  magnetic  field  resulting  from  a  given  current  configuration  is 
unique.  In  contrast,  the  inverse  problem  is  not  unique.  For  example,  Fukushima 
(1969)  considered  the  problem  of  a  vertical  line  current  along  the  z-axis 
entering  a  uniformly  conducting  plane  located  at  right  angles  to  the  z-axis. 

This  current  con f i gura t ion  -  vertical  line  current  plus  distributed  currents 
in  the  horizontal  plane  at  z=0  -  produces  no  magnetic  field  below  the  conducting 
plane  (i.e.,  at  z<0).  Thus,  if  the  ionosphere  is  approximated  by  a  plane 
conductor,  even  at  times  of  zero  magnetic  disturbance  at  the  surface  of  the 
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earth  there  may  or  may  not  be  overhead  currents  of  magnetospherlc  origin.  This 
simple  example  demonstrates  that  inferences  to  be  made  from  ground  magnetic 
observations  of  ionospheric  and/or  magnetospheric  currents  can  only  be 
made  within  the  context  of  a  model.  In  the  "equivalent  current’'  representation, 
the  ground  magnetic  signature  is  assumed  to  be  generated  by  currents  confined 
completely  to  the  ionosphere. 

The  development  of  realistic  models  to  represent  the  magnetic  signatures 
of  ground  observations  is  further  complicated  by  induction  effects  within 
the  earth.  These  induction  effects  "contaminate"  the  magnetic  record  in 
reflecting  the  external  current  system.  For  example,  Nopper  and  Hermance 
(1974)  have  shown  that,  for  substorm-like  events,  induction  effects  alter 
magnetic  component  disturbance  magnitudes  and  introduce  significant  phase 
shifts.  Early  methods  of  quantifying  effects  of  induction  were  the  construction 
of  "induction  arrows"  and  the  separation  of  magnetic  fluctuations  into 
internal  and  external  parts  using  analytic  techniques  (e.g..  Price,  1967). 

With  the  refinement  of  mathematical  techniques  and  the  introduction  of  arrays, 
more  quantitative  characterization  of  source  and  induced  fields  are  now 
available  (e.g.,  Schmucker,  1970a,  b;  Alabi  et  al . ,  1975).  Gregori  and 
Lanzerotti  (1980)  have  reviewed  this  subject.  Most  of  these  studies  have 
been  aimed  at  determining  earth  conductivity  structure  and  not  ionospheric 
current  configurations.  The  determination  and  removal  of  the  effects  of 
induction  from  ground  magnetic  observatory  records  constitutes  a  major 
field  of  research  in  itself. 

Other  recent  considerations  of  the  inverse  problem  are  those  of  Kamide  et 
al .  (1976)  and  Oldenburg  (1978).  Kamide  et  al.  (1976)  have  developed  an 
automated  technique  to  solve  Poisson's  equation  for  the  magnetic  potential 
by  a  finite-difference  method  to  obtain  an  equivalent  current  system. 


Oldenberg  (1978)  used  a  linear  Inverse  theory  to  construct  a  most  probable 
current  system  compatible  with  the  ground  disturbance  pattern.  Both  of  these 
efforts  are,  of  course,  subject  to  the  difficulties  mentioned  above,  namely, 
the  necessity  to  constrain  the  currents  to  flow  in  a  particular  configuration 
in  order  to  obtain  results. 

The  calculation  of  a  ground  magnetic  signature  from  a  given  current 
distribution  may  be  performed  without  the  inherent  ambiguity  of  the  inverse 
problem.  Thus,  magnetic  field  signatures  have  been  obtained  from  wire 
models  (Cummings ,  1966),  from  direct  integration  over  distributed  currents 
using  the  Biot-Savart  law  (Crooker  and  Siscoe,  1974;  Kisabeth ,  1972),  from 
wire-model  simulations  (Olson,  1974)  and  from  a  magnetic  charge  representation 
(Kisabeth,  1979);  these  computations  are  generally  straightforward  but  tedious. 
The  difficulty  with  these  direct  approaches  derives  from  our  lack  of  knowledge 
of  the  actual  spatial  configurations  of  the  contributing  currents,  which 
are  known  to  include  ionospheric,  field-aligned,  ring,  tail,  magnetopause, 
and  earth  components.  An  unfortunate  practical  difficulty  of  assumed  current 
models  is  that  several  of  these  components  of  a  complete  current  system  in 
a  given  event-model  contribute  on  a  roughly  equal  footing,  and  some  in  opposing 
senses  (e.g.,  Fukushima  and  Kamide,  1974;  R.  Wolf,  pvt.  communication; 

Alabi  et  al . ,  1975).  As  a  result,  the  calculated  ground  signature  is  quite 
sensitive  to  details  of  the  model  current  system,  whose  correspondence  with 
reality  is  poorly  known.  These  limitations  are  apparent,  in  context,  in  the 
simple  model  described  in  Chapter  5  of  the  ground  signature  for  the  Birkeland 
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5.  SIMPLE  ESTIMATE  OP  THE  GROUND  MAGNETIC  SIGNATURE  OF  THE 
BIRKELAND- IONOSPHERIC  CURRENT  SYSTEM  AT  LOW  LATITUDES 

A.  Introduction 

Inherent  difficulties  for  the  determination  of  the  ground  magnetic 
signature  for  the  complete  current  configuration  generated  by  Birkeland 
currents  were  discuseed  in  Chapter  4.  In  this  chapter,  a  simple  model  is 
developed  to  obtain  first  estimates  of  the  ground  magnetic  signature  of  the 
Birkeland-Ionospheric  current  system.  The  approach  taken  was  to  combine 
the  Fukushima-Vasyliunas  theorem  and  the  requirement  of  induced  shielding 
currents  within  the  earth  with  a  schematic  representation  of  the  main  Birkeland 
current  systems  to  reduce  the  problem  to  an  analytic  expression  of  great 
simplicity.  With  this  model,  we  show  that  the  main  characteristic  of  the  low 
latitude  disturbance  field  of  such  a  current  system  is  a  harmonic  variation 
in  the  H-component  of  the  field.  The  local  time  of  maximum  disturbance  depends 
on  the  nature  of  the  ionospheric  closure  of  the  Birkeland  current  system.  One 
suggestion  compatible  with  our  model  is  the  possibility  that  the  classical 
mid  latitude  disturbance  field  asymmetry  normally  associated  with  a  partial 
ring  current  is  instead  the  signature  of  the  Birkeland  current  system. 

The  systems  of  permanent  polar  Birkeland  current  first  identified  by 
Zmuda  and  Armstrong  (1974)  are  now  believed  to  be  the  main  electrodynamical 
links  connecting  the  ionosphere  to  the  magnetosphere  and  the  solar  wind.  The 
near-earth  distribution  of  these  current  systems  and  their  variations  in  time 
and  space  are  presently  the  subjects  of  intensive  investigation  (see  the  review 
by  Potemra,  1979).  Already  a  fairly  complete  picture  of  the  current  patterns 
has  been  inferred  statistically  from  data  on  encounters  with  the  current  systems 
made  by  low-altitude,  polar-orbiting  spacecraft  (e.g.,  Iljima  and  Potemra,  1976; 


lb/B).  Those-  patterns  wore  utilized  in  the  modelling  calculations  presented  in 
Chapter  4.  This  relatively  recent  emergence  of  a  new  and  fairly  well  defined 
global  current  pattern  comes  after  a  long  history  of  analysis  and  interpretation 
of  ground-based  magnetic  records.  The  question  thus  arises:  What  is  the 
ground  magnetic  signature  of  the  Birkeland  current  systems?  The  signature 
most  likely  would  be  distinctive  and  be  present  among  the  known  types  of 
magnetic  disturbances,  where  its  interpretation  would  remain  incomplete  or 
incorrect  until  it  is  identified. 

A  strong  motivation  for  identifying  the  ground  magnetic  signature  of 
the  Birkeland  current  systems,  for  our  purposes,  is  that  by  monitoring  the 
signature,  the  behavior  of  the  current  systems  themselves  can  be  inferred. 

Ground  based  magnetometer  observations  have  the  advantage  of  high  time  resolu¬ 
tion  (essentially  continuous)  and  easy  availability  compared  to  satellite 
observations.  For  these  reasons,  we  have  under  the  contract  taken  a  number  of 
steps  toward  the  goal  of  achieving  ground  based  magnetometer  characterization 
of  the  polar  Birkeland  current  systems.  Our  approach  was  to  develop  theoretical 
models  to  determine  the  types  of  ground  magnetic  signature  produced  by  given 
Birkeland  current  configurations.  Our  findings  are  presented  in  this  chapter. 

B.  Ionospheric  Shielding  Effects 

The  magnetic  signature  at  the  surface  of  earth  is  created  by  the  combined 
effects  of  internal  and  external  currents.  Internal  currents  account  for  the 
main  field  (essentially  dipolar)  and  induction  effects  due  to  transient 
disturbances  of  external  origin.  External  currents  are  primarily  of 
magnetospheric  and  ionospheric  origin  and  constitute  complicated  three  dimen¬ 
sional  arrays.  The  full  Birkeland  current  system  (with  its  full  three 
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dimensional  character)  is  one  such  external  current  system.  Under  certain 
conditions,  it  is  well  known  that  the  ionosphere  serves  to  shield  the  earth 
from  external  magnetic  effects.  The  physical  basis  for  the  shielding  is 
nicely  demonstrated  in  the  simple  theorems  stated  recently  by  Fukushima  (1969) 
(planar  geometry)  and  Vasyliunas  (1970)  (spherical  geometry).  We  begin  with 
a  statement  and  simple  derivation  of  these  theorems. 

Given:  An  infinite  plane  is  connected  to  infinity  by  straight  lines 
perpendicular  to  its  surface.  The  lines  are  all  in  the  same  half  space  (e.g., 
all  extending  upwards  from  a  horizontal  plane).  Along  these  lines,  electrical 
current  flows  to  and  from  infinity  and  closes  in  the  plane,  which  has  uniform, 
isotropic  conductivity.  There  is  no  restriction  on  the  number  of  lines  or 
their  location. 

Fukushima 's  Theorem:  There  is  no  magnetic  field  produced  by  this 
current  system  in  the  current -free  half  space  (i.e, ,  below  the 
horizontal  plane) . 

Given:  A  spherical  surface  with  uniform,  isotropic  electrical  conductivity 
is  connected  to  infinity  by  straight,  radial  lines.  The  lines  carry  electrical 
current  to  and  from  the  sphere,  and  the  current  closes  in  the  surface  of  the 
sphere.  There  is  no  restriction  on  the  number  or  the  disposition  of  the  lines. 
Vasyliuna^  Theorem:  There  is  no  magnetic  field  produced  by  this 
current  system  in  the  interior  of  the  sphere. 

It  is  understood  that  these  theorems  are  stated  in  the  context  of  magneto¬ 
statics  in  the  sense  that  signal  propagation  times  can  be  ignored.  We  prove 
the  theorems  for  one  connecting,  current-carrying  line  in  each  case.  Then, 
since  the  problem  is  linear,  the  theorem  follows  for  an  arbitrary  number  and 
disposition  of  lines  by  linear  superposition. 


•  »*  <Vf 
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In  the  case  of  a  single  line,  we  may  choose  the  line  to  be  the  axis 
of  a  cylindrical  polar  coordinate  system.  Since  the  only  source  for  the 
field  in  this  case  is  the  current  through  the  line,  the  fields  must  be  axi- 
symmetric  with  respect  to  the  line.  Also,  since  V_«B_  =  0,  axial  symmetry 
requires  that  the  magnetic  field  be  of  the  form 

B  *  B(r,z,t)<£_ 

where  r,z  and  4>  are  the  radial,  axial,  and  azimuthal  variables  of  the  chosen 
cylindrical  coordinate  system,  and  |  is  a  unit  vector  in  the  azimuthal  direction. 
The  magnetic  field  is  to  be  determined  from  Ampere's  law  given  in  MKS  notation 
by: 

V  x  B  =  P0J  |f  (1) 

Application  of  Stoke' s  Theorem  to  equation  (1)  yields 

B(r ,  z ,  t)  =  2^  (Ic  +  ID)  (2) 

where  1^  and  Ip  are  the  total  electric  and  displacement  currents  through  a 
circle  of  radius  r  centered  at  the  location  z  on  the  axis. 

Through  any  point,  P,  in  the  current-free  half  space  or  in  the  interior 
of  the  sphere  we  can  draw  an  axis-centered  circle.  Because  of  axisymmetry, 
the  circle  will  lie  wholly  in  the  current-free  region  in  both  cases.  Thus, 

1=0  for  any  such  circle  for  each  point  P. 

To  find  the  Ip  appropriate  to  the  circle  through  P,  we  need  to  know  the 
flux  of  electric  field  through  the  circle.  The  relevant  electric  field  for 
the  problem  is  generated  by  the  electrical  charges  that  are  carried  by  the 
current  between  infinity  and  the  conducting  surface  in  each  case,  Because  of 


the  symmetries  involved,  planar  in  one  case  and  spherical  in  the  other,  the 
electric  field  lies  wholly  between  the  surfaces  and  the  charges  at  infinity. 

We  must  keep  in  mind  that  the  charges  transferred  to  the  surfaces  are  compen¬ 
sated  for  by  equal  and  opposite  charges  at  infinity.  Thus  the  field  is 
confined  between  the  surfaces  and  infinity  in  the  manner  of  a  parallel  plate 
condenser  or  a  spherical  condenser.  Hence,  we  have  IQ  =  0  also. 

Since  we  have  shown  that  Ic  and  Ip  are  zero,  it  follows  that  B(r,z,t)  =  0 
for  any  point  in  the  current-free  half  plane  and  in  the  interior  of  a  conducting 
sphere.  Hence,  is  zero  everywhere  in  these  source  free  regions,  and  stated 
in  the  theorems. 

In  most  intended  applications,  more  than  one  current-carrying  line  would 
be  involved,  and  current  could  then  flow  in  both  directions.  In  this  way  the 
net  current  can  be  maintained  equal  to  zero.  Also  the  result  is  unaffected 
if  the  number  of  current-carrying  lines  is  proliferated  to  generate  distributed 
current  systems, 

C.  Model  Results  at  Low  Latitudes 

The  application  of  the  theorem  to  the  case  of  Birkeland  current  systems 
is  fairly  direct.  The  Birkeland  field-aligned  currents  would  correspond  to 
the  radial  currents  and  the  ionosphere  to  the  conducting,  spherical  surface. 
Field  lines  at  high  latitudes  are  approximately  radial.  For  example,  at  a 
geomagnetic  colatitude  of  20°,  the  field  is  inclined  10°  with  respect  to  the 
radial  direction,  the  main  contribution  to  the  surface  magnetic  field  produced 
by  a  field-aligned  current  is  made  by  the  near-earth,  quasi-radial  portion 
of  the  field  line  (Crooker  and  Siscoe,  1974).  Thus,  the  Birkeland  field- 
aligned  current  systems  approximate  the  radial  currents  for  which  the  theorem 


strictly  applies. 
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With  sufficient  accuracy  for  some  purposes  the  Birkeland  current  systems 
can  be  represented  as  four  separate  radial  line  currents  directed  toward  or 
away  from  the  ionosphere,  as  illustrated  in  Figure  5.1.  At  not  a  very  great 
depth  in  the  earth,  the  conductivity  is  large  enough  that  earth  shielding 
currents  are  induced  and  have  the  effect  of  cancelling  the  Birkeland  current 
field  at  greater  depths,  but  increasing  the  field  at  the  surface.  It  is 
easy  to  demonstrate  that  the  field  at  the  surface  produced  by  the  Birkeland 
currents  and  their  shielding  currents  is  approximately  the  same  as  the  field 
produced  by  extending  the  Birkeland  currents  through  the  earth,  as  suggested 
by  the  dashed  lines  on  the  figure.  (One  should,  in  fact,  also  retain  a 
current-free  gap  between  the  level  at  which  the  shielding  currents  flow  and 
the  ionosphere.  But  this  is  a  small  correction  at  low  latitudes.) 

The  largest  perturbation  in  the  H-component  at  low  latitudes  produced 
by  the  Birkeland  currents  systems  occurs  along  the  noon-midnight  meridian. 

The  value  at  point  A  on  the  noon  equator  is 

u 

Ha  =  “  sin  “i  ~  l2  sin  a2^  (3) 

e 

where  is  measured  positive  northward,  I^(i  =  1,2)  is  the  total  current  in 
the  region  i  current  system,  and  is  the  colatitude  of  the  region  i  current 
system.  By  symmetry: 


H 


B 


(4) 


The  ionospheric  closure  currents  provide  shielding  and  tend  to  reduce  the 


perturbation  magnetic  field  at  the  surface.  As  we  have  seen,  if  the  conduc¬ 
tivity  of  the  ionosphere  were  isotropic  and  uniform,  the  closure  current 
would  cancel  the  field  of  the  Birkeland  current  everywhere  under  the  ionosphere. 


The  conductivity  of  the  ionosphere  is  neither  isotropic  nor  uniform. 

The  presence  of  Hall  currents  and  of  spatial  gradients  in  the  conductivities 
violate  the  conditions  necessary  for  perfect  cancellation  of  the  Birkeland 
current  field  at  the  surface  of  the  earth.  There  are  two  limiting  cases  which 
illustrate  the  range  of  the  surface  fields  that  result.  In  the  direct  closure 
limit,  the  Birkeland  currents  close  in  narrow,  latitudinally  confined  electro¬ 
jets.  In  the  other  limit,  the  Pederson  conductivity  is  taken  to  be  uniform, 
so  that  the  field  of  the  Birkeland  current  plus  its  Pederson  current  closure 
is  zero  under  the  ionosphere.  However,  the  associated  Hall  current  produces 
a  non-zero  field  at  the  surface. 

The  direct  closure  case  is  represented  to  a  good  approximation  by  equations 
(3)  and  (4)  for  and  Hg.  The  field  produced  at  low  latitudes  by  the  direct 
closure  current  (subscript  D)  is  considerably  smaller  than  that  produced  by 
the  Birkeland  current  (Crooker  and  Siscoe,  1974).  Thus, 

(Vd  =  -  (Vd  *  -  rtT  (I1  sin  “1  '  T2  sin  “2)  (5) 

e 

Also,  by  symmetry,  referring  to  points  C  and  D  in  Figure  5.1, 

<Vo  ■  <Vd  ■  0  <6> 

To  obtain  similar  expressions  for  the  case  of  uniform  Pederson 
conductivity,  we  assume  that  the  ratio  of  the  Hall  to  Pederson  conductivities 
is  a  constant  n,  i.e.,  the  Hall  conductivity  is  also  taken  to  be  uniform. 
(However,  the  sign  of  n  must  change  across  the  equator  to  reflect  the  change 
in  the  sign  of  the  vertical  component  of  the  earth's  magnetic  field.)  The 
Hall  current,  J„,  is  perpendicular  to  the  Pederson  current  J  everywhere, 

“n 

in  a  direction  given  by  J  x  £  in  the  northern  hemisphere,  where  r  is  a  unit 
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vector  in  the  radial  direction.  J.,  flows  in  the  -J  x  r  direction  in  the 
southern  hemisphere.  The  Hall  current  streamline  pattern  is,  therefore, 
everywhere  orthogonal  to  the  streamline  pattern  of  the  Pederson  current.  It 

is  easy  to  demonstrate  by  means  of  a  sketch  of  the  current  streamlines  that  at 

middle  and  low  latitudes,  the  Hall  current  pattern  can  be  obtained  simply  by 

performing  a  90°  eastward  rotation  in  longitude  of  the  Pederson  current  pattern. 

At  high  latitudes,  the  two  patterns  deviate  even  after  rotation,  since  one 
closes  completely  in  the  ionosphere  whereas  the  other  terminates  at  the  iono¬ 
spheric  intersection  of  the  Birkeland  currents.  The  two  patterns  would  coincide 
exactly  after  rotation  in  the  limit  in  which  and  approach  zero  (c.f,, 

Figure  5.1). 

The  field  at  middle  and  low  latitudes  is  determined  mainly  by  the  middle 
and  low  latitude  part  of  the  current  pattern.  Hence,  the  Hall  current  field  is 
approximately  the  same  as  the  field  of  the  Pederson  current  appropriately 
rotated.  But  from  the  Fukushima-Vasyl iunas  theorem,  we  know  that  the  surface 
field  of  the  Pederson  current  is  equal  and  opposite  to  that  of  the  Birkeland 
current.  After  tracing  the  connection  from  the  Hall  current  field  to  the 
Pederson  and  Birkeland  current  fields,  remembering  that  the  Birkeland  and 
Pederson  fields  have  opposite  signs,  we  find  that  the  Hall  field  is  approximately 
the  same  as  the  Birkeland  field  obtained  by  rotating  the  current  pattern  in 
Figure  5.1  90°  westward.  We  must  also  multiply  the  Birkeland  field  by  the 
conductivity  ratio  n  to  obtain  the  strength  of  the  Hall  field.  Thus,  we  arrive 
at  the  results  for  uniform  conductivity  (subscript  u) : 


dusk-dawn : 


(8) 


In  both  limits,  the  quantity  (1^  sin  -  I2  s^n  a 2>  enters  and  must  be 
supplied  empirically.  From  the  review  article  by  Potemra  (1979),  we  extract 
the  statistical  values  shown  in  Table  1  for  quiet  and  disturbed  conditions. 

Table  5.1  Parameters  used  to  evaluate  antipodal  magnetic 
disturbances  at  the  equator. 


Disturbance 

Level 

Z*l 

I* 

2 

*1 

a2 

(Vd 

(HC>u 

Quiet  ( |AL| 

<  100  y) 

1.55 

1.10 

16° 

19° 

-4.3  y 

8.7 

Disturbed  ( 

|  AL |  >  100  y) 

2.70 

2.45 

18° 

21° 

2.7  Y 

-5.4 

Difference 
Antipodal  Difference 

7.0  y 
14.1  y 

-14.1 

-28.2 

*  in  units  of  106  amperes 
+  n  *  2  assumed 

Since  the  base  line  for  measuring  the  amplitude  of  disturbances  is 
established  during  quiet  conditions,  the  proper  measure  of  the  disturbance 
produced  by  the  Birkeland  current  systems  is  the  difference  between  quiet 
and  disturbed  conditions.  The  total  amplitude  of  the  global  disturbance  field 
is  the  difference  between  antipodal  points  on  the  equator. 

If  the  disturbance  field  is  represented  by  a  sine  wave  with  zero  average, 
the  Table  5.1  gives  for  its  amplitude  a  value  of  7  y  in  the  direct  closure  case 
and  14  y  in  the  uniform  conductivity  case.  The  Hall-to-Pederson  conductivity 
ratio  n,  was  taken  to  be  2  (Nopper  and  Carovillano,  1978),  The  low  latitude 
disturbance  field  in  the  direct  closure  case  takes  the  form  of  a  noon-midnight 
asymmetry  in  H,  with  a  positive  increment  at  noon.  In  the  uniform  conductivity 
case,  it  takes  the  form  of  a  dawn-dusk  asymmetry  in  H  with  a  negative  increment 
at  dusk. 

Returning  now  to  the  original  question  of  the  low  latitude  ground  magnetic 
signature  of  the  Birkeland  currents  system,  we  see  that  it  should  take  the  form 
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of  an  antipodal  asymmetry  in  the  H-component  of  the  disturbance  field,  but 
that  the  local  time  meridian  of  its  maximum  amplitude  depends  on  the  detailed 
nature  of  the  associated  ionospheric  currt=nt  systems.  The  most  notable 
representative  of  an  antipodal  asymmetry  at  low  latitudes  is  the  well-known 
dawn-dusk  magnetic  storm  disturbance  field  asymmetry  (Sugiura  and  Chapman, 

1960;  Crooker  and  Siscoe,  1971),  which  since  1966  has  been  attributed  to  a 
partial  ring  current  system  (Cummings,  1966;  Crooker  and  McPherron,  1972; 

Crooker  and  Siscoe,  1974).  This  disturbance  field  is  characterized  by  a 
negative  increment  in  H  in  the  dusk  meridian,  and  thus  conforms  to  the  case  of 
uniform  conductivity.  The  amplitude  found  in  our  model  is  based  on  averages 
of  numbers  representing  quiet  and  disturbed  conditions.  Individual  cases 
would  undoubtedly  yield  both  larger  and  smaller  amplitudes.  However,  it  is 
true  that  in  previous  studies  of  the  low  latitude  disturbance-field  asymmetry, 
an  amplitude  of  14  y  would  be  among  the  very  lowest  values  considered  (e.g. , 
Crooker  and  Siscoe,  1971).  The  calculation  of  the  asymmetry  for  values  of  Ij 
and  I 2  based  on  time  periods  for  which  | AL | >  500  y  or  1000  y,  would  likely 
represent  more  faithfully  the  conditions  under  which  the  low  latitude  disturbance 
field  asymmetry  is  most  pronounced. 

We  note  with  regard  to  the  amplitude  of  the  low  latitude  disturbance  field 
produced  by  the  Birkeland  current  system,  that  the  disturbance  fields  of  the 
two  components  of  this  system,  the  region  I  and  region  2  currents,  nearly 
cancel  each  other.  What  remains  is  a  small  residual  field.  For  example,  the 
amplitude  of  the  region  2  current  system  alone  in  the  uniform  conductivity, 
disturbed  condition  case  is  =  -95  y.  It  would  appear  as  if  the  two 

component  systems  are  coupled  In  such  a  way  that  if  they  act  in  direct  opposition 
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as  we  have  arranged  them  to  do,  they  nearly  cancel  each  other.  In  the  interest 
of  simplicity,  we  have  perhaps  been  too  strict  in  imposing  exact  opposition. 

The  results  of  this  study  were  reported  at  the  Geomagnetism  Workshop 
held  at  AFGL,  April  1979.  The  presentation  lead  R.A.  Wolf  and  coworkers  to  use 
the  Rice  University  magnetospheric  convection  model  to  test  the  general  hypo¬ 
thesis  that  the  low  latitude  disturbance  field  asymmetry  results  from  the 
Birkeland  current  systems  (Harel  et  al . ,  1979).  They  indeed  found  the  predicted 
asymmetry,  but  they  judged  that  the  main  factor  contributing  to  it  was  a  down¬ 
ward  current  into  the  polar  ionosphere  in  the  noon  meridian  and  an  upward 
current  in  the  midnight  meridian.  The  low  latitude  disturbance  field  was 
then  produced  in  the  direct  closure  mode  from  this  noon-midnight  current 
residual.  (These  currents  would  not  be  represented  in  our  models  as  presented 
in  Chapers  4  and  5.) 

A  noon-midnight  Birkeland  current  component  can  be  simulated  in  our 
two  component  model  by  means  of  a  longitudinal  offset  between  the  region  1 
and  region  2  current  systems.  The  region  2  system  must  be  rotated  westward 
with  respect  to  the  region  1  system. 

In  terms  of  the  quantities  already  defined,  if  we  add  an  offset  angle  B, 
then  the  direct  closure,  disturbance  field  resulting  from  an  uncompensated 
westward  rotation  of  the  region  2  current  system  is 

p 

Hc  S  I2  sin  a2  *ln  (10) 

e 

Using  disturbance  condition  parameters  and  B  =  45°,  we  find  (H„)  =  -34  v. 

C  o 

Once  again,  larger  values  would  result  if  currents  appropriate  to  more  dis¬ 
turbed  conditions  were  used. 
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D.  Reconmended  Approach 

For  the  reasons  indicated  in  Chapter  4,  while  magnetic  signatures 
corresponding  to  our  globally  determined  ionospheric  current  configurations 
could  easily  be  calculated,  the  utility  of  such  results  could  be  limited. 

Above  we  presented  within  the  context  of  a  simple  model  most  features  necessary 
to  consider  in  determining  the  ground  magnetic  signature.  For  the  purpose 
of  determination  of  such  magnetic  signatures,  the  approach  we  recommend  is 
to  build  upon  the  current  Report  to:  (1)  increase  the  applicability  of  our 
models  to  real  events  (see  section  4E) ;  (2)  develop  a  magnetospheric  extension 
of  our  current  models  to  provide  closure  of  currents  out  of  the  ionosphere 
(this  should  be  done  in  recognition  of  the  work  of  R.  Wolf  and  collaborators) ; 
(3)  take  earth  induction  effects  into  account  realistically. 
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6.  AURORAL  BOUNDARY  LOCATIONS 

Magnetospher tc  conditions  affect  both  the  occurrence  and  location 
of  auroras,  but  until  recently,  there  was  no  means  of  moniLoring  auroral 
activity  on  a  global  scale.  Images  from  DMST  satellites  now  provide  the 
opportunity  to  observe  auroras  every  1-2/3  hours  over  each  pole.  The 
images  sacrifice  time  resolution  compared  with  all  sky  camera  coverage, 
but  this  is  more  than  compensated  for  by  the  images'  wide  spatial  coverage, 
clockwork  regularity,  and  timely  reception  vital  for  practical  applications. 
With  proper  training,  DMSP  images  can  be  interpreted  reliably  without  very 
much  difficulty.  Also,  the  sensitivity  of  DMSP  Images  rivals  or  surpasses 
the  sensitivity  of  all  sky  cameras,  under  good  conditions  (Eather ,  1979). 

These  features  make  the  images  useful  tools  for  observing  auroral  activity 
and  determining  boundary  positions. 

The  auroral  boundary  marks  the  termination  of  particle  precipitation 
sufficient  to  excite  aurora  above  a  threshold  intensity,  which  depends  on 
the  means  of  detection.  Auroras  generally  expand  both  poleward  and  equator- 
ward  with  increasing  magnetic  activity,  but  the  connection  is  not  tight; 
although  particle  precipitation  may  carry  some,  as  yet  unknown,  fraction  of 
field  aligned  currents,  various  definitions  of  magnetic  activity  may  also 
reflect  electric  currents  linked  only  indirectly  to  auroras,  such  as  the  ring 
current  or  certain  ionospheric  currents. 

We  have  studied  the  relationship  between  the  equatorward  auroral  boundary 
and  two  familiar  indices,  Kp  and  AE  (S)ieehan  and  Carovillano,  1978).  A 
reprint  ol  the  article  is  included  at  the  end  of  this  chapter.  Rased  on  162 
DMSP  images  taken  over  the  north  pole,  the  Kp  dependence  of  the  boundary  near 
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local  midnight  is  consistent  with  earlier  findings  by  Lui  (1975),  who  used 
ISIS  2  satellite  images,  and  results  from  the  south  pole  DMSP  images  analyzed 
by  Nagata  et  al.  (1975).  In  Figure  6.1,  the  average  boundary  latitude  is 
seen  to  steadily  decrease  from  67°  CGL  (Kp<0)  to  63°  CGL  (Kp^4)  with  increasing 
magnetic  activity.  A  2°  CGL  range  of  scatter  above  and  below  the  average 
locations  can  be  attributed  partly  to  the  incompatibility  between  instantaneous 


Figure  6.1.  Average  corrected  geomagnetic  latitude  of  the  equatorward  auroral 
boundary  at  0100  CGLT  as  a  function  of  Kp.  Error  bars  demark 
one  standard  deviation.  From  Sheehan  and  Caro vi llano  (1978). 


auroral  determinations  and  the  3  hour  interval  defining  Kp,  but  similar 
scatter  also  applies  to  the  boundary  latitude  with  AE,  which  is  a  virtually 
instantaneous  index.  The  average  boundary  location  plotted  in  Figure  6.1  follow 
rather  well  the  ionospheric  mapping  of  an  injection  boundary  in  the  equatorial 
plane  determined  empirically  by  Mcllwain  (1974),  as  indicated  in  Figure  6.2. 


— 
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Figure  6.2.  Data  of  Figure  6.1  along  with  the  Mcllwain  injection  boundary 
mapped  to  the  auroral  zone.  Two  field  models  shown,  dipole 
(solid  line),  and  empirical  Mead-Fairf ield  (MF73).  From  Sheehan 
and  Carovillano,  1978. 


Thus,  the  equatorward  boundary  seen  in  DMSP  images  is  related  to  plasma 

injected  near  the  inner  edge  of  the  plasma  sheet  at  the  onset  of  substorms,  at  least 

statistically. 

Cussenhoven  et  al.  (1978)  carried  out  a  similar  analysis  of  auroral 
boundaries  in  DMSP  images  but  concentrated  on  the  poleward  boundary,  which 
is  more  difficult  to  clearly  determine  than  the  equatorward  one.  The  authors 
also  considered  D  behavior  (Figures  6.3  and  6. A,  at  two  local  times)  with 
the  intention  of  estimating  ring  current  influence  on  polar  cap  size.  They 
found  that  separating  Dst  affects  from  AE  affects  was  difficult  to  accomplish 
because  the  general  relationship  between  AE  and  Dst  hampered  the  accumulation 
of  good  statistics  at  extreme  values.  Nevertheless,  the  task  of  mapping  auroral 
locations  globally,  under  various  magnetic  conditions,  remains  of  great 
importance  in  studying  the  overall  state  of  the  magnetosphere. 


MOST  EQUATQRWARD  AND  MOST  POLEWARD  BOUNDARIES  IN  CGL 
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Characteristics  of  the  Equatorward  Auroral  Boundary  Near  Midnight 
Determined  From  DMSP  Images 

R.  E.  Sheehan  and  R.  L.  Carovillano 

Department  of  Physics.  Boston  College.  Chestnut  Hill,  Massachusetts  02lf>7 

The  latitude  of  the  equatorward  auroral  boundary  near  local  midnight  has  been  determined  for  162 
Defense  Meteorological  Satellite  Program  (DMSP)  images  in  November-December  1972  When  grouped 
according  to  Kp  and  AE.  these  observations  show  approximate  linear  decreases  in  the  average  boundary 
latitude  with  increasing  values  of  these  magnetic  indices  There  appears  also  to  be  a  slight  diurnal 
variation  in  the  boundary  location.  Mapping  of  the  appropriate  Mcllwain  injection  boundaries  to  auroral 
latitudes  shows  good  agreement  with  the  average  DMSP  equatorward  auroral  boundary  latitude.  Similar 
analyses  at  2000  and  2200  COLT  (corrected  geomagnetic  local  time)  using  a  different  set  of  DMSP  images 
yield  similar  results,  with  somewhat  pooicr  agreement  under  quiet  conditions. 


Introduction 

Auroral  images  from  satellites  have  been  especially  useful  in 
studies  of  global  auroral  patterns  during  substorms  [Pike  and 
Whalen,  1974;  Snyder  el  at..  1974;  Akasofu.  1974;  Carovillano 
et  al  .  1974]  Isis  satellites  provide  valuable  spectral  informa- 
Itot.  | /.mi  et  ul .  I975u.  6],  while  U  S  Air  Force  satellites  in  the 
Defense  Meteorological  Satellite  Program  (DMSP)  series  pro¬ 
duce  htgh-resolulion  images  in  a  broad  band  extending  from 
long  visible  lo  near-infrared  wavelengths.  Some  studies  are 
based  on  a  small  number  of  images  selected  lo  complement 
other  data  Work  involving  large  groups  of  DMSP  images 
| Hot:  north  and  Meng.  1975;  Carovillano.  1975;  Berkey  and 
Ktwtide.  1976]  requires  extensive  gridding  and  analysis, 
straightforward  but  tedious  by  hand  and  expensive  by  com¬ 
puter  I  sen  so.  statistical  studies  exploring  the  gamut  of  mag- 
netosphene  conditions  become  more  inviting  and  possible  as 
large  quantities  of  images  become  available. 

The  primary  measurement  reported  in  this  paper  is  the 
latitude  of  the  equatorward  auroral  boundary  as  seen  in 
DMSP  images  Orbiting  in  circular  paths  with  102-min  peri¬ 
ods.  DMSP  satellites  travel  in  either  the  noon-midnight  or  the 
dawn-dusk  meridian  plane  A  radiometer  pointing  toward 
earth  builds  an  image  by  rapidly  scanning  narrow  strips  per¬ 
pendicular  to  the  satellite  path  projected  to  the  earth's  surface 
(the  ground  track)  The  apparent  intensity  of  auroras  seen  in 
DMSP  images  depends  on  several  factors;  absolute  radiometer 
sensitivity,  original  image  and  copy  quality  (in  our  case.  35 
mm  transparencies  from  the  World  Data  Center),  oblique 
viewing  angles  away  from  the  center  of  the  ground  track  (Van 
Rhijn  effect ).  and  uncertain  background  albedo  from  ice  and 
clouds  (see  Carovillano  (1975]  for  a  discussion  of  these  effects) 

Our  data  base  consists  of  162  DMSP  images  taken  by  a 
noon-midnight  satellite  during  the  period  November  14  to 
December  15.  1972  A  DMSP  satellite  can  produce  a  maxi¬ 
mum  of  about  430  images  in  a  month,  but  gaps  in  Ihe  data  arc 
rather  common,  and  some  images  are  of  doubtful  quality  for 
reasons  related  to  the  problems  mentioned  above.  In  addition, 
exposure  problems  and  interference  by  the  full  moon  com¬ 
monly  result  in  poor-quality  images  The  data  base  in  this 
study  was  chosen  because  it  contained  a  good  sample  of  high- 
quality  images  Our  procedure  was  to  grid  each  image  with  a 
corrected  geomagnetic  (CG)  coordinate  system  [Whalen.  1970] 
and  to  measure  the  latitude  of  the  equatorward  auroral  bound¬ 
ary  at  the  0100  COLT  (CG  local  time)  meridian  This  meri- 
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dian  is  convenient  because  of  its  proximity  to  the  ground  track 
for  images  from  a  satellite  orbiting  in  the  noon-midnight 
plane.  Because  the  satellite  is  stabilized  in  a  fixed  plane  relative 
to  the  earth-sun  line,  the  magnetic  coordinates  in  a  series  of 
images  shift  diurnally  as  the  earth  rotates  below,  and  this  effect 
is  taken  into  account  by  the  gridding  process. 

Generally,  diffuse-appearing  aurora  served  to  define  the 
measured  equatorward  boundary.  However,  DMSP  images 
cannot  accurately  portray  the  entire  range  of  possible  auroral 
intensities.  Saturation  occurs  at  relatively  low  brightness  levels 
and  often  hinders  discrimination  between  diffuse  and  discrete 
aurora,  based  on  structure  alone.  When  a  distinction  is  pos¬ 
sible,  it  is  based  solely  on  the  actual  DMSP  images  and  not  on 
complementary  information  When  comparisons  have  been 
possible,  our  results  have  been  quite  consistent  with  features 
seen  in  all-sky  camera  pictures  Whether  diffuse  or  discrete  in 
appearance,  the  boundary  is  usually  quite  sharp  and  is  mea¬ 
sured  within  i°  of  latitude 

Histograms  With  Kp,  AE.  and  UT  as  Paramltlrs 

Figures  la,  16.  and  Ir  contain  histograms  showing  the  dis¬ 
tribution  of  DMSP  images  by  equatorward  auroral  boundary 
latitude  when  grouped  parametrically  according  to  the  Kp 
index,  the  AE  index,  and  image  universal  time,  respectively.  In 
this  format,  each  image  has  four  quantities  associated  with  it: 
equatorward  boundary  latitude,  universal  time,  and  the  appro¬ 
priate  Kp  and  AE  indices  For  increasing  values  of  Kp  and  AE 
the  histograms  display  shifts  in  the  distributions  from  about 
67°  CGL  (quiet  conditions)  to  about  62°  CGL  (active  condi¬ 
tions)  (CGL  is  corrected  geomagnetic  latitude)  There  are  only 
a  small  number  of  cases  corresponding  to  large  values  of  Kp  or 
AE.  but  the  trends  are  both  unmistakable  and  expected  from 
the  definitions  of  the  indices.  The  distributions  for  Kp  and  AE 
are  sharply  peaked  at  67°  CGI.  during  very  quiet  conditions, 
but  in  all  cases  where  there  are  more  than  a  few  images  Ihe 
width  of  each  distribution  spans  about  5°  of  latitude  Large 
scatter  is  understandable  when  Kp  is  used  as  a  parameter 
because  it  is  a  3-hour-averaged  mid-latitude  magnetic  index 
and  therefore  may  not  represent  the  wide  degree  of  auroral 
variability  possible  within  a  3-hour  interval  The  AE  index, 
however,  more  nearly  provides  an  instantaneous  indication  of 
magnetic  activity  somewhere  in  the  auroral  zone.  The  scatter 
tn  the  distributions  for  the  boundary  latitude  for  50-y  AE 
ranges  means  lhat  there  is  a  general,  but  by  no  means  unique, 
connection  between  general  high-latitude  magnetic  conditions 
and  the  auroral  boundary  location 
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II  auroral  activity  occurred  randomly  at  any  lime  of  the  day, 
(he  distribution  of  images  by  boundary  latitude  Mould  be 
broad  and  flat  when  grouping  was  done  according  to  the 
universal  time  of  the  image.  Although  the  small  number  of 
cases  in  each  UT  interval  in  Figure  l<  makes  it  difficult  to  draw 
firm  conclusions,  the  data  in  our  sample  tend  to  suggest  that 
the  auroral  equatorward  boundary  more  commonly  extends  to 
low  latitudes  in  the  latter  half  of  the  UT  day.  Images  taken 
around  OtltK)  UT  may  be  biased,  however,  because  at  this  time 
the  tlltKl  t'Gl.T  meridian  is  in  a  position  prone  to  Van  Rhijn 
enhancement.  If  real,  the  small  diurnal  effect  exhibited  in 
t  igurc  U  could  be  a  reflection  of  a  diurnal  change  in  the 
average  polar  cap  si/c  or  possibly  of  a  misalignment  of  the 
corrected  geomagnetic  coordinate  system  with  the  average 
uuroiul  oval  The  si/e  of  the  polar  cap  cannot  be  determined 
accurately  from  DMSP  images  except  in  rare  cases  because  ol 
their  limited  spatial  coverage,  but  similar  trends  with  no  phase 
lags  at  several  other  local  times  would  indicate  an  overall 
expansion  and  contraction  of  the  average  nightside  auroral 
oval  If  svstcni.uic  phase  differences  among  the  diurnal  trends 
at  several  local  times  were  observed  this  behavior  would  vug- 
pc  a  111  it  a  const  aril  average  ov  at  pertains  which  sin  Its  location 
relative  to  an  eccentric  land  perhaps  nonlinear)  coordinate 
sv stem  as  the  earih  ro  ales 


ilMOHAI  Bill  SIIAKV  NlAK  MlOSU.III 

Variation  iihhi  Avntxra  Fijcatorw aro  Bm  miars 
Laiitldi  With  Kp  and  At 

Figure  2  displays  the  average  and  standard  dev  union  of  the 
equatorward  auroral  boundary  latitude  distributions  at  0100 
CGI  T  versus  Kp  Here  integer  Kp  values  encompass  neighbor¬ 
ing  fractional  cases,  and  the  Kp  -  I  data  point  also  includes 
images  with  Kp  =  0  The  decrease  in  ihc  average  latitude  of  ihe 
equatorward  boundary  varies  approximately  linearly  wiih  Kp 
when  I  £  Kp  5  When  the.,  results  are  compared  wuh  the 
results  of  Lui  cl  al  11975/))  from  Isis  2  data,  a  general  agree¬ 
ment  between  the  two  types  of  imaging  is  found  (figure  .'I 
Also,  the  equatorward  trend  over  the  Kp  range  I  •  '  is  com¬ 
parable  to  that  shown  by  Starkor  [1969)  based  on  all-sky 
camera  data  for  Q  index  values  ranging  from  2  to  6  On  the 
other  hand,  accord  is  not  as  good  with  an  analysis  ol  south 
pole  DMSP  images  during  Mav-June  1973  by  Sagala  el  al. 
[1975)  As  shown  in  Figure  4.  the  trends  are  similar,  but  our 
average  latitudes  arc  l°-2“  lower  for  Kp  =  2.  3.  and  4  These 
differences  could  arise  from  the  particular  sets  of  images  used 
in  each  study  or  physically  from  a  lack  of  conjugacy  between 
the  northern  and  the  southern  hemisphere  in  the  corrected 
geomagnetic  system.  The  change  in  the  average  equatorial 
latitude  versus  At'  (Figure  Alcan  be  approximated  roughly  by 
a  linear  decrease  of  about  1°  (  GL/  100  y 

Only  on  exceptional  occasions  Joes  the  aurora  extend  below 
60°  CGI  For  example,  during  the  great  magnetic  storm  pe¬ 
riod  in  early  August  1972  the  equatorward  boundary  ap¬ 
proached  55°  CGI  only  during  an  extremely  disturbed  inter¬ 
val  when  A/i  reached  8  A  linear  extrapolation  of  Figure  2  to 
Kp  >  5  leads  to  a  similar  value 

Pi  ASMA  INJI  nttlN  PltKNOMliNA 

Observations  of  magnetosphcric  plasma  at  sy  nchronous  or¬ 
bit  ('6.6  R, )  have  been  reported  by  DeFonsi  and  MclUam 
[1971  ]  Those  authors  exhibit  their  data  in  the  form  of  energy- 
time  spectrograms  and  give  fluxes  of  electrons  and  protons  in  a 
range  of  energies  as  a  function  of  time.  Characteristic  dis¬ 
persion  structures  are  seen  in  these  spectrograms  that  indicate 
that  the  satellite  encounters  protons  and  electrons  of  different 
energies  at  various  points  along  its  orbit  As  the  satellite  moves 


!  tp  2  Avcf.ipi  ior»cv»«ril  ytnm-iynctu.  latitude  of  On  >  tj'i.tinf  ss.ud 
.idMir.il  Hound  trv  .il  0100  (  ( il  I  .is  .1  fun*  lion  nl  A /»,  usniy  On;  viinr 
data  h;i\c  as  *.is  used  in  I  ijfurr  I  I  rror  O.itn  iknurk  oik  standard 
dc*  i.itmn 
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I  ,(!  <  Daij  of  I  i|turc  7  (large  dots)  superimposed  on  the  results 
al  /  j ,i  i-t  jl  [  |d-\x/>|  using  Isis  7  images  Inch  small  dot  represents  one 
Isti  2  image  A  indicates  invariant  latitude,  which  is  very  similar  to 
corrected  geomagnetic  latitude 


Irortt  evening  local  time  toward  midnight,  high-energy  protons 
are  encountered  Itrst.  followed  by  successively  lower  energy 
protons  until  there  occurs  a  simultaneous  encounter  with  near- 
/cro  energy  protons  and  electrons  (the  inner  edge  of  the 
plasma  sheet)  Occasionally,  protons  and  electrons  of  all 
energies  show  enhanced  fluxes  simultaneously,  as  if  a  cloud  of 
heated  plasma  suddenly  appeared  or  was  created  at  the  satel¬ 
lite's  location.  This  is  called  an  injection  event. 

By  calculating  the  particle  drift  motion  caused  by  magnetic 
field  curvature  and  gradient  and  adding  the  convection  of  a 
postulated  electric  field.  DeForest  and  Mcllwatn  were  able  to 
trace  particle  trajectories  backward  in  time  from  their  loca¬ 
tions  at  the  satellite  encounter  at  6.6  Rt.  If  protons  and  elec¬ 
tions  of  all  energies  are  traced  back  in  the  equatorial  plane  to  a 
common  lime  of  origin,  (heir  locus  falls  on  a  smooth  spiral 
curving  toward  the  earth  front  dusk  to  midnight.  The  locus 
(the  so-called  injection  boundary )  may  he  interpreted  to  be  the 
earthward  boundary  of  a  common  source  region  for  particles 
of  all  energies  at  the  time  of  injection.  Mcllwain  (1974)  derived 
a  formula  for  Ihe  statistical  injection  boundary  given  by  R,  = 
il?7  10  Kp)  (*-  7  1)  K,  is  the  injection  boundary  position 

in  units  of  earth  r.tdn  at  ihe  equatorial  plane  as  j  function  of 
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fig  5  Same  as  I  igurc  2  except  fur  At.  based  on  data  of  f  igurc  16 


Kp  and  magnetic  local  lime  in  hours,  e  (1800  S  <  2400 

COLT).  Mauk  und  Mcllwatn  1 1974)  analyzed  their  ATS  5  data 
in  this  way  using  the  value  of  Kp  at  the  lime  of  encounter  with 
cold  plasma  and  regard  /?„  to  be  useful  in  describing  the 
injection  boundary's  location  when  a  substorm  occurs 
Coordinated  ground-based  photometer  and  ATS  5  studies  by 
Rather  rt  al  [1976)  give  further  significance  to  the  injection 
boundary  concept 

DMSP  Auroral  Boundary  and  Tilt  Melt  wain  Injection 
Boundary  Near  Local  Midnight 

In  Figure  6  a  comparison  is  made  between  the  DMSP 
boundary  location  and  the  Mcllwain  injection  boundary  /?„  at 
0000  COLT  as  a  function  of  Kp.  The  formula  for  /?„  is  valid 
only  up  to  local  midnight,  and  its  extension  past  local  mid¬ 
night  is  questionable  Although  our  DMSP  data  apply  to  0100 
instead  of  0000  CO  I  T.  ibis  should  be  of  little  practical  signifi¬ 
cance,  since  the  auroral  boundary  is  generally  parallel  to  cir¬ 
cles  of  constant  latitude  in  the  vicinity  of  local  midnight  A 
simple  dipole  model  (solid  line)  and  an  empirical  field  model 
(dashed  line)  are  used  to  map  the  injection  boundary  from  the 
equatorial  plane  to  the  auroral  zone.  The  realistic  magnetic- 
fields  used  in  Figure  6  are  the  Ml  Q  (quiet.  Kp  <  2)  and  MFD 
(disturbed,  Kp  ?  2)  models  described  by  h'air/teld  and  Mead 
1 1975)  In  a  review  of  magnetic  field  models.  Walker  (1976) 
considered  the  Mead-1  airfield  models  (M FQ  and  MI  D)  to  be 
among  the  best,  principally  because  they  are  based  on  satellite 
data 

In  the  equatorial  plane  the  most  accurate  region  of  the 
Mead-Fairlield  model  is  -5  Rr,  2  X  ?  - 12  Rf,  where  .V  is  the 
distance  from  the  earth's  center  along  the  earth-sun  line,  posi- 


I  ig  b  Data  nf  1  igurc  2  along  wiih 'lie  Mcllwain  in|tct"<h  bon  d 
ary  mapped  In  the  auroral  /nnc  using  ilipnte  magnetic  field  nu.dtl 
luil  It)  line)  and  'he  crnpirn  al  Mead  I  a  it  eld  7  '  I M  t  J'i  "  ir"-' 
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\  "  I  holograms  showing  the  number  ol  PMSP  images  with 

.  \  .ird  ii.i  oral  boumf.ii  ics  m  the  indicated  I  0  interval  of  C'( i I  ol 

;  U  s  low  i!  tunes  r..ngmg  limn  |WMi  to  OtMMi  (  til  1  Images  are 
vil  .lAi'r^liiu'  to  p. evading  A”  I  fiangics  mark  the  latitude  ol  the 
I  ciii'ieui  si.  1  *•; ^  d  oar.-r.tl  oval  at  the  unhealed  local  lime  lor  various 
(»  v  v  alt.es  |  II  ii.s.'t  u  (J  i ■  :«!«•  x  *  all-  s  1 .  1  and  5  (and  “  for  Kp 

» •  a i isp« unt  to  Hu  tele-  ivrc  bom  h'glu-st  to  lowest  latitude  tin 


]  lively  directed  toward  the  sun  At  0000  (  ( j  |  T.  for 0  r  Kp  \  5. 
she  range  oj  the  Me  1 1  wain  boundary  in  the  mghtside  cqua 
tori.J  plane  is  4  J  R,  '•  K*.  *  h  7  Ht.  which  falls  m.nnlv 
]  in i* 1  the  field  i m Kiel  s  most  .uait.ilc  range  It  is  evident  frt>m 
^  I  ie  ^  b  th.it  there  is  good  agreement  between  the  average 
■  «  OMNI*  .ioror.it  boundary  location  .md  the  mapping  ol  H„  by 
.5  the  vmmncal  Mcad-I  airfield  magnetic  held  model  1  he  undiv 
3  lotted  dipole  m.jppjne  ol  K*  also  compares  quite  well  at  all  Kp 
X  -dues  It  should  be  noted  that  the  formula  lor  the  Mcllwam 
nicvOt-n  boundary  defines  a  curve  at  the  instant  o/  infection 
\::cr  micciion  the  bound.irv  position  shills  m  response  to 
chaiiK  1  chls  and  particle  drill  motions  It  a  OMNI*  image  is 

•  iKca  some  time  atur  an  t  meet  ion  event  occurs,  the  auroral 
houm!ar\  position  should  reflect  this  change  Thus  to  the 
extern  that  the  auroral  boumJjrv  is  related  lo  the  infection 
bound.irv  the  scatter  depicted  in  I  igurc  b  could  result  from 
tw-»  -onrees  ( I  >  statistical  variations  m  the  initial  injection 
boundary  location  and  i  2 )  modifications  of  the  boundary  I  oca 

•  .  after  the  mjCvtu  *i  event  Naturally,  one  would  expect  to 

see  larger  modifications  in  the  hmindary  as  the  time  tag  be- 
i  a  ecu  ,o:  nuectioii  c'CMt  and  a  subsequent  l)MSI*  image  in- 
‘ r  •'  ■'«  '  lor  i  large  lumber  of  cases  under  given  conditions 


the  average  tunc  lag  would  depend  on  the  relative  Wcqocnc)  of 
DMSP  images  and  injection  events  Our  observations  show 
that  in  spue  of  these  factors  the  average  equator*  rd  auroral 
houndary  relates  well  io  the  empirical  boundary  found  by 
Mcllwam 

A  study  by  Ku/ntJc  and  B tnninyhani  vm  the  low- 

latitude  houmlarv  of  *  100-cV  electrons  compared  ihe  shape 
and  location  of  these  boundaries  protected  onto  (lie  equatorial 
plane  with  b  mudaries  of  calculated  Alfvcn  layers  An  Alfvcn 
laser  dchnea  cs  the  closest  approach  to  earth  ol  plasma  mov¬ 
ing  earthward  irom  the  tail  and  depends  on  the  mpressed  and 
corotation  electric  fields  as  well  as  particle  energy  and  mag¬ 
netic  moment  I  he  low-energy  electron  boundaries  tracked 
well  with  'he  \llven  layer  for  the  appropriate  electric  fields 
that  were  estimated  from  h  -urlv  averaged  li.  ol  the  inter¬ 
planetary  magnetic  field  1  his  lact  is  m  accord  wth  our  results 
using  the  emfnrical  boundary  that  depends  on  Kp  1  he  above 
authors  also  commented  that  their  findings  compared  favor¬ 
ably  with  those  of  l  ui  <  t  al  |  h\  cited  earlier,  except  that 

the  particle  boundaries  were  distributed  somewhat  lower  in 
latitude  than  the  optical  boundaries  seen  by  l*is 

I)  MS  I*  Al  KURM  Bo  l  Ml  SHY  AM)  TIK  M(  1 1  VS  AIN  Is  Jl  »  T  ION 
Bo  I  NDAKY  I’Kl  NUDNK.HT 

Because  the  known  injection  boundary  extends  from  dusk  !o 
midnight,  a  comparison  with  the  auroral  equator  ward  bound¬ 
ary  can  also  be  made  in  that  region  to  sec  whether  the  agree¬ 
ment  found  above  extends  into  the  evening  sector  Histograms 
of  the  DMSP  equatorward  auroral  boundary  latitude,  as  func¬ 
tions  of  both  Kp  and  local  times  ranging  from  IKOO  to  2400 
C  Cil  T.  are  given  in  I  igurc  7  A  DMSP  satellite  in  a  dawn- 
dusk  orbital  plane  provided  images  during  the  same  time 
period  as  was  used  above  (November  14  to  December  15, 
1972).  Its  orbit  allows  nearly  daylong  coverage  ol  the  evening 
local  time  sector  from  2000  to  7 TOO  t  GL  T  Triangles  at  de¬ 
creasing  latitudes  in  f  igure  1  locate  the  equatorial  boundary  ol 
the  statistical  I  cidstein  oval  as  given  hv  H  halvn  (1970)  for  Q 
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I .  *  .ind  5  {)  -  7  is  jlso  included  for  hp  H  i  The  Q  index  is  a 
I'  nun  .iuror.il  /one  magnetic  index  employed  by  Russian 
s*i'rkers  and  is  rs'ughis  cquualeni  lo  Kp,  especially  during 
quid  eoiidil.ons  [father  artj  MenJe  1971)  Although  the  Q 
index  does  not  correspond  cxjclls  to  hp.  it  is  obvious  that  the 
I7MSI*  bound.irs  distributions  h.ise  a  local  tunc  dependence 
similar  to  that  of  the  auroral  oral 

I  igurc  K  goes  the  latitude  of  the  auroral  and  Mcllwam 
injection  boundaries  at  7000  and  2700  C  trl  T  in  the  same 
format  as  ssas  used  in  I  igurc  6  In  Kigure  8  the  small  kink  seen 
in  the  Mcad  Knrheld  7 3  (Ml  7Jy  field  mapping  results  from 
the  transition  from  the  Ml  'MQ  to  the  Mb  73D  model  at  hp  = 
7  I  or  hp  •  3  the  DMSP  auroral  boundary  agrees  ssell  with 
the  field-mapped  Mclmam  boundary  When  data  from  the 
two  DMSI*  satellites  are  combined,  xsc  conclude  from  f  igures 
(•  and  s  that  for  model  ate  and  high  hp  values  the  average 
equ  Itorward  auroral  boundarv  maps  closely  onto  the  statisti¬ 
cal  Mcllxxuin  injection  bound.irs 

"s  also  observe  in  figure  S  that  agreement  between  the 
statistical  micvtion  and  auroral  boundaries  worsens  toward 
low  hp  > alues  i  hp  d  2 1  Aside  Irom  purely  statistical  errors, 
there  aic  1“,.  possible  explanations  lor  this  trend  ( I  )  Some  of 
the  low  hp  images  represented  in  figure  8  (about  20‘7  l  have  a 
reduced  sensitivity .  as  ev  idenced  hv  la  in  (ness  of  city  lights  and 
background  illuminaii  >n  Such  images  would  miss  faint  aur¬ 
oras  from  precipitation  occu'rmg  several  hours  alter  the  last 
mieclion  and  therefore  would  create  a  bias  under  low  Kp 
vondiiions  (?)  During  intervals  of  low  Kp,  injections  arc  rare, 
theiclorc  onlv  ci  lew  DMSP  images  ma>  show  auroral  condi¬ 
tions  ditccilv  connected  with  an  injection  To  compare  DMSP 
results  with  the  Mcllwam  boundarv  under  quiet  conditions, 
"c  should  select  i>n In  those  im  igcs  taken  soon  after  the  infre¬ 
quent  iniCCtnMis  that  do  occur  In  tact  this  procedure  should 
be  lived  under  ad  conditions,  bat  because  high  Kp  is  gencruMv 
associated  with  suhstorms  and  micctions.  most  of  the  images 
.it  such  tunes  will  auto'ualicallv  be  within  several  hours  o!  the 
mi'si  icweni  iniedior  event 

(  •  *\1 !  VKISDN  v\  1 1 1#  (  >KiM  ND-Bvsl  1)  Pllol  OStl  I  I  K  STIDIIS 

Observations  relating  the  mieclion  boundarv  to  auroral  pre- 
vipit.ihon  have  been  tarried  out  proiouslv  lor  example, 
ground  based  photometric  measurements,  made  Irom  the  foot 
ol  the  A  1 S  s  tield  line  bv  t.uthvr  el  j l  |!^76|.  have  been 
compared  with  particle  fluxes  observed  bv  A7S  5  at  syrehro- 
noiis  orbit  I  heir  auroral  observations  span  latitudes  above 
nut  below  the  \  IS  *  held  line  projection  and  provide  informa- 
i  on  ibowt  m agnctosphcric  regions  remote  from  the  ATS  satel¬ 
lite  !  hev  observe  that  infections  occurring  at  ATS  5  corrc 
spond  I**  auroral  enhancements  at  the  foot  of  the  ATS  5  field 
tme  and  amoral  enhancements  at  other  latitudes  are  injections 
wh.vti  x  in  be  interred  from  the  dispersion  in  particle  arrivals 
Steady  we  ill  auroras  correspond  to  electron  precipitation 
bom  the  inter  edge  of  the  plasm  a  sheet  On  the  }  quiet  days  in 
iheti  suidv  period,  f  ather  el  a!  observed  a  premidnight  au- 
forai  boundary  ranging  from  to  69“  at  2000  ??00  i'G I  1 
)  igu’e  "  shows  that  DMSP  images  exhibited  mans  cases  of 
bound  0 irs  n  this  are  -  when  Kp  0  or  I .  although  there  arc 
>1-  is.  >  we  poleward  I  hose  oioral  boundaries  which 
•v  v  ;■  .  :  ;  ides  greater  ;*■  ■  ■  "n  l  (il  are  subst antiaUv 

p  dew  >t .t  ot  the  MclUa  m  in-cvt.on  tvi'iinil.irx  regardless  ot 

'lu  '  ■  o  >  "  uld  u.cd  Such  jiiH't.iv  >v  i  hi  Id  tall  outside  of.  or 

.  -  Hr’  it.  'he  1  ii  'll.'  n  vh  ot  the  ground  Wised 

;  '  n-  -t.'-t  >  '  <i,l 

I  '  •  i..''.."i.  |IMM’ I';-  ,•  di'U'vi  j,.ri,r,i.  v„rtr 


spending  to  steady  plasma  sheet  precipitation  could  result 
from  the  exposure  problems  discussed  earlier.  However,  an 
explanation  is  still  needed  for  the  presence  of  intense  arcs  (by 
DMSP  standards)  observed  poleward  of  the  projected  Mcll- 
wain  boundary  If  such  instances  are  related  to  injections,  then 
they  may  be  examples  of  the  so-called  ‘contracted  oval  sub- 
storms'  ( hamttlc  et  a /  .  1 975.  Akasofu.  1 975  J  suggested  to  origi¬ 
nate  well  back  in  the  plasma  sheet  This  proposed  type  of 
subslorm  occurs  well  poleward  of  most  auroral  /one  observa¬ 
tories  and  is  reported  to  coincide  with  a  northward  pointing 
interplanetary  magnetic  bold  Poleward  injections  giving  rise 
to  contracted  oval  suhstorms  do  not  contradict  the  observa¬ 
tions  of  father  ft  til  |I9T(>].  because  they  are  presumed  to 
occur  within  the  preexisting  plasma  sheet,  although  further 
down  the  tail  than  usual 

Si  MMAHX 

Measurements  of  the  eqtiatorward  auroral  boundary  lati¬ 
tude  near  local  midnight  from  162  DMSP  images  show  that 
the  average  latitude  varies  linearly  with  hp  or  Ah,  decreasing 
from  67°  CGI  to  about  61°  CGI  for  hp  <  5  or  AE  <,  600 
The  rather  large  scatter  around  the  average  positions  indicates 
that  the  auroral  boundary  is  only  roughly  related  to  ground 
magnetic  conditions  However,  a  mapping  of  the  appropriate 
Mcllw.un  injection  boundary  from  the  equatorial  plane  to 
auroral  latitudes  produces  a  curve  that  coincides  fairly  well 
with  average  aiooral  boundaries  at  local  midnight  for  hp 
ranging  from  I  to  5  At  2000  and  2200  CGI.T  the  agreement  is 
good  for  hp  -  3-5.  but  the  auroral  boundary  averages  about 
2°  poleward  of  the  Mcllwatn  boundary  during  quieter  condi¬ 
tions  ( hp  0-2) 
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DMSP  POLEWARD  AURORAL  BOUNDARIES  AND 
MOVEMENTS  WITH  MAGNETIC  ACTIVITY 

M.  S.  Gussenhoven,  M.  D'Anton,  R.  L.  Carovillano  (all  at: 
Department  of  Physics,  Boston  College,  Chestnut  Hill,  MA 
02167) 

The  DMSP  images  for  October,  1974,  are  used  to  study  the 
variation  of  the  poleward  edge  of  auroral  displays  and  their 
latitudinal  extent  with  a  variety  of  ground  magnetic  indices. 
The  variation  with  Dst  relates  to  the  influence  of  the  ring 
current  on  high  latitude  accessibility  of  auroral  particles. 
Consecutive  DMSP  images  fo:  this  period  (which  includes  a 
magnetic  storm  and  several  disturbed  days)  are  gridded  in 
corrected  geomagnetic  coordinates  and  local  time.  Two  DMSP 
satellites  (8531  and  9532)  were  recording  at  this  time.  The 
most  poleward  and  the  most  equatorward  edge  of  the  aurora  shown 
for  each  hourly  interval  of  local  time  from  1800  to  0600  hours 
through  the  midnight  sector  are  recorded.  The  well-established 
equatorward  motion  of  the  equatorward  edge  of  the  aurora  near 
midnight,  as  magnetic  activity  increases,  is  shown  also  to  be 
the  case  at  other  local  times  studied.  The  poleward  motion  is 
somewhat  more  complex  and  less  regular;  however,  a  significant 
movement  of  the  auroral  boundary  towards  higher  latitude  with 
increasing  Dst  is  found  at  all  local  times. 
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7.  MAGNETOSPHERIC  DETERMINANTS  OF  THE  GLOBAL 
DISTRIBUTION  OF  PARTICLE  PRECIPITATION 

A.  Introduction 

The  goal  of  this  research  has  been  to  identify  and  quantify  the 
magnetospher ic  factors  that  determine  where  the  particle  precipitation 
characteristics  of  the  auroral  zone  occur  at  any  given  time.  New  observa¬ 
tional  characteristics  of  the  auroral  zone  were  presented  in  Chapter  6.  In 
addition,  the  statistical  or  average  auroral  zone  location  is  well  known, 
as  is  the  fact  that  displacements  from  the  average  are  common.  The  magneto- 
spheric  processes  that  govern  the  displacements  are  the  subjects  of  concern 
here . 

In  the  simplest  descriptions,  the  geometrical  (convertible  to  geographical) 
quantities  of  interest  are:  I)  the  opening  angle  of  the  poleward  limit  of  the 
precipitation  band,  2)  the  opening  angle  of  the  equatorward  limit  of  the 
precipitation  band,  and  3)  the  displacement  of  the  center  of  the  band  (in  the 
sense  of  the  center  of  a  circle)  from  the  dipole  axis.  These  are  essentially 
latitudinal  variations.  Variations  in  longitude  were  not  included  in  this  first 
(three  part)  phase  of  the  investigation  except  as  they  arise  from  the  displace¬ 
ments  of  the  centers  of  the  auroral  bands  from  the  geographical  poles. 

The  first  two  quantities  listed  are  believed  to  be  controlled  by  the 
orientation  and  strength  of  the  interplanetary  magnetic  field  and  by  the 
magnetospher ic  ring  current.  The  ring  current  contribution  is  particularly 
important  in  accounting  for  the  largest  displacements  of  the  auroral  precipitation 
from  its  average  location. 

A  model  describing  how  the  ring  current  controls  the  opening  angle  of  the 


poleward  edge  of  the  precipitation  band  has  been  proposed  under  this  contract. 
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It  requires  for  its  use  knowledge  of  the  magnetic  moment  of  the  ring  current 
(M^) .  in  this  model,  the  location  of  the  equatorward  edge  of  the  auroral  zone 
is  identified  with  the  inner  edge  of  the  ring  current  (fO. 

Existing  ring  current  models  are  not  of  much  use  for  our  purposes,  because 
they  require  that  the  entire  ring  current  particle  distribution  be  specified 
as  an  input  function.  That  function  must  be  found  empirically,  by  jin  situ 
spacecraft  measurements,  or  by  other  means  that  are  extraneous  to  the  models 
themselves.  It  was  necessary,  therefore,  as  a  first  stage  of  this  investigation 
to  develop  a  ring  current  model  in  which  the  ring  current  evolved  self-consistently 
to  a  given  (presumably  storm-time)  state  from  a  prior  (presumably  prestorm)  state. 

h.  Quasi-Self-Consistent  Ring  Current  Model 

In  this  section,  wo  summarize  the  quasi-self-consistent,  high  beta,  injec¬ 
tion  model  for  the  storm-time  ring  current  that  has  been  successfully  developed 
by  SjLsc9c  (1979a).  It  has  the  advantage  over  previous  models  of  completely 
determining  all  ring  current  parameters  in  terms  of  two  observables:  1)  the 
geomagnetic  Ost  index,  and  2)  the  beta  of  the  ring  current  plasma  at  a  single 
(and  arbitrary)  distance.  The  ring  current  beta  can  be  determined  at  any  time 
prior  to  the  ring  current  decay  phase-even  prior  to  the  development  of  the 
storm.  If  this  parameter  (beta)  is  the  same  for  all  storms,  then  only  Pst  is 
needed  to  specify  all  ring  current  parameters.  This  possibility  can  he  investi¬ 
gated  empirically. 

Some  results  obtained  from  t  lie  model  are  shown  in  Figures  7.1  to  7.3. 

Figure  7.1  shows  the  calculated  relationship  between  list  and  the  geocentric 
distance  to  the  inner  edge  of  the  ring  current  (Rr)  in  units  of  earth  radii 
(R  ) .  The  family  parameter,  R,  is  the  second  of  the  two  model  parameters.  it 


--  .•  “■ 
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The  dependence  of  the  ring  current  disturbance  field 
at  earth  (H)  on  the  distance  to  the  inner  edge  of  the 
ring  current  (R)  for  different  "entropy"  plasmas 
characterized  by  k,  the  distance  at  which  the  beta  of 
the  plasma  is  unity. 


Figure  7.1 


The  magnetic  moment  of  the  ring 
current  (MRg)  in  units  of  the  geo¬ 
magnetic  dipole  moment  (Mg)  as  a 
function  of  the  disturbance  field 
(H)  for  different  values  of  k. 


A  comparison  between  the  model- 
predicted  (R  =  4,  k  =  8)  disturbance 
field  profile  and  storm-time  distur¬ 
bance  field  data  from  the  S3-A 
satellite  for  the  storm  of  17 
December  1971  (Cahill,  1973). 


Figure  7.2 


Figure  7.3 
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Is  yj/i-n  by  -  (R  ,  which  I  ft  a  constant  of  the  model.  Figure  7.1  indicates 

that  when  K ^  moves  into  tlie  range  of  3  to  4  (R^) ,  storm-time  Dst  values  result, 

in  accord  with  available  observations.  Figure  7.2  shows  the  dependence  of  the 

ring-current  magnetic  moment  (normalized  to  units  of  the  earth's  dipole  moment) 

on  Dst  and  the  parameter  k.  Figure  7.3  gives  a  comparison  between  the  ring 

current  magnetic  field  calculated  from  the  model  and  the  only  published  example 

for  which  the  values  of  k  and  R^_  can  be  empirically  and  simultaneously  determined 

(Hoffman,  1973).  The  agreement  is  quite  good.  The  deviations  seen  at  r  <  3.5 

(R  )  can  be  attributed  to  time  variations,  which  were  not  taken  into  account 
e 

in  the  calculated  curve. 

As  these  figures  show,  given  observed  values  of  Dst  and  k  (or  in  the  case 
of  k  perhaps  a  standard  value  will  be  adequate) ,  values  of  and  Rr  can  be 
determined.  These  are  the  quantities  which  then  govern  the  values  of  the  first 
of  the  listed  geometrical  features  of  the  particle  precipitation. 

C.  Ring  Current  Effect  on  Polar  Cap  Size 

The  Dst  versus  relationship  displayed  in  Figure  7.2  has  been  used 
to  assess  the  ring  current  contribution  to  the  equatorward  expansion  of  the 
auroral  oval  during  geomagnetic  storms  (Siscoe,  1979b).  The  principle  upon 
which  the  contribution  is  determined  is  simply  stated.  The  ring  current  magnetic 
moment  is  a  vector  which  points  in  the  same  direction  as  the  geomagnetic 
dipole  moment  M^,.  The  solar  wind  interacts  with  the  coupled  earth-ring  current 
system,  which  as  a  combined  dipole  moment  >  Mg.  Now  we  simply  scale 

up  the  size  of  the  magnetosphere  from  its  size  when  there  is  no  storm  and  when 

,  using  the  well-known  dipole  scaling  relations.  A  larger  magnetosphere 
necessarily  has  a  larger  amount  of  magnetic  flux  forming  its  tail.  The  magnetic 
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flux  in  the  tail  maps  into  the  polar  caps  of  the  earth,  which  demarcate  the 
polar  boundaries  of  the  auroral  ovals.  Thus,  the  poleward  border  of  the  ovals 
move  equatorward  as  increases.  The  formula  expressing  this  relationship 
(Siscoe,  1979b)  is: 

.  ,  /,  .  V\l/3/ Pst\  1/12 

p  p  V  v  \p*t/ 


where  is  the  magnetic  latitude  of  the  poleward  edge  of  the  oval,  p  ^  is 

the  solar  wind  stagnation  pressure,  and  X*  =  70°  and  p*  =  1.2  x  10~8  dvne 

p  st 

cm-2  are  fiducial  (actually  about  average)  values.  Since  we  are  interested 


here  in  the  contribution  from  M^,  we  set  P  t  =  P*{-»  i-e-»  we  assume  an  average 

solar  wind  stagnation  pressure.  Then  we  can  plot  the  equatorward  displacement 

AX  =  X*  -  X  as  a  function  of  M^M^.  This  is  shown  here  as  Figure  7.4.  We 

see  from  this  figure  that  a  displacement  of  5°  requires  a  magnetic  moment 

ratio  near  unity,  and  a  displacement  of  10°,  a  ratio  of  approximately  two. 

The  curve  in  Figure  7.4  and  the  cuves  in  Figure  7.2  can  be  combined  to 

2 

provide  a  plot  of  the  displacement,  AX,  as  a  function  of  H(=  -j  Dst)  for 

different  assumed  values  of  the  parameter  k.  This  is  shown  here  as  Figure  7.5. 

We  see,  as  expected,  AX  increases  with  increasing  | Dst |  and  with  increasing  k. 

As  an  example,  consider  the  curve  for  which  k  *  8  and  for  the  purpose  of 
3 

illustration  let  Dst  1  y  H.  The  actual  proportionality  factor  between  Dst  and 
H  depends  on  the  strength  of  the  currents  induced  in  the  earth  by  the  storm 
and  can  vary  from  a  minimum  of  unity  to  a  maximum  of  3/2,  the  value  assumed. 


Thus,  in  the  case  of  a  large  storm  that  achieved  a  disturbance  level  given 
by  Dst  ■  -500  nT,  the  corresponding  value  of  H  would  be  -333  nT.  Using  the 
k  ■  8  curve,  we  then  find  AX  =  5.75°,  or  that  the  poleward  edge  of  the  auroral 
oval  would  be  located  at  a  geomagnetic  latitude  near  64°. 
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The  dependence  of  the  latitudinal  shift  of  the 
poleward  border  of  the  auroral  oval  on  the 
magnetic  moment  of  the  ring  current, 

(measured  here  in  units  of  the  geomagnetic 
dipole  moment,  M^) . 


The  two  cases  shown  correspond  to  initial  (FL  =  0) 
locations  of  the  border  at  \*  =  60°  and  70° 
geomagnetic  latitudes.  P 


Figure  7.4 


The  dependence  of  the  latitudinal  shift  of  the  poleward  border 
of  the  auroral  oval  on  the  ring  current  magnetic  field  (Dst~1.5  H)  , 
as  determined  from  the  1*  =  70°  curve  of  Fig.  7.4  and  the  previous 
derived  relation  hetweenPM^  and  H  (Siscoe,  1979a)  in  Fig.  7.2. 


F igure  7 . 5 
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Against  the  prediction  given  here,  we  can  compare  the  observations 
relating  to  the  storm  of  11  February  1958  published  by  Akasofu  and  Chapman 
(1962).  This  particular  storm  was  characterized  by  a  large  negative  value 
of  Dst,  which  reached  -500  nT,  and  by  a  large  equatorward  motion  of  the 
auroral  zone.  All  sky  camera  observations  stations  located  in  the  mainland 
United  States,  Canada  and  Alaska  revealed  that  during  one  interval  of  about 
20  minutes  duration,  the  northern  edge  of  the  auroral  oval  was  located  at 
approximately  55°  geomagnetic  latitude. 

Clearly  it  would  be  desirable  to  perform  similar  comparisons  for  other 
storms.  However,  the  important  conclusion  to  be  drawn  from  this  one  storm 
for  which  the  comparison  can  be  readily  made  is  that  the  effect  of  the  ring 
current  magnetic  moment  can  account  for  about  30  percent  of  the  maximum 
achieved  equatorward  displacement  of  the  poleward  edge  of  the  oval. 

We  know  already  that  the  size  of  the  oval  is  modulated  by  the  direction 
and  magnitude  of  the  north-south  component  of  the  interplanetary  magnetic 
field  (IMF).  The  IMF  control  is  independent  of  the  ring  current  control 
discussed  here.  The  study  by  Ho lz worth  and  Meng  (1975)  shows  that  the  oval 
expands  by  about  0.9°  for  every  1  nT  increase  in  the  southward  component  of 
the  IMF.  Thus  a  10°  displacement  requires  a  southward  component  of  only  11  nT, 
which  is  not  a  very  large  value  to  be  associated  with  a  major  geomagnetic  storm. 
It  is,  therefore,  not  unreasonable  to  conclude  that  in  the  case  of  the  11 
February  1958  storm  the  combined  effects  of  the  ring  current  and  the  IMF 
could  have  brought  the  poleward  edge  of  the  auroral  oval  southward  of  the 


Canadian  border. 
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D-  Polar  Cap  Off-Set 

The  radius  of  the  auroral  circle  varies  typically  between  15  and  23  degrees 
in  response  primarily  to  changes  in  the  orientation  of  the  interplanetary 
magnetic  field  (Holzworth  and  Meng,  1975).  The  auroral  circle  is  displaced 
on  the  average  about  4  degrees  in  the  anti-solar  direction  from  the  geomagnetic 
pole,  but  at  any  time  the  offset  varies  typically  within  3  degrees  from  the 
average  value  (Meng  et  al. ,  1977). 

The  geometry  is  shown  in  Figure  7.6.  The  equations  on  the  figure  give 

the  variation  in  both  the  radius  6  and  the  offset  fi  as  a  function  of  solar 

P 

wind  dynamic  pressure  p  ,  as  predicted  by  the  standard  magnetospheric  scaling 
relations.  The  weak  dependence  of  0^  and  on  p^t  implies  that  the  commonly 
observed  variations  in  0^  and  6  cannot  be  caused  by  corresponding  variations 
in  p 

st 

in  the  course  of  this  work,  we  have  demonstrated  that  the  interplanetary 
magnetic  field  does  not  exert  its  main  control  over  the  polar  cap  size  through 
the  modulation  of  the  magnetospheric  ring  current.  We  determined  that  the  ring 
current  required  to  produce  a  change  in  polar  cap  size  comparable  to  that  which 
is  commonly  obsctved  would  also  produce  a  large  geomagnetic  storm  (Dst ~-200y) , 
which  is  not  commonly  observed.  (In  fact,  the  calculation  was  performed 
to  show  that  the  ring  current  effect  can  be  important  during  magnetic  storms.) 
This  result,  not  unexpectedly,  supports  the  general  belief  that  the  variation 
in  polar  cap  size  results  from  a  net  transfer  of  magnetic  flux  into  the  tail 
from  the  dayside,  with  an  associated  change  in  the  shape  of  the  magnetosphere. 

We  note  here  that  the  presently  existing  magnetospheric  scaling  relations, 
which  can  be  used  to  investigate  the  effect  of  the  ring  current  on  the  polar 
cap  size,  are  not  able  to  treat  phenomena  that  entail  a  change  in  magnetospheric 
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shape.  This  result  motivated  the  investigation  into  two-parameter  scaling 
relations,  which  can  handle  such  problems,  as  is  described  below. 

The  anti-sunward  displacement  of  the  center  of  the  polar  cap  with  respect 
to  the  geomagnetic  pole  evidently  has  its  ultimate  cause  in  the  sunward- 
tailward  asymmetry  in  the  magnetosphere  itself,  which  is  imposed  by  the  sense 
of  the  solar  wind  flow. 

In  thinking  of  magnetic  field  lines  in  mechanical  terms,  field  lines  are 
given  the  attribute  of  tension.  The  force  bv  which  the  solar  wind  drags  out 
the  magnetosphere  to  form  the  geomagnetic  tail  is  then  envisioned  to  be 
communicated  to  the  earth  by  the  tension  in  the  field  lines  that  connect  the 
earth  and  the  tail.  These  are  the  polar  cap  field  lines.  In  light  of  this 
conception,  which  has  its  basis  in  the  theory  of  plastic  media,  it  is  natural 
to  suggest  that  the  ta ilward  displacement  in  question  is  a  direct  consequence 
of  the  tension  operating  through  the  polar  cap  field  lines. 

The  mathematical  statement  of  the  hypothesis  that  the  tailward  offset 
of  the  polar  cap  results  from  the  force  acting  between  the  earth  and  the  tail 
was  treated  analytically  for  the  case  of  a  two-dimensional  interaction.  An 
offset  of  the  order  of  0.4  degrees  was  obtained  this  way.  The  calculated  two- 
dimensional  offset  is  too  small  to  account  for  the  observations,  but  the 
differences  between  the  two  and  three-dimensional  cases  were  sufficient  to 
offer  hope  that  the  three  dimensional  offset  could  be  substantially  larger. 

The  three-dimensional  case  was  then  Integrated  numerically.  The  result 
did  not  differ  materially  from  the  previous  two-dimensional  case.  If  the  tall 
force  only  is  considered,  the  offset  is  an  order  of  magnitude  too  small.  Inter¬ 
estingly,  if  the  dayside  force  (the  Chapman-Ferraro  force)  is  used,  the 
resulting  offset  is  slightly  over  2  degrees,  or  about  50  percent  of  the  observed 
displacement.  Since  the.  dayside  force  is  dominated  by  the  dynamic  pressure  of 
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the  soLar  wind,  there  is  here  a  justification  for  ascertaining  the  degree  of 
correlation  between  the  polar  cap  offset  and  the  solar  wind  dynamic  pressure. 

Even  should  the  total  solar  wind  force  on  the  boundary  account  for 
roughly  half  of  the  polar  cap  offset,  we  would  conclude  the  magnetospher ic 
distortion  not  directly  related  to  communicating  the  solar  wind  force  to  the 
earth  must  contribute  substantially,  if  not  predominantly,  to  the  offset.  This 
result,  like  the  one  relating  to  the  polar  cap  size,  indicates  that  the  magneto¬ 
sphere  regularly  changes  its  shape  (i.e.,  degree  of  distortion)  3nd  not  merely 
its  size,  in  response  to  change  in  solar  wind  conditions.  The  need  r  generate 
magnetospheric  scaling  relations  to  account  for  changes  in  shape  is  thus  once 
again  emphasized. 

E,  Two-Parameter  Scaling  Relations 

The  analysis  described  above  provided  the  motivation  for  proceeding  to  a 

specific  extension  of  the  magnetospheric  scaling  relations.  We  will  here  briefly 

outline  the  progression  this  problem  made  under  this  contract.  Wo  will  illustrate 

that  progress  by  means  of  an  example  which  shows  that  the  new  approach  is  well 

chosen  and  that  it  provides  answers  to  a  new  class  of  problems. 

When  the  set  of  independent  variables  is  restricted  to  the  solar  wind 

stagnation  pressure  (i.e.,  only  the  normal  component  of  the  solar  wind  force 

on  the  magnetopause)  and  the  planetary  magnetic  moment  (including  the  ring 

current  contribution),  there  is  only  one  combination  of  independent  parameters 

that  has  the  dimension  of  length.  The  length  can  be  chosen  to  be  the  planeto- 

centric  distance  to  the  subsolar  stagnation  point,  R  .  It  is  evident  that  if 

m 

a  geometry  is  described  by  a  single  scale  length,  then  change  in  size  but  not 
shape  will  occur  as  R  changes.  This  is  the  underlying  basis  for  the  "constant 

TT1 

shape"  assumption  of  the  existing  scaling  relations. 


Many  previous  studies  support  the  principle  that  the  amount  of  flux 

in  the  tail  is  determined  primarily  by  the  tangential  component  of  the  solar 

wind  force  on  the  magnetopause,  a  neglected  component  in  the  prior  scaling 

iormalism.  Whereas  in  the  prior  formalism  the  independent  variable  could 

be  taken  to  be  the  solar  wind  stagnation  pressure  or  the  scale  size  R  ,  since 

m 

one  is  uniquely  given  in  terms  of  the  other,  in  the  present  formalism  one  may 

choose  a  solar  wind  pair  of  independent  variables,  namely:  the  normal  and 

tangential  drags  on  the  magnetopause,  or  an  internal  pair  of  variables,  that, 

in  principle,  mav  be  derived  from  the  solar  wind  pair.  We  choose  R  and  the 
'  m 

magnetic  fLux  in  the  tail,  F  ,  to  be  the  pair  of  independent  variables.  This 
choice  will  serve  to  extend  the  previous  one-parameter  set  of  scaling  relations 
to  a  two-parameter  set.  The  method  proceeds  as  before,  only  now  there  are  more 
parameters  out  of  which  to  construct  the  dimensions  of  the  magnetospher ic 
quantities  of  interest. 

A  detailed  dimensional  argument,  the  final  realization  of  which  permits 

the  easy  application  of  the  method,  shows  that  the  new  scaling  relations  take 

the  same  form  as  the  old  ones,  except  that  the  old  scaling  constants  become 

functions  of  the  dimensionless  ratio  F„R  /u  M  ,  where  M  is  the  effective 

T  m  o  p  p 

magnetic  moment  of  the  planet. 

To  illustrate  the  result,  consider  the  magnetospher ic  quantity  B  , 

which  is  the  magnetic  field  strength  inside  the  magnetopause  at  the  subsolar 

stagnation  point.  in  the  old  formalism  R  was  fixed  simply  by  B?  /2u  =  p  , 

st  '  J  3  st  o  1 st 

the  stagnation  pressure.  Or  in  terms  of  R  ,  it  was  given  bv  B  =  a  /R  \ 

m  st  o  m 

where  a  is  known  from  theoretical  calculations  performed  in  the  early  1960's 

by  Mead  and  Beard  and  by  Midgley  and  Davis.  In  the  new  formalism  B  must  be 

a ( F_R  /u  M  )  St 

T  m  o  p  ,  /  \  r  , 

- R~3 -  '  where  n(x)  is  some  as  yet  unspecified 

m 


expressed  in  the  form  B 
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function  of  x.  Now,  typical  values  of  the  independent  parameters  for  the 

case  of  earth  give  F  R  /y  MgSs:  Thus,  in  the  spirit  of  going  from  a  prior 

formulation  characterized,  in  effect,  by  the  condition  F„R  /p  Mr,  =  0  to  one 

T  m  o  E 

for  which  the  ratio  is  non-zero  but  small,  we  are  justified  in  claiming 
improvement  if  we  merely  make  a  Taylor  series  expansion  of  the  function  a(x) 
and  terminate  the  series  at  the  linear  term.  That  is,  we  take 


a 


a  +  a, 
o  3 


V\ 


where  a  is  the  old  constant  and  a,  is  a  constant  that  must  be  determined 
o  1 

by  theory  or  experiment,  but  once  it  is  fixed  it  applies  to  all  values  of 
the  independent  parameters.  From  available  observations,  we  estimate  th'  '’ew 
term  to  be  about  30  percent  of  the  old  one,  and  with  opposite  sign.  For 
example , 


3  3  F  R 

“st  ■  inr  +  ITEif*  -  70  Y  -  20  „  -  50  Y 
m  o  e  m 

One  of  the  goals  of  the  present  investigation  is  to  improve  the  determination 
of  the  value  of  the  new  constant,  both  by  appropriately  tailored  data  analysis 
and  by  specific  model  calculations.  In  the  latter,  the  computer  model  of  Olson 
and  Pfitzer  will  be  especially  valuable. 

The  result  just  obtained  can  be  used  to  answer  one  of  the  new  class  of 
questions  involving  changes  in  magnetospheric  shape:  for  constant  solar  wind 
stagnation  pressure  (hence  Bg{_  =  const.),  in  a  flux  transfer  event,  by  how 
much  does  R  change  for  unit  flux  transferred  to  the  tail?  Mathematically 
the  answer  is 


3R 

n 


a,  R 


1 

1  m 


B  =const 
st 


3\J  M  a  +  2a.  F  K 
o  E  o  1  1  m 
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For  the  estimated  value  of  and  typical  values  of  the  other  parameters, 
this  equation  predicts 


AR 


8.5 


This  result  is  very  different  from  that  obtaired  by  merely  equating  the 

flux  missing  from  the  dayside  as  a  result  of  a  decrease  in  R  to  the  increase 

m 

in  the  flux  in  the  tail.  The  proportionality  factor  derived  based  on  this 

purely  kinematical  principle  is  less  than  2.  The  dynamical  forces  that  are 

in  fact  determining  the  constitution  of  the  tail  require  a  much  larger  increase 

in  tail  flux  for  a  given  decrease  in  R  .  The  additional  flux  comes  from  the 

m 

nightside  of  the  magnetosphere  and  is  associated  with  the  observable  phenomenon 

of  the  earthward  movement  of  the  plasma  sheet  when  the  tail  flux  increases. 

This  aspect  of  the  total  tail  phenomenology  was  explicitly  incorporated  in 

the  model  of  Coroniti  and  Kennel  (1972).  kTe  have  here  tied  down  the  many  model 

parameters  that  result  if  one  tries  to  construct  the  various  pieces  of  the 

dynamical  picture,  independently  and  then  assemble  them. 

AR 

The  observed  typical  range  in  -----  that  results  from  flux  transfer  has 

R 

m 

been  shown  hv  Mac z aw a  (1974)  and  by  Holzer  and  Slavin  (1978)  to  be  0.1.  The 
above  formula  then  predicts  a  typical  range  for  AF  /F  of  0.85,  which  corresponds 
to  a  polar  cap  angular  radius  range  from  18°  to  25°,  as  is  observed. 

The  method  of  the  two-parameter  scaling  analysis  is  now  ready  for  general 
exploitation  to  address  questions  concerning  the  energy  input  rate  into  the 
magnetosphere,  and  for  replacing  the  old  ad  hoc  electric  field  scaling  relations 
by  equations  properly  derived  according  to  the  principles  described  above.  The 
result  will  be  a  set  of  relations  connecting  mngnetospher ic  features  and  processes 
to  the  independent  solar  wind  parameters  that  determine  the  normal  and  tangential 
drags  on  the  magnetopause. 
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